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Summary. Various problems of geometry, topology and dynamical systems on sur- 
faces as well as some questions concerning one-dimensional dynamical systems lead 
to the study of closed surfaces endowed with a flat metric with several cone-type 
singularities. Such flat surfaces are naturally organized into families which appear 
to be isomorphic to moduli spaces of holomorphic one-forms. 

One can obtain much information about the geometry and dynamics of an in- 
dividual flat surface by studying both its orbit under the Teichmiiller geodesic flow 
and under the linear group action. In particular, the Teichmiiller geodesic flow plays 
the role of a time acceleration machine (renormalization procedure) which allows to 
study the asymptotic behavior of interval exchange transformations and of surface 
foliations. 

This survey is an attempt to present some selected ideas, concepts and facts in 
Teichmiiller dynamics in a playful way. 
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1 Introduction 

These notes correspond to lectures given first at Les Houches and later, in an 
extended version, at ICTP (Trieste). As a result they keep all blemishes of 
oral presentations. I rush to announce important theorems as facts, and then 
I deduce from them numerous corollaries (which in reality are used to prove 
these very keystone theorems). I omit proofs or replace them by conceptual 
ideas hiding under the carpet all technicalities (which sometimes constitute 
the main value of the proof). Even in the choice of the subjects I poach the 
most fascinating issues, ignoring those which are difficult to present no matter 
how important the latter ones are. These notes also contain some philosophical 
discussions and hopes which some emotional speakers like me include in their 
talks and which one, normally, never dares to put into a written text. 

I am telling all this to warn the reader that this playful survey of some 
selected ideas, concepts and facts in this area cannot replace any serious in- 
troduction in the subject and should be taken with reservation. 

As a much more serious accessible introduction I can recommend a collec- 
tion of introductory surveys of A. Eskin [E], G. Forni [For2], P. Hubert and 
T. Schmidt [HuSdt5] and H. Masur [Ma7], organized as a chapter of the Hand- 
book of Dynamical Systems. I also recommend recent surveys of H. Masur and 
S. Tabachnikov [MaT] and of J. Smillie [S]. The part concerning renormal- 
ization and interval exchange transformations is presented in the article of 
J.-C. Yoccoz [Y] of the current volume in a much more responsible way than 
my introductory exposition in Sec. 5. 
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1.1 Flat Surfaces 

There is a common prejudice which makes us implicitly associate a metric of 
constant positive curvature to a sphere, a metric of constant zero curvature 
to a torus, and a metric of constant negative curvature to a surface of higher 
genus. Actually, any surface can be endowed with a flat metric, no matter 
what the genus of this surface is... with the only reservation that this flat 
metric will have several singular points. Imagine that our surface is made 
from plastic. Then we can flatten it from the sides pushing all curvature to 
some small domains; making these domains smaller and smaller we can finally 
concentrate all curvature at several points. 

Consider the surface of a cube. It gives an example of a perfectly flat sphere 
with eight conical singularities corresponding to eight vertices of the cube. 
Note that our metric is nonsingular on edges: taking a small neighborhood of 
an interior point of an edge and unfolding it we get a domain in a Euclidean 
plane, see Fig. 1. The illusion of degeneration of the metric on the edges comes 
from the singularity of the embedding of our flat sphere into the Euclidean 
space R 3 . 




Fig. 1. The surface of the cube represents a flat sphere with eight conical singular- 
ities. The metric does not have singularities on the edges. After parallel transport 
around a conical singularity a vector comes back pointing to a direction different 
from the initial one, so this flat metric has nontrivial holonomy. 

However, the vertices of the cube correspond to actual conical singularities 
of the metric. Taking a small neighborhood of a vertex we see that it is iso- 
metric to a neighborhood of the vertex of a cone. A flat cone is characterized 
by the cone angle: we can cut the cone along a straight ray with an origin at 
the vertex of the cone, place the resulting flat pattern in the Euclidean plane 
and measure the angle between the boundaries, see Fig. 1. Say, any vertex of 
the cube has cone angle 37r/2 which is easy to see since there are three squares 
adjacent to any vertex, so a neighborhood of a vertex is glued from three right 
angles. 
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Having a manifold (which is in our case just a surface) endowed with a 
metric it is quite natural to study geodesies, which in a fiat metric are locally 
isometric to straight lines. 

General Problem. Describe the behavior of a generic geodesic on a flat sur- 
face. Prove (or disprove) that the geodesic flow is ergodic 1 on a typical (in any 
reasonable sense) flat surface. 

Does any (almost any) fiat surface has at least one closed geodesic which 
does not pass through singular points? 

If yes, are there many closed geodesies like that? Namely, find the asymp- 
totics for the number of closed geodesies shorter than L as the bound L goes 
to infinity. 

Believe it or not there has been no (even partial) advance in solving this 
problem. The problem remains open even in the simplest case, when a sur- 
face is a sphere with only three conical singularities; in particular, it is not 
known, whether any (or even almost any) such flat sphere has at least one 
closed geodesic. Note that in this particular case, when a flat surface is a flat 
sphere with three conical singularities the problem is a reformulation of the 
corresponding billiard problem which we shall discuss in Sect. 2.1. 

1.2 Very Flat Surfaces 

A general flat surface with conical singularities much more resembles a general 
Riemannian manifold than a flat torus. The reason is that it has nontrivial 
holonomy. 

Locally a flat surface is isometric to a Euclidean plane which defines a 
parallel transport along paths on the surface with punctured conical points. A 
parallel transport along a path homotopic to a trivial path on this punctured 
surface brings a vector tangent to the surface to itself. However, if the path is 
not homotopic to a trivial one, the resulting vector turns by some angle. Say, 
a parallel transport along a small closed path around a conical singularity 
makes a vector turn exactly by the cone angle, sec Fig. 1. (Exercise: perform 
a parallel transport of a vector around a vertex of a cube.) 

Nontrivial linear holonomy forces a generic geodesic to come back and 
to intersect itself again and again in different directions; geodesies on a flat 
torus (which has trivial linear holonomy) exhibit radically different behavior. 
Having chosen a direction to the North, we can transport it to any other point 
of the torus; the result would not depend on the path. A geodesic on the torus 
emitted in some direction will forever keep going in this direction. It will either 
close up producing a regular closed geodesic, or will never intersect itself. In 
the latter case it will produce a dense irrational winding line on the torus. 

1 In this context "ergodic" means that a typical geodesic will visit any region in 
the phase space and, moreover, that in average it will spend a time proportional to 
the volume of this region; see Appendix A for details. 
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Fortunately, the class of flat surfaces with trivial linear holonomy is not 
reduced to flat tori. Since we cannot advance in the General Problem from 
the previous section, from now on we confine ourselves to the study of these 
very flat surfaces (often called translation surfaces): that is, to closed ori- 
cntable surfaces endowed with a flat metric having a finite number of conical 
singularities and having trivial linear holonomy. 

Triviality of linear holonomy implies, in particular, that all cone angles 
at conical singularities are integer multiples of 2ir. Locally a neighborhood of 
such a conical point looks like a "monkey saddle" , see Fig. 2. 




Fig. 2. A neighborhood of a conical point with a cone angle 6n can be glued from 
six metric half discs 

As a first example of a nontrivial very flat surface consider a regular 
octagon with identified opposite sides. Since identifications of the sides are 
isometrics, we get a well-defined flat metric on the resulting surface. Since 
in our identifications we used only parallel translations (and no rotations), 
we, actually, get a very flat (translation) surface, ft is easy to see (check it!) 
that our gluing rules identify all vertices of the octagon producing a single 
conical singularity. The cone angle at this singularity is equal to the sum of 
the interior angles of the octagon, that is to 6ir. 

Figure 3 is an attempt to convince the reader that the resulting surface has 
genus two. We first identify the vertical sides and the horizontal sides of the 
octagon obtaining a torus with a hole of the form of a square. To simplify the 
drawing we slightly cheat: namely, we consider another torus with a hole of 
the form of a square, but the new square hole is turned by 7r/4 with respect to 
the initial one. Identifying a pair of horizontal sides of the hole by an isometry 
we get a torus with two holes (corresponding to the remaining pair of sides, 
which are still not identified). Finally, isometrically identifying the pair of 
holes we get a surface of genus two. 
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Fig. 3. Gluing a pretzel from a regular octagon 



Convention 1. From now on by a flat surface we mean a closed oriented surface 
with a flat metric having a finite number of conical singularities, such that the 
metric has trivial linear holonomy. Moreover, we always assume that the flat 
surface is endowed with a distinguished direction; we refer to this direction as 
the "direction to the North" or as the "vertical direction" . 

The convention above implies, in particular, that if we rotate the octagon 
from Fig. 3 (which changes the "direction to the North") and glue a flat 
surface from this rotated octagon, this will give us a different flat surface. 

We make three exceptions to Convention 1 in this paper: billiards in general 
polygons considered at the beginning Sec. 2.1 give rise to flat metrics with 
nontrivial linear holonomy. In Sec. 3.2 we consider flat tori forgetting the 
direction to the North. 

Finally, in Sec. 8.1 we consider half-translation surfaces corresponding to 
flat metrics with holonomy group Z/2Z. Such flat metric is a slight general- 
ization of a very flat metric: a parallel transport along a loop may change the 
direction of a vector, that is a vector v might return as — v after a parallel 
transport. 

1.3 Synopsis and Reader's Guide 

These lectures are an attempt to give some idea of what is known (and what is 
not known) about flat surfaces, and to show what an amazing and marvellous 
object a flat surface is: problems from dynamical systems, from solid state 
physics, from complex analysis, from algebraic geometry, from combinatorics, 
from number theory, ... (the list can be considerably extended) lead to the 
study of flat surfaces. 
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Section 2. Motivations 

To give an idea of how flat surfaces appear in different guises we give some 
motivations in Sec. 2. Namely, we consider billiards in polygons, and, in partic- 
ular, billiards in rational polygons (Sec. 2.1) and show that the consideration 
of billiard trajectories is equivalent to the consideration of geodesies on the 
corresponding flat surface. As another motivation we show in Sec. 2.2 how the 
electron transport on Fermi-surfaces leads to study of foliation defined by a 
closed 1-form on a surface. In Sec. 2.3 we show that under some conditions 
on the closed 1-form such a foliation can be "straightened out" into an appro- 
priate flat metric. Similarly, a Hamiltonian flow defined by the corresponding 
multivalued Hamiltonian on a surface follows geodesies in an appropriate flat 
metric-. 

Section 3. Basic Facts 

A reader who is not interested in motivations can proceed directly to Sec. 3 
which describes the basic facts concerning flat surfaces. For most of applica- 
tions it is important to consider not only an individual flat surface, but an 
entire family of flat surfaces sharing the same topology: genus, number and 
types of conical singularities. In Sec. 3.1 we discuss deformations of flat metric 
inside such families. As a model example we consider in Sec. 3.2 the family 
of flat tori. In Sec. 3.3 we show that a flat structure naturally determines a 
complex structure on the surface and a holomorphic one-form. Reciprocally, 
a holomorphic one-form naturally determines a flat structure. The dictionary 
establishing correspondence between geometric language (in terms of the flat 
metrics) and complex-analytic language (in terms of holomorphic one-forms) 
is very important for the entire presentation; it makes Sec. 3.3 more charged 
than an average one. In Sec. 3.4 we continue establishing correspondence be- 
tween families of flat surfaces and strata of moduli spaces of holomorphic 
one-forms. In Sec. 3.5 we describe the action of the linear group SL(2,M.) on 
flat surfaces - another key issue of this theory. 

We complete Sec. 3 with an attempt to present the following general prin- 
ciple in the study of flat surfaces. In order to get some information about 
an individual flat surface it is often very convenient to find (the closure of) 
the orbit of corresponding element in the family of flat surfaces under the 
action of the group SX(2,R) (or, sometimes, under the action of its diago- 
nal subgroup) . In many cases the structure of this orbit gives comprehensive 
information about the initial flat surface; moreover, this information might 
be not accessible by a direct approach. These ideas are expressed in Sec. 3.6. 
This general principle is illustrated in Sec. 3.7 presenting Masur's criterion 
of unique ergodicity of the directional flow on a flat surface. (A reader not 
familiar with the crgodic theorem can either skip this last section or read an 
elementary presentation of ergodicity in Appendix A.) 
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Section 4- Topological Dynamics of Generic Geodesies 

This section is independent from the others; a reader can pass directly to any 
of the further ones. However, it gives a strong motivation for rcnormalization 
discussed in Sec. 5 and in the lectures by J.-C. Yoccoz [Y] in this volume. It 
can also be used as a formalism for the study of electron transport mentioned 
in Sec. 2.2. 

In Sec. 4.1 we discuss the notion of asymptotic cycle generalizing the ro- 
tation number of an irrational winding line on a torus. It describes how an 
"irrational winding line" on a surface of higher genus winds around a flat 
surface in average. In Sec. 4.2 we heuristically describe the further terms of 
approximation, and we complete with a formulation of the corresponding re- 
sult in Sec. 4.3. 

In fact, this description is equivalent to the description of the deviation of a 
directional flow on a flat surface from the ergodic mean. Sec. 4.3 involves some 
background in ergodic theory; it can be either omitted in the first reading, or 
can be read accompanied by Appendix B presenting the multiplicative ergodic 
theorem. 

Section 5. Renormalization 

This section describes the relation between the Teichmiiller geodesic flow and 
renormalization for interval exchange transformations discussed in the lec- 
tures of J.-C. Yoccoz [Y] in this volume. It is slightly more technical than 
other sections and can be omitted by a reader who is not interested in the 
proof of the Theorem from Sec. 4.3. 

In Sec. 5.1 we show that interval exchange transformations naturally arise 
as the first return map of a directional flow on a flat surface to a transversal 
segment. In Sec. 5.2 we perform an explicit computation of the asymptotic 
cycle (defined in Sec. 4.1) using interval exchange transformations. In Sec. 5.3 
we present a conceptual idea of renormalization, a powerful technique of ac- 
celeration of motion along trajectories of the directional flow. This idea is 
illustrated in Sec. 5.4 in the simplest case where we interpret the Euclidean 
algorithm as a renormalization procedure for rotation of a circle. 

We develop these ideas in Sec. 5.5 describing a concrete geometric renor- 
malization procedure (called Rauzy-Veech induction) applicable to general 
flat surfaces (and general interval exchange transformations). We continue in 
Sec. 5.6 with the elementary formalism of multiplicative cocycles (see also Ap- 
pendix B). Following W. Vccch we describe in Sec. 5.7 zippered rectangles co- 
ordinates in a family of flat surfaces and describe the action of the Teichmiiller 
geodesic flow in a fundamental domain in these coordinates. We show that the 
first return map of the Teichmiiller geodesic flow to the boundary of the fun- 
damental domain corresponds to the Rauzy-Veech induction. In Sec. 5.8 we 
present a short overview of recent results of G. Forni, M. Kontsevich, A. Avila 
and M. Viana concerning the spectrum of Lyapunov exponents of the corre- 
sponding cocycle (completing the proof of the Theorem from Sec. 4.3). As 
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an application of the technique developed in Sec. 5 we show in Sec. 5.9 that 
in the simplest case of tori it gives the well-known encoding of a continued 
fraction by a cutting sequence of a geodesic on the upper half-plane. 

Section 6. Closed geodesies 

This section is basically independent of other sections; it describes the rela- 
tion between closed geodesies on individual flat surfaces and "cusps" on the 
corresponding moduli spaces. It might be useful for those who are interested 
in the global structure of the moduli spaces. 

Following A. Eskin and H. Masur we formalize in Sec. 6.1 the counting 
problems for closed geodesies and for saddle connections of bounded length 
on an individual flat surface. In Sec. 6.2 we present the Siegel-Veech Formula 
and explain a relation between the counting problem and evaluation of the 
volume of a tubular neighborhood of a "cusp" in the corresponding moduli 
space. In Sec. 6.3 we describe the structure of a simplest cusp. We describe the 
structure of general "cusps" (the structure of principal boundary of the moduli 
space) in Sec. 6.4. As an illustration of possible applications we consider in 
Sec. 6.5 billiards in rectangular polygons. 

Section 7 Volume of the Moduli Space 

In Sec. 7.1 we consider very special flat surfaces, so called square-tiled surfaces 
which play a role of integer points in the moduli space. In Sec. 7.2 we present 
the technique of A. Eskin and A. Okounkov who have found an asymptotic 
formula for the number of square-tiled surfaces glued from a bounded number 
of squares and applied these results to evaluation of volumes of moduli spaces. 

As usual, Sec. 7 is independent of others; however, the notion of a square- 
tiled surface appears later in the discussion of Veech surfaces in Sec. 9.5-9.8. 

Section 8. Crash Course in Teichmuller Theory 

We proceed in Sec. 8 with a very brief overview of some elementary background 
in Teichmiillcr theory. Namely, we discuss in Sec. 8.1 the extremal quasicon- 
formal map and formulate the Teichmuller theorem, which we use in Sec. 8.2 
we to define the distance between complex structures (Teichmuller metric). 
We finally explain why the action of the diagonal subgroup in SL(2; R) on the 
space of flat surfaces should be interpreted as the Teichmuller geodesic flow. 

Section 9. Main Conjecture and Recent Results 

In this last section we discuss one of the central problems in the area - a conjec- 
tural structure of all orbits of GL + (2, R). The main hope is that the closure of 
any such orbit is a nice complex subvariety, and that in this sense the moduli 
spaces of holomorphic 1-forms and the moduli spaces of quadratic differentials 
resemble homogeneous spaces under an action of a unipotent group. In this 
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section we also present a brief survey of some very recent results related to 
this conjecture obtained by K. Calta, C. McMullen and others. 

We start in Sec. 9.1 with a geometric description of the GL + (2, R)-action 
in Sec. 9.1 and show why the projections of the orbits (so-called Teichmiiller 
discs) to the moduli space Ai g of complex structures should be considered as 
complex geodesies. In Sec. 9.2 we present some results telling that analogous 
"complex geodesies" in a homogeneous space have a very nice behavior. It is 
known that the moduli spaces are not homogeneous spaces. Nevertheless, in 
Sec. 9.3 we announce one of the main hopes in this field telling that in the 
context of the closures of "complex geodesies" the moduli spaces behave as if 
they were. 

We continue with a discussion of two extremal examples of GL + (2,M)- 
invariant submanifolds. In Sec. 9.4 we describe the "largest" ones: the con- 
nected components of the strata. In Sec. 9.5 we consider flat surfaces S (called 
Veech surfaces) with the "smallest" possible orbits: the ones which are closed. 
Since recently the list of known Veech surfaces was very short. However, 
K. Calta and C. McMullen have discovered an infinite family of Veech surfaces 
in genus two and have classified them. Developing these results C. McMullen 
has proved the main conjecture in genus two. These results of K. Calta and 
C. McMullen are discussed in Sec. 9.7. Finally, we consider in Sec. 9.8 the 
classification of Teichmiiller discs in H.(2) due to P. Hubert, S. Lelievre and 
to C. McMullen. 

Section 10. Open Problems 

In this section we collect open problems dispersed through the text. 
Appendix A. Ergodic Theorem 

In Appendix A we suggest a two-pages exposition of some key facts and con- 
structions in ergodic theory. 

Appendix B. Multiplicative Ergodic Theorem 

Finally, in Appendix B we discuss the Multiplicative Ergodic Theorem which 
is mentioned in Sec. 4 and used in Sec. 5. 

We start with some elementary linear-algebraic motivations in Sec. B.l 
which we apply in Sec. B.2 to the simplest case of a "linear" map of a multi- 
dimensional torus. This examples give us intuition necessary to formulate in 
Sec. B.3 the multiplicative ergodic theorem. Morally, we associate to an ergodic 
dynamical system a matrix of mean differential (or of mean monodromy in 
some cases). We complete this section with a discussion of some basic proper- 
ties of Lyapunov exponents playing a role of logarithms of eigenvalues of the 
"mean differential" ( "mean monodromy" ) of the dynamical system. 
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2 Eclectic Motivations 

In this section we show how flat surfaces appear in different guises: we con- 
sider billiards in polygons, and, in particular, billiards in rational polygons. In 
Sec. 2.1 we show that consideration of billiard trajectories is equivalent to the 
consideration of geodesies on the corresponding flat surface. As another moti- 
vation we show in Sec. 2.2 how the electron transport on Fermi-surfaces leads 
to the study of foliation defined by a closed 1-form on a surface. In Sec. 2.3 we 
show, that under some conditions on the closed 1-form such foliation can be 
"straightened up" in an appropriate flat metric. Similarly, a Hamiltonian flow 
defined by the corresponding multivalued Hamiltonian on a surface follows 
geodesies in an appropriate flat metric. 

2.1 Billiards in Polygons 

Billiards in General Polygons 

Consider a polygonal billiard table and an ideal billiard ball which reflects 
from the walls of the table by the "optical" rule: the angle of incidence equals 
the angle after the reflection. We assume that the mass of our ideal ball is 
concentrated at one point; there is no friction, no spin. 
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We mostly consider regular trajectories, which do not pass through the 
corners of the polygon. However, one can also study trajectories emitted from 
one corner and trapped after several reflections in some other (or the same) 
corner. Such trajectories are called the generalized diagonals. 

To simplify the problem let us start our consideration with billiards in 
triangles. A triangular billiard table is defined by angles a, (3, 7 (proportional 
rescaling of the triangle does not change the dynamics of the billiard). Since 
a + (3 + 7 = 7r the family of triangular billiard tables is described by two real 
parameters. 

It is difficult to believe that the following Problem is open for many 
decades. 

Problem (Billiard in a Polygon). 

1. Describe the behavior of a generic regular billiard trajectory in a generic 
triangle, in particular, prove (or disprove) that the billiard flow is ergodic 2 ; 

2. Does any (almost any) billiard table has at least one regular periodic tra- 
jectory? If the answer is affirmative, does this trajectory survive under 
deformations of the billiard table ? 

3. If a periodic trajectory exists, are there many periodic trajectories like 
that? Namely, find the asymptotics for the number of periodic trajectories 
of length shorter than L as the bound L goes to infinity. 

It is easy to find a special closed regular trajectory in an acute triangle: sec 
the left picture at Fig. 4 presenting the Fagnano trajectory. This periodic tra- 
jectory is known for at least two centuries. However, it is not known whether 
any (or at least almost any) obtuse triangle has a periodic billiard trajectory. 




Fig. 4. Fagnano trajectory and Fox-Kershner construction 



Obtuse triangles with the angles a < (3 < 7 can be parameterized by a 
point of a "simplex" A defined as a + (3 < it/2, a, (3 > 0. For some obtuse tri- 
angles the existence of a regular periodic trajectory is known. Moreover, some 

2 On behalf of the Center for Dynamics and Geometry of Penn State University, 
A. Katok promised a prize of 10.000 euros for a solution of this problem. 
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of these periodic trajectories, called stable periodic trajectories, survive under 
small deformations of the triangle, which proves existence of periodic trajecto- 
ries for some regions in the parameter space A (see the works of G. Galperin, 
A. M. Stepin, Ya. Vorobcts and A. Zcmliakov [GaStVbl], [GaStVb2], [GaZe], 
[Vbl]). It remains to prove that such regions cover the entire parameter space 
A. Currently R. Schwartz is in progress of extensive computer search of sta- 
ble periodic trajectories hoping to cover A with corresponding computer- 
generated regions. 

Now, following Fox and Kershner [FxKr], let us see how billiards in poly- 
gons lead naturally to geodesies on flat surfaces. 

Place two copies of a polygonal billiard table one atop the other. Launch 
a billiard trajectory on one of the copies and let it jump from one copy to 
the other after each reflection (see the right picture at Fig. 4). Identifying the 
boundaries of the two copies of the polygon we get a connected path p on 
the corresponding topological sphere. Projecting this path to any of the two 
"polygonal hemispheres" we get the initial billiard trajectory. 

It remains to note that our topological sphere is endowed with a flat metric 
(coming from the polygon) . Analogously to the flat metric on the surface of a 
cube which is nonsingular on the edges of the cube (see Sec. 1.1 and Fig. 1), the 
flat metric on our sphere is nonsingular on the "equator" obtained from the 
identified boundaries of the two equal "polygonal hemispheres" . Moreover, let 
a; be a point where the path p crosses the "equator" . Unfolding a neighborhood 
of a point x on the "equator" we see that the corresponding fragment of the 
path p unfolds to a straight segment in the flat metric. In other words, the 
path p is a geodesic in the corresponding flat metric. 

The resulting flat metric is not very flat: it has nontrivial linear holonomy. 
The conical singularities of the flat metric correspond to the vertices of the 
polygon; the cone angle of a singularity is twice the angle at the corresponding 
vertex of the polygon. 

We have proved that every geodesic on our flat sphere projects to a billiard 
trajectory and every billiard trajectory lifts to a geodesic. This is why General 
Problem from Sec. 1.1 is so closely related to the Billiard Problem. 

Two Beads on a Rod and Billiard in a Triangle 

It would be unfair not to mention that billiards in polygons attracted a lot of 
attention as (what initially seemed to be) a simple model of a Boltzman gas. 
To give a flavor of this correspondence we consider a system of two elastic 
beads confined to a rod placed between two walls, see Fig. 5. (Up to the best 
of my knowledge this construction originates in lectures of Ya. G. Sinai [Sin].) 

The beads have different masses mi and ni2 they collide between them- 
selves, and also with the walls. Assuming that the size of the beads is negli- 
gible we can describe the configuration space of our system using coordinates 
< x± < X2 < a of the beads, where a is the distance between the walls. 
Rescaling the coordinates as 
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X\ = y/miXi 
X2 = 

we see that the configuration space in the new coordinates is given by a right 
triangle A, see Fig. 5. Consider now a trajectory of our dynamical system. We 
leave to the reader the pleasure to prove the following elementary Lemma: 

Lemma. In coordinates (xi,X2) trajectories of the system of two beads on a 
rod correspond to billiard trajectories in the triangle A. 

Billiards in Rational Polygons 

We have seen that taking two copies of a polygon we can reduce the study of 
a billiard in a general polygon to the study of geodesies on the corresponding 
flat surface. However, the resulting flat surface has nontrivial linear holonomy, 
it is not "very flat" . 

Nevertheless, a more restricted class of billiards, namely, billiards in ratio- 
nal polygons, lead to "very flat" (translation) surfaces. 

A polygon is called rational if all its angles are rational multiples of tt. 
A billiard trajectory emitted in some direction will change direction after the 
first reflection, then will change direction once more after the second reflection, 
etc. However, for any given billiard trajectory in a rational billiard the set 
of possible directions is finite, which make billiards in rational polygons so 
different from general ones. 

As a basic example consider a billiard in a rectangle. In this case a generic 
trajectory at any moment goes in one of four possible directions. Developing 
the idea with the general polygon we can take four (instead of two) copies 
of our billiard table (one copy for each direction). As soon as our trajectory 
hits the wall and changes the direction we make it jump to the corresponding 
copy of the billiard - the one representing the corresponding direction. 

By construction each side of every copy of the billiard is identified with 
exactly one side of another copy of the billiard. Upon these identifications the 
four copies of the billiard produce a closed surface and the unfolded billiard 
trajectory produces a connected line on this surface. We suggest to the reader 
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to check that the resulting surface is a torus and the unfolded trajectory is a 
geodesic on this fiat torus, see Fig. 6. 




Fig. 6. Billiard in a rectangle corresponds to directional flow on a flat torus. Billiard 
in a right triangle (n/8, 3tt/8, tt/2) leads to directional flow on a flat surface obtained 
from the regular octagon. 

A similar unfolding construction (often called Katok-Zemliakov construc- 
tion) works for a billiard in any rational polygon. Say, for a billiard in a right 
triangle with angles (tt/8, 3tt/8, 7t/2) one has to take 16 copies (corresponding 
to 16 possible directions of a given billiard trajectory). Appropriate identifi- 
cations of these 16 copies produce a regular octagon with identified opposite 
sides (see Fig. 6). We know from Sec. 1.2 and from Fig. 3 that the correspond- 
ing flat surface is a "very flat" surface of genus two having a single conical 
singularity with the cone angle 67r. 

Exercise. What is the genus of the surface obtained by Katok-Zcmliakov 
construction from an isosceles triangle (37r/8, 3tt/8, 7r/4)? How many coni- 
cal points does it have? What are the cone angles at these points? Hint: this 
surface can not be glued from a regular octagon. 

It is quite common to unfold a rational billiard in two steps. We first unfold 
the billiard table to a polygon, and then identify the appropriate pairs of sides 
of the resulting polygon. Note that the polygon obtained in this intermediate 
step is not canonical. 

Show that a generic billiard trajectory in the right triangle with angles 
(7r/2, 7t/5, 37r/10) has 20 directions. Show that both polygons at Fig. 7 can 
be obtained by Katok-Zemliakov construction from 20 copies of this triangle. 
Verify that after identification of parallel sides of these polygons we obtain 
isometric very flat surfaces (see also [HuSdt5]). What genus, and what conical 
points do they have? What are the cone angles at these points? 
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Fig. 7. We can unfold the billiard in the right triangle (tt/2, tt/5, 37r/10) into 
different polygons. However, the resulting very flat surfaces are the same (see 
also [HuSdto]). 

Note that in comparison with the initial construction, where we had only 
two copies of the billiard table we get a more complicated surface. However, 
what we gain is that in this new construction our flat surface is actually "very 
flat" : it has trivial linear holonomy. It has a lot of consequences; say, due to a 
Theorem of H. Masur [Ma4] it is possible to find a regular periodic geodesic on 
any "very flat" surface. If the flat surface was constructed from a billiard, the 
corresponding closed geodesic projects to a regular periodic trajectory of the 
corresponding billiard which solves part of the Billiard Problem for billiards 
in rational polygons. 

We did not intend to present in this section any comprehensive information 
about billiards, our goal was just to give a motivation for the study of flat 
surfaces. A reader interested in billiards can get a good idea on the subject 
from a very nice book of S. Tabachnikov [T]. Details about billiards in polygons 
(especially rational polygons) can be found in the surveys of E. Gutkin [Gul], 
P. Hubert and T. Schmidt [HuSdt5], H. Masur and S. Tabachnikov [MaT] and 
J. Smillic [S]. 

2.2 Electron Transport on Fermi- Surfaces 

Consider a periodic surface M in R 3 (i.e. a surface invariant under trans- 
lations by any integer vector in Z 3 ). Such a surface can be constructed in a 
fundamental domain of a cubic lattice, see Fig. 8, and then reproduced re- 
peatedly in the lattice. Choose now an aflme plane in R 3 and consider an 
intersection line of the surface by the plane. This intersection line might have 
some closed components and it may also have some unbounded components. 
The question is how does an unbounded component propagate in R 3 ? 

The study of this subject was suggested by S. P. Novikov about 1980 
(see [N]) as a mathematical formulation of the corresponding problem concern- 
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Fig. 8. Riemann surface of genus 3 embedded into a torus T 3 

ing electron transport in metals. A periodic surface represents a Fermi- surface, 
affinc plane is a plane orthogonal to a magnetic field, and the intersection line 
is a trajectory of an electron in the so-called inverse lattice. 




Fig. 9. Fermi surfaces of tin, iron and gold correspond to Riemann surfaces of high 
genera. (Reproduced from [AzLK] which cites [AGLP] and [WYa] as the source) 

It was known since extensive experimental research in the 50s and 60s 
that Fermi-surfaces may have fairly complicated shape, see Fig. 9; it was also 
known that open (i.e. unbounded) trajectories exist, see Fig. 10, however, up 
to the beginning of the 80s there were no general results in this area. 

In particular, it was not known whether open trajectories follow (in a large 
scale) the same direction, whether there might be some scattering (trajectory 
comes from infinity in one direction and then after some scattering goes to 
infinity in some other direction, whether the trajectories may even exhibit 
some chaotic behavior? 

Let us see now how this problem is related to flat surfaces. 

First note that passing to a quotient M 3 /Z 3 = T 3 we get a closed orientable 
surface M 2 C T 3 from the initial periodic surface M 2 . Say, identifying the 
opposite sides of a unit cube at Fig. 8 we get a closed surface M 2 of genus 
.9 = 3. 
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Fig. 10. Stereographic projection of the magnetic field directions (shaded regions 
and continuous curves) which give rise to open trajectories for some Fermi-surfaces 
(experimental results in [AzLK]). 

We are interested in plane sections of the initial periodic surface M 2 . This 
plane sections can be viewed as level curves of a linear function f(x,y, z) = 
ax + by + cz restricted to M 2 . 

Consider now a closed differential 1-form Co = a dx + bdy + cdz in R 
and its restriction to M 2 . A closed 1-form defines a codimension-one foliation 
on a manifold: locally one can represent a closed one-form as a differential 
of a function, Co = df; the foliation is defined by the levels of the function 
/. We prefer to use the 1-form Co = adx + bdy + cdz to the linear function 
f(x, y, z) = ax + by + cz because wc cannot push the function f(x, y, z) into a 
torus T 3 while the 1-form a-> = a dx+b dy+ cdz is well-defined inT 3 . Moreover, 
after passing to a quotient over the lattice R 3 — > E 3 /Z 3 the plane sections of 
M 2 project to the leaves of the foliation defined by restriction of the closed 
1-form u> in T 3 to the surface M 2 . 

Thus, our initial problem can be reformulated as follows. 

Problem (Novikov's Problem on Electron Transport). Consider a fo- 
liation defined by a linear closed 1-form to = adx + bdy + cdz on a closed 
surface M 2 C T 3 embedded into a three-dimensional torus. How do the leaves 
of this foliation get unfolded, when we unfold the torus T 3 to its universal 
cover R 3 ? 

The foliation defined by a closed 1-form on a surface is a subject of discus- 
sion of the next section. We shall see that under some natural conditions such 
a foliation can be "straightened up" to a geodesic foliation in an appropriate 
flat metric. 

The way in which a geodesic on a flat surface gets unfolded in the universal 
Abelian cover is discussed in detail in Sec. 4. 

To be honest, we should admit that 1-forms as in the Problem above usu- 
ally do not satisfy these conditions. However, a surface as in this Problem can 
be decomposed into several components, which (after some surgery) already 
satisfy the necessary requirements. 



20 Anton Zorich 



References for Details 

There is a lot of progress in this area, basically due to S. Novikov's school, and 
especially to I. Dynnikov, and in many cases Novikov's Problem on electron 
transport is solved. 

In [Zol] the author proved that for a given Fermi-surface and for an open 
dense set of directions of planes any open trajectory is bounded by a pair of 
parallel lines inside the corresponding plane. 

In a series of papers I. Dynnikov applied a different approach: he fixed the 
direction of the plane and deformed a Fermi-surface inside a family of level 
surfaces of a periodic function in R 3 . He proved that for all but at most one 
level any open trajectory is also bounded between two lines. 

However, I. Dynnikov has constructed a series of highly elaborated exam- 
ples showing that in some cases an open trajectory can "fill" the plane. In 
particular, the following question is still open. Consider the set of directions 
of those hypcrplanes which give "nontypical" open trajectories. Is it true that 
this set has measure zero in the space KP 2 of all possible directions? What 
can be said about Hausdorff dimension of this set? 

For more details we address the reader to papers [Dl], [D2] and [NM]. 

2.3 Flows on Surfaces and Surface Foliations 

Consider a closed 1-form on a closed oricntablc surface. Locally a closed 1- 
form to can be represented as the differential of a function ui = df . The level 
curves of the function / locally define the leaves of the closed 1-form u. (The 
fact that the function / is defined only up to a constant does not affect the 
structure of the level curves.) We get a foliation on a surface. 

In this section we present a necessary and sufficient condition which tells 
when one can find an appropriate flat metric such that the foliation defined 
by the closed 1-form becomes a foliation of parallel geodesies going in some 
fixed direction on the surface. This criterion was given in different context by 
different authors: [Clb], [Katl], [HbMa]. We present here one more formulation 
of the criterion. Morally, it says that the foliation defined by a closed 1-form 
lo can be "straightened up" in an appropriate flat metric if and only if the 
form ui does not have closed leaves homologous to zero. In the remaining part 
of this section we present a rigorous formulation of this statement. 

Note that a closed 1-form w on a closed surface necessarily has some critical 
points: the points where the function / serving as a "local antiderivative" 
lo = df has critical points. 

The first obstruction for "straightening" is the presence of minima and 
maxima: such critical points should be forbidden. Suppose now that the closed 
1-form has only isolated critical points and all of them are "saddles" (i.e. to 
does not have minima and maxima). Say, a form defined in local coordinates 
as df, where / = x 3 +y 3 has a saddle point in the origin (0, 0) of our coordinate 
chart, see Fig. 11. 
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Fig. 11. Horizontal foliation in a neighborhood of a saddle point. Topological 
(nonmetric) picture 

There are several singular leaves of the foliation landing at each saddle; 
say, a saddle point from Fig. 11 has six prongs (separatrices) representing 
the critical leaves. Sometimes a critical leaf emitted from one saddle can land 
to another (or even to the same) saddle point. In this case we say that the 
foliation has a saddle connection. 

Note that the foliation defined by a closed 1-form on an oriented surface 
gets natural orientation (defined by grad(/) and by the orientation of the 
surface) . Now we are ready to present a rigorous formulation of the criterion. 

Theorem. Consider a foliation defined on a closed orientable surface by a 
closed 1-form tu. Assume that to does not have neither minima nor maxima but 
only isolated saddle points. The foliation defined by u> can be represented as a 
geodesic foliation in an appropriate fiat metric if and only if any cycle obtained 
as a union of closed paths following in the positive direction a sequence of 
saddle connections is not homologous to zero. 

In particular, if there are no saddle connections at all (provided there are 
no minima and maxima) it can always be straightened up. (In slightly different 
terms it was proved in [Katl] and in [HbMa].) 

Note that saddle points of the closed 1-form u> correspond to conical points 
of the resulting flat metric. 

One can consider a closed 1-form a; as a multivalued Hamiltonian and con- 
sider corresponding Hamiltonian flow along the leaves of the foliation defined 
by ui. On the torus T 2 is was studied by V. I. Arnold [Ald2] and by K. Khanin 
and Ya. G. Sinai [KhSin]. 

3 Families of Flat Surfaces and Moduli Spaces of 
Abelian Differentials 

In this section we present the generalities on flat surfaces. We start in Sec. 3.1 
with an elementary construction of a flat surface from a polygonal pattern. 
This construction explicitly shows that any flat surface can be deformed in- 
side an appropriate family of flat surfaces. As a model example we consider 
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in Sec. 3.2 the family of flat tori. In Sec. 3.3 we show that a flat structure 
naturally determines a complex structure on the surface and a holomorphic 
one-form. Reciprocally, a holomorphic one-form naturally determines a flat 
structure. The dictionary establishing correspondence between geometric lan- 
guage (in terms of the flat metrics) and complex-analytic language (in terms 
of holomorphic one-forms) is very important for the entire presentation; it 
makes Sec. 3.3 more charged than an average one. In Sec. 3.4 we continue 
with establishing correspondence between families of flat surfaces and strata 
of moduli spaces of holomorphic one-forms. In Sec. 3.5 we describe the action 
of the linear group SL(2, R) on flat surfaces - another key issue of this theory. 

We complete Sec. 3 with an attempt to present the following general prin- 
ciple in the study of flat surfaces. In order to get some information about an 
individual flat surface it is often very convenient to find (the closure of) the 
orbit of the corresponding clement in the family of flat surfaces under the 
action of the group SL(2,R) (or, sometimes, under the action of its diagonal 
subgroup). In many cases the structure of this orbit gives a comprehensive 
information about the initial flat surface; moreover, this information might 
be not accessible by a direct approach. These ideas are expressed in Sec. 3.6. 
This general principle is illustrated in Sec. 3.7 presenting Masur's criterion of 
unique ergodicity of the directional flow on a flat surface. (A reader not famil- 
iar with ergodic theorem can either skip this last section or read an elementary 
presentation of ergodicity in Appendix A.) 

3.1 Families of Flat Surfaces 

In this section we present a construction which allows to obtain a large variety 
of flat surfaces, and, moreover, allows to continuously deform the resulting flat 
structure. Later on we shall see that this construction is even more general 
than it may seem at the beginning: it allows to get almost all flat surfaces in 
any family of flat surfaces sharing the same geometry (i.e. genus, number and 
types of conical singularities). The construction is strongly motivated by an 
analogous construction in the paper of H. Masur [Ma3]. 

Consider a collection of vectors Vi, . . . , v„ in M 2 and construct from these 
vectors a broken line in a natural way: a j-th edge of the broken line is 
represented by the vector Vj. Construct another broken line starting at the 
same point as the initial one by taking the same vectors but this time in the 
order tV(i), ■ • ■ , ?V(n)j where it is some permutation of n elements. 

By construction the two broken lines share the same endpoints; suppose 
that they bound a polygon as in Fig. 12. Identifying the pairs of sides cor- 
responding to the same vectors «j, j = 1, . . . ,n, by parallel translations we 
obtain a flat surface. 

The polygon in our construction depends continuously on the vectors Vj. 
This means that the topology of the resulting flat surface (its genus g, the 
number m and the types of the resulting conical singularities) do not change 
under small deformations of the vectors Vj. Say, we suggest to the reader to 
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Fig. 12. Identifying corresponding pairs of sides of this polygon by parallel trans- 
lations we obtain a flat surface. 

check that the fiat surface obtained from the polygon presented in Fig. 12 has 
genus two and a single conical singularity with cone angle 6tt. 

3.2 Toy Example: Family of Flat Tori 

In the previous section we have seen that a flat structure can be deformed. 
This allows to consider a flat surface as an element of a family of flat surfaces 
sharing a common geometry (genus, number of conical points). In this section 
we study the simplest example of such family: we study the family of flat tori. 
This time we consider the family of flat surfaces globally. We shall see that 
it has a natural structure of a noncompact complex-analytic manifold (to be 
more honest - orbifold). This "baby family" of flat surfaces, actually, exhibits 
all principal features of any other family of flat surfaces, except that the family 
of flat tori constitutes a homogeneous space endowed with a nice hyperbolic 
metric, while general families of flat surfaces do not have the structure of a 
homogeneous space. 

To simplify consideration of flat tori as much as possible we make two 
exceptions from the usual way in which we consider flat surfaces. Temporarily 
(only in this section) we forget about the choice of the direction to North: in 
this section two isometric flat tori define the same element of the family of all 
flat tori. Another exception concerns normalization. Almost everywhere below 
we consider the area of any flat surface to be normalized to one (which can be 
achieved by a simple homothcty) . In this section it would be more convenient 
for us to apply homothety in the way that the shortest closed geodesic on our 
flat torus would have length 1 . Find the closed geodesic which is next after the 
shortest one in the length spectrum. Measure the angle <f>, where < 4> < ir 
between these two geodesies; measure the length r of the second geodesic and 
mark a point with polar coordinates (r, </>) on the upper half-plane. This point 
encodes our torus. 
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neighborhood 
of a cusp 





Fig. 13. Space of fiat tori 



Reciprocally, any point of the upper half-plane defines a flat torus in the 
following way. (Following the tradition we consider the upper half-plane as 
a complex one.) A point x + iy defines a parallelogram generated by vectors 
vi = (1,0) and V2 = (x + iy), see Fig. 13. Identifying the opposite sides of 
the parallelogram we get a flat torus. 

To make this correspondence bijective we have to be sure that the vector 
vi represents the shortest closed geodesic. This means that the point (x + 
iy) representing V2 cannot be inside the unit disc. The condition that V2 
represents the geodesic which is next after the shortest one in the length 
spectrum implies that —1/2 < x < 1/2, sec Fig. 13. Having mentioned these 
two hints we suggest to the reader to prove the following Lemma. 

Lemma 1. The family of flat tori is parametrized by the shadowed fundamen- 
tal domain from Fig. 13, where the parts of the boundary of the fundamental 
domain symmetric with respect to the vertical axis (0,iy) are identified. 

Note that topologically we obtain a sphere punctured at one point: the 
resulting surface has a cusp. Tori represented by points close to the cusp are 
"disproportional" : they are very narrow and very long. In other words they 
have an abnormally short geodesic. 

Note also that there are two special points on our modular curve: they 
correspond to points with coordinates (0 + i) and ±1/2 + i\fi/2. The corre- 
sponding tori has extra symmetry, they can be represented by a square and by 
a regular hexagon with identified opposite sides correspondingly The surface 
glued from the fundamental domains has "corners" at these two points. 

There is an alternative more algebraic approach to our problem. Actually, 
the fundamental domain constructed above is known as a modular curve, it 
parameterizes the space of lattices of area one (which is isomorphic to the 
space of flat tori). It can be seen as a double quotient 
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3.3 Dictionary of Complex- Analytic Language 

We have seen in Sec. 3.1 how to construct a flat surfaces from a polygon as 
on Fig. 12. 

Note that the polygon is embedded into a complex plane C, where the 
embedding is defined up to a parallel translation. (A rotation of the polygon 
changes the vertical direction and hence, according to Convention 1 it changes 
the corresponding flat surface.) 

Consider the natural coordinate z in the complex plane. In this coordinate 
the parallel translations which we use to identify the sides of the polygon are 
represented as 

z = z + const (1) 

Since this correspondence is holomorphic, it means that our flat surface S with 
punctured conical points inherits the complex structure. It is an exercise in 
complex analysis to check that the complex structure extends to the punctured 
points. 

Consider now a holomorphic 1-form dz in the initial complex plane. When 
we pass to the surface S the coordinate z is not globally defined anymore. 
However, since the changes of local coordinates are defined by the rule (1) we 
see that dz = dz' . Thus, the holomorphic 1-form dz on C defines a holomorphic 
1-form lj on S which in local coordinates has the form u> = dz. Another 
exercise in complex analysis shows that the form o> has zeroes exactly at 
those points of S where the flat structure has conical singularities. 

In an appropriate local coordinate w in a neighborhood of zero (different 
from the initial local coordinate z) a holomorphic 1-form can be represented 
as w d dw, where d is called the degree of zero. The form ui has a zero of degree 
d at a conical point with cone angle 2n(d + 1). 

Recall the formula for the sum of degrees of zeroes of a holomorphic 1-form 
on a Ricmann surface of genus g: 

ni 

Y,d j = 2g-2 

3 = 1 

This relation can be interpreted as the formula of Gauss-Bonnet for the flat 
metric-. 

Vectors Vj representing the sides of the polygon can be considered as 
complex numbers. Let Vj be joining vertices Pj and Pj+\ of the polygon. 
Denote by pj the resulting path on S joining the points Pj,Pj +i e S. Our 
interpretation of Vj as of a complex number implies the following obvious 
relation: 

Vj : = / dz= uj (2) 
J Pi J Pi 

Note that the path pj represents a relative cycle: an element of the relative 
homology group Hi(S, {Pi, . . . , P m }', Z) of the surface S relative to the finite 
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collection of conical points {Pi, . . . , P m }. Relation (2) means that Vj repre- 
sents a period of u>: an integral of u> over a relative cycle pj. 

Note also that the flat area of the surface S equals the area of the original 
polygon, which can be measured as an integral of dx A dy over the polygon. 
Since in the complex coordinate z we have dx A dy = jdz A dz we get the 
following formula for the flat i'll'Oi'l of S: 

■ p ■ a 

area(S) = -Ju>A(D = - ^{AjBj - A 3 B 3 ) (3) 

s j—i 

Here we also used the Riemann bilinear relation which expresses the integral 
j s ui A ui in terms of absolute periods Aj, Bj of to, where the absolute periods 
Aj, Bj are the integrals of w with respect to some symplcctic basis of cycles. 

An individual flat surface defines a pair: (complex structure, holomorphic 
1-form). A family of flat surfaces (where the flat surfaces are as usual endowed 
with a choice of the vertical direction) corresponds to a stratum H(di, . . . , d m ) 
in the moduli space of holomorphic 1-forms. Points of the stratum are rep- 
resented by pairs (point in the moduli space of complex structures, holomor- 
phic 1-form in the corresponding complex structure having zeroes of degrees 
di, ■ ■ ■ , dm). 

The notion "stratum" has the following origin. The moduli space of pairs 
(holomorphic 1-form, complex structure) forms a natural vector bundle over 
the moduli space M g of complex structures. A fiber of this vector bundle 
is a vector space C 9 of holomorphic 1-forms in a given complex structure. 
We already mentioned that the sum of degrees of zeroes of a holomorphic 
1-form on a Riemann surface of genus g equals 2g — 2. Thus, the total space 
Tt g of our vector bundle is stratified by subspaces of those forms which have 
zeroes of degrees exactly d\, . . . , d m , where d\ + ■ ■ ■ + d rn = 2g — 2. Say, for 
g = 2 we have only two partitions of number 2, so we get two strata H(2) 
and H(l,l). For g = 3 we have five partitions, and correspondingly five strata 
W(4),W(3, 1),W(2,2),W(2, 1, 1),W(1, 1, 1, 1). 

H(di, . . . , d m ) C H g 

i (4) 

Mg 

Every stratum TL{d\, . . . , d m ) is a complex-analytic orbifold of dimension 

dime H(di, . . . , d m ) =2g + m—l (5) 

Note, that an individual stratum 7i{d\, . . . , d m ) does not form a fiber bundle 
over Mg- For example, according to our formula, dime 7~L(2g — 2) = 2g, while 
dime M g = 3g — 3. 

We showed how the geometric structures related to a flat surface define 
their complex-analytic counterparts. Actually, this correspondence goes in two 
directions. We suggest to the reader to make the inverse translation: to start 
with complex-analytic structure and to see how it defines the geometric one. 
This correspondence can be summarized in the following dictionary. 
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Table 1. Correspondence of geometric and complex-analytic notions 



Geometric language 


Complex-analytic language 


flat structure (including a choice 
of the vertical direction) 


complex structure + a choice 
of a holomorphic f-form to 


conical point 
with a cone angle 2-7r(d + f ) 


zero of degree d 
of the holomorphic 1-form lo 
(in local coordinates lo = w d dw) 


side Vj of a polygon 


relative period Jp J+1 w = f v , lo 
of the 1-form lo 


area of the flat surface S 


\ J s u A O = 


family of flat surfaces sharing the same 
types 27r(di + 1), . . . , 2Tv(d m + 1) 
of cone angles 


stratum Ti(di , . . . , d m ) in the 
moduli space of Abelian differentials 


coordinates in the family: 
vectors Vj 
defining the polygon 


coordinates in H(di, . . . , d m ) : 
collection of relative periods of lo, 
i.e. cohomology class 
M e H 1 (S,{P 1 ,...,Pm};C) 



3.4 Volume Element in the Moduli Space of Holomorphic 
One-Forms 

In the previous section we have considered vectors vi,...,v„ determining 
the polygon from which we glue a flat surface 5 as on Fig. 12. We have 
identified these vectors Vj € R 2 ~ C with complex numbers and claimed 
(without proof) that under this identification Vi, . . . ,v„ provide us with lo- 
cal coordinates in the corresponding family of flat surfaces. We identify ev- 
ery such family with a stratum TL{d\, . . . , d m ) in the moduli space of holo- 
morphic 1-forms. In complex-analytic language we have locally identified a 
neighborhood of a "point" (complex structure, holomorphic 1-form lo) in 
the corresponding stratum with a neighborhood of the cohomology class 
[lo] g H 1 (S,{P u ...,P m };C). 

Note that the cohomology space H 1 (S, {Pi, . . . , P m }', C) contains a natu- 
ral integer lattice -ff 1 (S l , {Pi, . . . , P m }', Z © \/^TZ). Consider a linear volume 
element dv in the vector space H l (S, {Pi, . . . , Pn}; C) normalized in such a 
way that the volume of the fundamental domain in the "cubic" lattice 

P x (5, {P, . . . , P m }; Z © V^TZ) c H^S, {Pi, . . . , P m }; C) 



is equal to one. In other terms 
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1 1 



dvidvi . . . <7v„<7v 



J (2v^T) 



1 1 



where J is the determinant of a change of the basis {vi, . . . , v„} considered 
as a basis in the first relative homology to some "symplectic" basis in the first 
relative homology. 

Consider now the real hypersurface 



defined by the equation area(S) = 1. Taking into consideration formula (3) 
for the function area(S) we see that the hypersurface Hi(dx, ■ ■ . , d m ) defined 
as area(S) = 1 can be interpreted as a "unit hyperboloid" defined in local 
coordinates as a level of the indefinite quadratic form (3) . 

The volume element dv can be naturally restricted to a hypcrplane defined 
as a level hypersurface of a function. We denote the corresponding volume 
element on H.i(di, . . . , dm) by dv\. 

Theorem (H. Masur. W. A. Veech). The total volume 



of every stratum is finite. 

The values of these volumes were computed only recently by A. Eskin 
and A. Okounkov [EOk], twenty years after the Theorem above was proved 
in [Ma3], [Ve3] and [Ve8]. We discuss this computation in Sec. 7. 

3.5 Action of SX(2,R) on the Moduli Space 

In this section we discuss a property of flat surfaces which is, probably, the 
most important in our study: we show that the linear group acts on every 
family of flat surfaces, and, moreover, acts ergodically (see Append. A for 
discussion of the notion of ergodicity). This enables us to apply tools from 
dynamical systems and from ergodic theory. 

Consider a flat surface S and consider a polygonal pattern obtained by 
unwrapping it along some geodesic cuts. For example, one can assume that 
our flat surface S is glued from a polygon 77 C R 2 as on Fig. 12. Consider a 
linear transformation g £ G7 + (2,R) of the plane R 2 . It changes the shape of 
the polygon. However, the sides of the new polygon gll are again arranged 
into pairs, where the sides in each pair are parallel and have equal length 
(different from initial one), see Fig. 14. Thus, identifying the sides in each 
pair by a parallel translation we obtain a new flat surface gS. 

It is easy to check that the surface gS does not depend on the way in 
which S was unwrapped to a polygonal pattern 77. It is clear that all topolog- 
ical characteristics of the new flat surface gS (like genus, number and types 



Hi(di, . . . , 



d m ) c H{dx, . . . , 
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Fig. 14. Action of the linear group on flat surfaces 



of conical singularities) are the same as those of the initial flat surface S. 
Hence, we get a continuous action of the group GL + (2,M.) on each stratum 
H(di, . . . , d m ). 

Considering the subgroup SX(2,R) of area preserving linear transforma- 
tions we get the action of SL(2, R) on the "unit hyperboloid" Hi(di, . . . , d m ). 

(e t \ 

Considering the diagonal subgroup I ^ _ t I C SL(2, R) we get a continuous 

action of this one-parameter subgroup on each stratum r H{d\ 1 . . . ,d m ). This 
action induces a natural flow on the stratum, which is called the Teichmiiller 
geodesic flow. 

Key Theorem (H. Masur. W. A. Veech). The action of the groups 



SL(2,] 



and 



e 



t I preserves the measure dv\. Both actions are ergodic 



with respect to this measure on each connected component of every stratum 
Hi(dx,..., d m ). 

This theorem might seem quite surprising. Consider almost any flat surface 
S as in Fig. 12. "Almost any flat surface" is understood as "corresponding to 
a set of parameters vi, . . . , V4 of full measure; here the vectors v, define the 
polygon 77 from Fig. 12. 

Now start contracting the polygon 77 it in the vertical direction and ex- 
panding it in the horizontal direction with the same coefficient e*. The theorem 
says, in particular, that for an appropriate i £ 1 the deformed polygon will 
produce a flat surface gtS which would be arbitrary close to the flat surface 
So obtained from the regular octagon as on Fig. 3 since a trajectory of al- 
most any point under an ergodic flow is everywhere dense (and even "well 
distributed"). However, it is absolutely clear that acting on our initial poly- 
gon 77 from Fig. 12 with expansion-contraction we never get close to a regular 
octagon... Is there a contradiction?.. 

There is no contradiction since the statement of the theorem concerns fiat 
surfaces and not polygons. In practice this means that we can apply expansion- 
contraction to the polygon 77, which does not change too much the shape of 
the polygon, but radically changes the flat structure. Then we can change 
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the way in which we unwrap the flat surface g t S (see Fig. 15). This radically 
changes the shape of the polygon, but does not change at all the flat structure! 




Fig. 15. The first modification of the polygon changes the flat structure while the 
second one just changes the way in which we unwrap the flat surface 



3.6 General Philosophy 

Now we are ready to describe informally the basic idea of our approach to the 
study of flat surfaces. Of course it is not universal; however, in many cases it 
appears to be surprisingly powerful. 

Suppose that we need some information about geometry or dynamics of an 
individual flat surface S. Consider the element S in the corresponding family of 
flat surfaces H(di, d m ). Denote by C{S) = GL+(2,R)S C H(di,. ...dm) 
the closure of the GL + (2, R)-orbit of S in Ti.(d\, . . . , d m ). In numerous cases 
knowledge about the structure ofC(S) gives a comprehensive information about 
geometry and dynamics of the initial flat surface S. Moreover, some delicate 
numerical characteristics of S can be expressed as averages of simpler char- 
acteristics over C(S). 

The remaining part of this survey is an attempt to show some implementa- 
tions of this general philosophy. The first two illustrations would be presented 
in the next section. 

We have to confess that we do not tell all the truth in the formulation 
above. Actually, there is a hope that this philosophy extends much further. 
A closure of an orbit of an abstract dynamical system might have extremely 
complicated structure. According to the most optimistic hopes, the closure 
C(S) of the GL + (2, R)-orbit of any flat surface S is a nice complex-analytic 
variety. Moreover, according to the most daring conjecture it would be possible 
to classify all these GL + (2, M)-invariant subvarieties. For genus two the latter 
statements were recently proved by C. McMullcn (see [McM2] and [McM3]) 
and partly by K. Calta [Clt]. 
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Wc discuss this hope in more details in Sec. 9, in particular, in Sec. 9.3. 
We complete this section by a Theorem which supports the hope for some 
nice and simple description of orbit closures. 

Theorem (M. Kontsevich) . Suppose that a closure C{S) in 7i(di, . . . , d m ) 
of a GL + (2,R)-orbit of some flat surface S is a complex- analytic subvariety. 
Then in cohomological coordinates i? 1 (S l , {-Pi, ■ ■ ■ , P m }; C) it is represented by 
an afflne subspace. 

3.7 Implementation of General Philosophy 

In this section we present two illustrations showing how the "general philos- 
ophy" works in practice. 

Consider a directional flow on a flat surface S. It is called minimal when the 
closure of any trajectory gives the entire surface. When a directional flow on a 
flat torus is minimal, it is necessarily ergodic, in particular, any trajectory in 
average spends in any subset fcT 2 a time proportional to the area (measure) 
of the subset U. Surprisingly, for surfaces of higher genera a directional flow 
can be minimal but not ergodic! Sometimes it is possible to find some special 
direction with the following properties. The flow in this direction is minimal. 
However, the flat surface S might be decomposed into a disjoint union of 
several subsets Vi of positive measure in such a way that some trajectories 
of the directional flow prefer one subset to the others. In other words, the 
average time spent by a trajectory in the subset Vi is not proportional to 
the area of Vi anymore. (The original ideas of such examples appear in [Vcl], 
[Katl], [Sat] [Kea2]; see also [MaT] and especially [Ma7] for a very accessible 
presentation of such examples.) 

Suppose that we managed to find a direction on the initial surface Sq such 
that the flow in this direction is minimal but not ergodic (with respect to 
the natural Lebesgue measure). Let us apply a rotation to So which would 
make the corresponding direction vertical. Consider the resulting flat surface 
S (see Convention 1 in Sec. 1.2). Consider the corresponding "point" S G 
H(dx, . . . , d m ) and the orbit {gtS}teM of S under the action of the diagonal 



Recall that the stratum (or, more precisely, the corresponding "unit hyper- 
boloid") H\(di, . . . , d m ) is never compact, it always contains "cusps": regions 
where the corresponding flat surfaces have very short saddle connections or 
very short closed geodesies (see Sec. 3.2). 

Theorem (H. Masur). Consider a flat surface S. If the vertical flow is 
minimal but not ergodic with respect to the natural Lebesgue measure on the 
flat surface then the trajectory g t S of the Teichmuller geodesic flow is diver- 
gent, i.e. it eventually leaves any fixed compact subset K C 7i\{di, . . . , d m ) in 
the stratum. 



subgroup gt 
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Actually, this theorem has an even stronger form. 

A stratum . . . , d m ) has "cusps" of two different origins. A flat sur- 

face may have two distinct zeroes get very close to each other. In this case S 
has a short saddle connection (or, what is the same, a short relative period). 
However, the corresponding Ricmann surface is far from being degenerate. 
The cusps of this type correspond to "simple noncompactness" : any stratum 
Hi(di, . . . , d m ) is adjacent to all "smaller" strata TLi(d\ + e?2, e?3, ■ ■ ■ , dm), 

Another type of degeneration of a fiat surface is the appearance of a short 
closed geodesic. In this case the underlying Riemann surface is close to a 
degenerate one; the cusps of this second type correspond to "essential non- 
compactness" . 

To formulate a stronger version of the above Theorem consider the natural 
projection of the stratum 7i{d\, . . . , d m ) to the moduli space M. g of complex 
structures (see (4) in Sec. 3.3). Consider the image of the orbit {gtS}tem. in 
M g under this natural projection. By the reasons which we explain in Sec. 8 
it is natural to call this image a Teichmiiller geodesic. 

Theorem (H. Masur). Consider a flat surface S. If the vertical flow is 
minimal but not ergodic with respect to the natural Lebesgue measure on the 
flat surface then the "Teichmiiller geodesic" gtS is divergent, i.e. it eventu- 
ally leaves any fixed compact subset K C M. in the moduli space of complex 
structures and never visits it again. 

This statement (in a slightly different formulation) is usually called Ma- 
sur 's criterion of unique ergodicity (see Sec. A for discussion of the notion 
unique ergodicity). 

As a second illustrations of the "general philosophy" we present a combi- 
nation of Veech criterion and of a Theorem of J. Smillic. 

Recall that closed regular geodesies on a flat surface appear in families of 
parallel closed geodesies. When the flat surface is a flat torus, any such fam- 
ily covers all the torus. However, for surfaces of higher genera such families 
usually cover a cylinder filled with parallel closed geodesic of equal length. 
Each boundary of such a cylinder contains a conical point. Usually a geodesic 
emitted in the same direction from a point outside of the cylinder is dense in 
the complement to the cylinder or at least in some nontrivial part of the com- 
plement. However, in some rare cases, it may happen that the entire surface 
decomposes into several cylinders filled with parallel closed geodesies going 
in some fixed direction. This is the case for the vertical or for the horizontal 
direction on the flat surface glued from a regular octagon, see Fig. 3 (please 
check). Such direction is called completely periodic. 

Theorem (J. Smillie; W. A. Veech). Consider a flat surface S. If its 
GL + (2, M.)-orbit is closed in TL{d\, . . . , d m ) then a directional flow in any di- 
rection on S is either completely periodic or uniquely ergodic. 

(see Sec. A for the notion of unique ergodicity). Note that unique ergodicity 
implies, in particular, that any orbit which is not a saddle connection, (i.e. 
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which does not hit the singularity both in forward and in backward direction) 
is everywhere dense. We shall return to this Theorem in Sec. 9.5 where we 
discuss Veech surfaces. 

4 How Do Generic Geodesies Wind Around Flat 
Surfaces 

In this section we study geodesies on a flat surface S going in generic directions 
on S. Such geodesies are dense on S; moreover, it is possible to show that they 
wind around S in a relatively regular manner. Namely, it is possible to find a 
cycle c e H\(S;M.) such that in some sense a long piece of geodesic pretends 
to wind around 5* repeatedly following this asymptotic cycle c. 

In Sec. 4.1 we study the model case of the torus and give a rigorous defini- 
tion of the asymptotic cycle. Then we study the asymptotic cycles on general 
flat surfaces. In Sec. 4.2 we study "further terms of approximation". The 
asymptotic cycle describes the way in which a geodesic winds around the 
surface in average. In Sec. 4.2 we present an empirical description of the de- 
viation from average. The corresponding rigorous statements are formulated 
in Sec. 4.3. Some ideas of the proof of this statement are presented in Sec. 5. 

4.1 Asymptotic Cycle 

Asymptotic Cycle on a Torus 

As usual we start from the model case of the torus. We assume that our 
flat torus is glued from a square in the natural way. Consider an irrational 
direction on the torus; any geodesic going in this direction is dense in the 
torus. 

Fix a point xq on the torus and emit a geodesic in the chosen direction. 
Wait till it winds for some time around the torus and gets close to the initial 
point xq. Join the endpoints of the resulting piece of geodesic by a short path. 
We get a closed loop on the torus which defines a cycle c\ in the first homology 
group f/i(T 2 ;Z) of the torus. Now let the initial geodesic wind around the 
torus for some longer time; wait till it get close enough to the initial point xq 
and join the endpoints of the longer piece of geodesic by a short path. We get 
a new cycle C2 € i?i(T 2 ; Z). Considering longer and longer geodesic segments 
we get a sequence of cycles Ci £ Hi(T 2 ; Z). 

For example, we can choose a short segment X going through xq orthogo- 
nal (or just transversal) to the direction of the geodesic. Each time when the 
geodesic crosses X we join the crossing point with the point xq along X ob- 
taining a closed loop. Consecutive return points x\,X2, ■ ■ ■ define a sequence 
of cycles c\, c%, . . . , see Fig. 16. 

For the torus case we can naturally identify the universal covering space 
M 2 -> T 2 with the first homology group #i(T 2 ;]R) ~ R 2 . Our irrational 
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Fig. 16. A sequence of cycles approximating a dense geodesic on a torus 



geodesic unfolds to an irrational straight line Vi in IR 2 and the sequence of 
cycles ci, C2, . . . becomes a sequence of integer vectors Vi, V2, • • • € C R 2 
approximating Vi, see Fig. 16. 

In particular, it is not surprising that there exists the limit 

lim 71 — T = c ( 6 ) 

Tl^OO || C „|| 

Under our identification Hi(T 2 ;M.) ~ M 2 the cycle c represents a unit vector 
in direction Vi . 

Let the area of the torus be normalized to one. Let the interval X, which 
we use to construct the sequence c±, C2, ■ • ■ , be orthogonal to the direction of 
the geodesic. Denote by \X\ its length. The following limit also exists and is 
proportional to the previous one: 

lim i Cn = -|- ■ c (7) 

The cycles obtained as limits (6) and (7) are called asymptotic cycles. They 
show how the corresponding irrational geodesic winds around the torus in 
average. It is easy to see that they do not depend on the starting point x. 
The notion "asymptotic cycle" was introduced by S. Schwartzman [Schw]. 

Asymptotic Cycle on a Surface of Higher Genus 

We can apply the same construction to a geodesic on a flat surface S of 
higher genus. Having a geodesic segment X C S and some point x G X we 
emit from x a geodesic orthogonal to X. From time to time the geodesic 
would intersect X. Denote the corresponding points as x\, x%, ■ ■ ■ ■ Closing up 
the corresponding pieces of the geodesic by joining the endpoints Xo,Xj with 
a path going along X we again get a sequence of cycles c\ , c<i , . . . . 

Proposition 1. For any flat surface S of area one and for almost any di- 
rection a on it any geodesic going in direction a is dense on S and has an 
asymptotic cycle which depends only on a. 

In other words, for almost any direction the limit 
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exists and the corresponding asymptotic cycle c does not depend on the starting 
point x G S. 

This proposition is an elementary corollary from the following theorem of 
S. Kerckhoff, H. Masur and J. Smillie [KMaS] (which is another Key Theorem 
in this area). 

Theorem (S. Kerckhoff, H. Masur, J. Smillie). For any flat surface S 
the directional flow in almost any direction is ergodic. 

In this case the asymptotic cycle has the same dynamical interpretation 
as for the torus: it shows how a geodesic going in the chosen direction winds 
around the surface S in average. 

Remark. Note that the asymptotic cycle c S Hi(S, R) also has a topological 
interpretation. Assume for simplicity that c corresponds to the vertical direc- 
tion. Let u> be the holomorphic 1-form corresponding to the flat structure on 
S (see Sec. 3.3). Then the closed 1-form ujo = Re(w) defines the vertical folia- 
tion and c = -D[a>o] is Poincare dual to the cohomology class of loq. Choosing 
other ergodic directions on the flat surface S we get asymptotic cycles in the 
two-dimensional subspace (D[u!o\, D[lui])t^ c H\(S, R) spanned by homology 
classes dual to cocycles loq = Re(o->) and u\ = Im(cj). 

4.2 Deviation from Asymptotic Cycle 

We have seen in the previous section that a sequence of cycles c\ , c% , . . . ap- 
proximating long pieces of an "irrational" geodesic on a flat torus T 2 and on 
a flat surface S of higher genus exhibit similar behavior: their norm grows 
(approximately) linearly in n and their direction approaches the direction of 
the asymptotic cycle c. Note, however, that for the torus the cycles c„ live 
in the two-dimensional space iii(T 2 ;R) ~ R 2 , while for the surface of higher 
genus g > 2 the cycles live in the larger space i?i(5;R) ~ R 2fl . In particular, 
they have "more room" for deviation from the asymptotic direction. 

Namely, observing the right part of Fig. 16 we see that all vectors c n 
follow the line Vi spanned by the asymptotic cycle c rather close: the norm 
of projection of c n to the line orthogonal to Vi is uniformly bounded (with 
respect to n and to the choice of the starting point xo). 

The situation is different for surfaces of higher genera. Choose a hyperplane 
S _L c in Hi(S, R) as a screen orthogonal (transversal) to the asymptotic cycle 
c and consider a projection to this screen parallel to c. Projections of cycles 
c„, would not be uniformly bounded anymore. There is no contradiction since 
if the norms of these projections grow sublinearly, then the directions of the 
cycles c n still tend to direction of the asymptotic cycle c. 

Let us observe how the projections are distributed in the screen S. Fig- 
ure 17 shows results of numerical experiments where we take a projection of 
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Fig. 17. Projection of a broken line joining the endpoints of ci, C2, ciooooo to a 
screen orthogonal to the asymptotic cycle. Genus g = 3 

a broken line joining the endpoints of c\, c%, . . . , ciooooo and we take a two- 
dimensional screen orthogonal to c to make the picture more explicit. 

We see that the distribution of projections of cycles c ra in the screen S is 
anisotropic: the projections accumulate along some line. This means that in 
the original space M. 2g the vectors c„ deviate from the asymptotic direction Vi 
not arbitrarily but along some two-dimensional subspace Vi D Vi, see Fig. 18. 




.'•2 



Fig. 18. Deviation from the asymptotic direction 

Moreover, measuring the norms of the projections proj(c n ) to the screen 
S orthogonal to L\ = (c)r, we get 

log||proj(c„)|| 
km sup = Vi < 1 

n~*oo lOgn 
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In other words the vector c n is located approximately in the subspace V2, and 
the distance from its endpoint to the line V\ C V2 is bounded by const - 1 Ctt. 1 1 ^ 2 , 
see Fig. 18. 

Consider now a new screen £2 -L V2 orthogonal to the plane V2- Now the 
screen S 2 has codimension two in H\{S, R) ~ R 2s . Considering the projections 
of c„ to £2 we eliminate the asymptotic directions Vi and V2 and we see how 
do the vectors c n deviate from Vi. On the screen 1S2 we see the same picture 
as in Fig. 17: the projections are located along a one-dimensional subspace. 

Coming back to the ambient space Hi(S, K) ~ R 2ff , this means that in 
the first term of approximation all vectors c n are aligned along the one- 
dimensional subspace Vi spanned by the asymptotic cycle. In the second term 
of approximation, they can deviate from Vi, but the deviation occurs mostly 
in the two-dimensional subspace V2, and has order 1 1 c | j ^ 2 where V2 < 1. In the 
third term of approximation we see that the vector c n may deviate from the 
plane V2 , but the deviation occurs mostly in a three-dimensional space V3 and 
has order \\c\\ U3 where 1^3 < v^. 

Going on we get further terms of approximation. However, getting to a 
subspace V g which has half the dimension of the ambient space we shall see 
that, in a sense, there is no more deviation from V g : the distance from any c n 
to V g is uniformly bounded. 

Note that the intersection form endows the space Hi(S,M) ~ M. 2g with 
a natural symplectic structure. It can be checked that the resulting g- 
dimcnsional subspace V g is a Lagrangian subspace for this symplectic form. 

4.3 Asymptotic Flag and "Dynamical Hodge Decomposition" 

A rigorous formulation of phenomena described in the previous section is given 
by the following Theorem proved 3 by the author in [Zo3] and [Zo4] . 

Following Convention 1 we always consider a flat surface together with 
a choice of direction which by convention is called the vertical direction, or 
direction to the North. Using an appropriate homothety we normalize the area 
of S to one, so that S G Hi(dx,... ,d m ). 

We chose a point xq € S and a horizontal segment X passing through x$; 
by \X\ we denote the length of X. We consider a geodesic ray 7 emitted from 
Xo in the vertical direction. (If xq is a saddle point, there are several outgoing 
vertical geodesic rays; choose any of them.) Each time when 7 intersects X we 
join the point x n of intersection and the starting point xo along X producing 
a closed path. We denote the homology class of the corresponding loop by c n . 

Let w be the holomorphic 1-form representing S; let g be genus of S. 
Choose some Euclidean metric in -ff^S 1 ; R) ~ R 2fl which would allow to mea- 
sure a distance from a vector to a subspace. Let by convention log(O) = —00. 

3 Actually, the theorem was initially proved under certain hypothesis on the Lya- 
punov exponents of the Teichmiiller geodesic flow. These conjectures were proved 
later by G. Forni and in the most complete form by A. Avila and M. Viana; see the 
end of this section and especially Sec. 5.8 
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Theorem. For almost any fiat surface S in any stratum TLi (d\ , . . . , d m ) there 
exists a flag of subspaces 

Vi C V 2 c ••• C V g C £Ti(S;R) 

in the first homology group of the surface with the following properties. 

Choose any starting point xq € X in the horizontal segment X . Consider 
the corresponding sequence ci, c%, . . . of cycles. 

— The following limit exists 

\X\ lim - c„ = c, 

where the nonzero asymptotic cycle c G Hi(Mg,M) is Poincare dual to the 
cohomology class of loq — Re[w], and the one- dimensional subspace Vi = (c)r 
is spanned by c. 

— For any j = 1, . . . , g — 1 one /ias 

log o!ist(c„,V i ) 

ImiSUp : — = Vj+l 

n-*oo logn 

and 

\dist(c n , V g )\ < const, 

where the constant depends only on S and on the choice of the Euclidean 
structure in the homology space. 

The numbers 2, 1 + 1/2, . . . , 1 + V g are the top g Lyapunov exponents of the 
Teichmuller geodesic flow on the corresponding connected component of the 
stratum H(di, . . . ,d m ); in particular, they do not depend on the individual 
generic flat surface S in the connected component. 

A reader who is not familiar with Lyapunov exponents can either read 
about them in Appendix B or just consider the numbers Vj as some abstract 
constants which depend only on the connected component TC comp (di, . . . ,d m ) 
containing the flat surface S. 

It should be stressed, that the theorem above was initially formulated 
in [Zo4] as a conditional statement: under the conjecture that v g > there 
exist a Lagrangian subspace V g such that the cycles are in a bounded distance 
from V g ; under the further conjecture that all the exponents fj, for j = 
1,2, . . . , g, arc distinct, there is a complete Lagrangian flag (i.e. the dimensions 
of the subspaces Vj, where j = l,2,...,g, rise each time by one). These 
two conjectures were later proved by G. Forni [Fori] and by A. Avila and 
M. Viana [AvVi]. We discuss their theorems in Sec. 5.8. 

Another remark concerns the choice of the horizontal segment X . By con- 
vention it is chosen in such way that the trajectories emitted in the vertical 
direction (in direction to the North) from the endpoints of X hit the conical 
points before the first return to X. Usually we just place the left endpoint of 
X at the conical point. 
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Omitting this condition and considering a continuous family of horizon- 
tal subintervals X t of variable length (say, moving continuously one of the 
endpoints), the theorem stays valid for a subset of X t of full measure. 

5 Renormalization for Interval Exchange 
Transformations. Rauzy— Veech Induction 

In this section we elaborate a powerful time acceleration machine which allows 
to study the asymptotic cycles described in Sec. 4. Following the spirit of this 
survey we put emphasis on geometric ideas and omit proofs. This section can 
be considered as a geometric counterpart of the article of J.-C. Yoccoz [Y] in 
the current volume. 

I use this opportunity to thank M. Kontsevich for numerous ideas and 
conjectures which were absolutely crucial for my impact in this theory: without 
numerous discussions with M. Kontsevich papers [Zo2] and [Zo3], probably, 
would be never written. 

5.1 First Return Maps and Interval Exchange Transformations 

Our goal is to study cycles obtained from long pieces of "irrational" geodesic 
on a flat surface by joining their endpoints along a transversal segment X. 
To perform this study we elaborate some simple machine which generates the 
cycles, and then we accelerate this machine to obtain very long cycles in a 
rather short time. 

Consider all geodesies emitted from a transverse segment X in the same 
generic direction and let each of them come back to X for the first time. We 
get a first return map T : X — ► X which is interesting by itself and which 
deserves a separate discussion. Its properties play a crucial role in our study. 
(See Appendix A for general properties of the first return map.) 

As usual let us start with a model case of a flat torus. Take a meridian 
of the torus as a transversal X and emit from X a directional flow. Every 
geodesic comes back to X inducing the first return map T : X — > X, which in 
this case isometrically rotates the meridian X along itself by an angle which 
depends on the direction of the flow, see Fig. 19. 

Assume that our flat torus T 2 is glued from a unit square. Let us replace 
now a meridian of the torus by a generic geodesic segment X orthogonal to the 
direction of the flow. From every point of X we emit a geodesic in direction 
orthogonal to X and wait till it hits X for the first time. We again obtain 
a first return map T : X — ► X, but this time the map T is slightly more 
complicated. 

To study this map it is convenient to unfold the torus into a plane. The 
map T is presented at Fig. 20. It chops X into three pieces and then shuffles 
them sending the left subinterval to the right, the right subinterval to the left 
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Fig. 19. The first return map of a meridian to itself induced by a directional flow 
is just a twist 

and keeping the middle one in the middle but shifting it a bit. The map T 
gives an example of an interval exchange transformation. 




Fig. 20. Directional flow on a torus. The first return map of a segment to itself is 
an interval exchange transformation of three subintervals 

Note that when the direction of the flow is irrational, the geodesies emit- 
ted from X again cover the entire torus T 2 before coming back to X. The 
torus is get ripped into three rectangles based on the three subintervals in 
which T chops X. The corresponding building of three rectangles gives a new 
fundamental domain representing the torus: one can see at Fig. 20 that it tiles 
the plane. Initially we glued our flat torus from a square; the building under 
consideration gives another way to unwrap T 2 into a polygon. We recommend 
to the reader to check that identifying the two pairs of corresponding vertical 
sides of the building and gluing the top horizontal sides of the rectangles to 
the bottom of X as prescribed by the interval exchange transformation T we 
get the initial torus. 
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Consider now a flat surface of genus higher than one. Say, consider a fiat 
surface of genus g = 2 as on Fig. 21. We suggest to the reader to check 
that this flat surface has a single conical singularity with a cone angle 6tt 
(see Fig. 2). To study a directional flow choose as before a geodesic segment 
X C S orthogonal to the direction of the flow and consider the first return map 
T : X — > X induced by the flow; see Fig. 21. We see that X is chopped into 
a larger number of subintervals (in comparison with the torus case), namely 
for our choice of X it is chopped into four subintervals. 

Now we observe a new phenomenon: trajectories emitted from some points 
of X hit the conical point and our directional flow splits at this point. Since 
in our particular case the cone angle at the conical point is 6tt = 3 • 2tt there 
are three trajectories in direction v which hit it. The corresponding points at 
which X is chopped are marked with bold dots. The remaining discontinuity 
point of X corresponds to a trajectory which hits the endpoint of X. 

Our construction with a segment X transversal to the flow and with tra- 
jectories of the flow emitted from X and followed till their first return to X 
trims a braid from the flow. Conical points play the role of a comb which 
splits the flow into several locks and then trims them in a different order. 
Note, however, that if we follow the flow till the second return to X it will 
pass through the comb twice, and thus will be generically split already into 
seven locks. (If you are interested in details, think why this second return has 
seven and not eight locks and what sort of genericity we need). 

Similarly, the interval exchange transformation T of the base interval X 
can be compared to a shuffling machine. Imagine that X represents a stock of 
cards. We split the stock into n parts of fixed widths and shuffle the parts in 
a different order (given by some permutation tt of n elements). At the second 
iteration we again split the new stock in the parts of the same widths as 
before and shuffle the parts according to the same permutation tt, etc. Note, 
that even if the permutation tt is such that tt 2 = id, say, tt = (4, 3, 2, 1), the 
second iteration T 2 is not an identical transformation provided the widths 
Ai, . . . , A4 are not symmetric: for a generic choice of Ai, . . . , A4 the interval 
exchange transformation T 2 has 6 discontinuities (and hence 7 subintervals 
under exchange). 

Exercise. Consider an interval exchange transformation T(A, tt) corresponding 
to the permutation tt = (4, 3, 2, 1). Choose some generic values of the lengths 
Ai, . . . , A4 of subintervals and construct T 2 and T 3 . 

5.2 Evaluation of the Asymptotic Cycle Using an Interval 
Exchange Transformation 

Now we can return to our original problem. We want to study long pieces 
of leaves of the vertical foliation. Fix a horizontal segment X and emit a 
vertical trajectory from some point x € X . When the trajectory intersects X 
for the first time join the corresponding point T(x) to the original point x 
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Fig. 21. The first return map T : X — > X of a geodesic segment X defined 
by a directional flow decomposes the surface into four rectangles "zippered" along 
singular trajectories 

along X to obtain a closed loop. Here T : X — > X denotes the first return 
map to the transversal X induced by the vertical flow. Denote by c(x, 1) the 
corresponding cycle in Hi(S;Wi). Following the vertical trajectory further on 
we shall return to X once again. Joining x and the point T(T(x)) of the second 
return to X along X we obtain the second cycle c(x,2). We want to describe 
the cycle c(x, N) obtained after a very large number N of returns. 

Actually, we prefer to close up a piece of trajectory going from x G X to the 
first return point T(x) £lina slightly different way. Instead of completing 
the path joining the endpoints it is more convenient to close this piece of 
trajectory joining both points x and T(x) to the left endpoint of X along X 
(see Fig. 22). This modified path defines the same homology cycle c(x, I) as 
the closed path for which the points x and T(x) arc joined directly. 

Consider now the "first return cycle" c(x, 1) as a function c(x) = c{x, I) of 
the starting point x € X . Let the interval exchange transformation T : X — > X 
decompose X into n subintervals Xi U • • • U X n . It is easy to see that the 
function c(x) is piecewise constant: looking at Fig. 22 one can immediately 
verify that if two points x\ and xi are not separated by a discontinuity point 
(i.e. if they belong to the same subinterval Xj) they determine homologous 
(and even homotopic) cycles c(x\) = c{x2). Each subinterval Xj determines 
its own cycle c(Xj). 
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Fig. 22. Decomposition of a long cycle into the sum of basic cycles. (We have 
unfolded the fiat surface along vertical trajectory emitted from the point x) 



The vertical trajectory emitted from a point x <G X returns to X at 
the points T(x), T 2 (x), . . . , T N (x). It is easy to see that the cycle c(x, 2) 
corresponding to the second return can be represented as a sum c(x, 2) = 
c(x) + c(T(x)), see Fig. 22. Similarly the cycle c(x,N) obtained by closing 
up a long piece of vertical trajectory emitted from x <G X and followed up to 
N-tYi return to X can be represented as a sum 

c(x, N) = c(x) + c(T{x)) + ■■■ + c{T N -\x)) (8) 

According to the fundamental Theorem of S. Kcrckhoff, H. Masur and 
J. Smillie [KMaS], for any flat surface the directional flow in almost any di- 
rection is ergodic, and even uniquely ergodic. Hence, the same is true for the 
corresponding interval exchange transformation. Applying the ergodic theo- 
rem (see Appendix A) to the sum (8) and taking into consideration that c(x) 
is a piecewise-constant function we get 

1 f 1 

c{x, N) ~ N ■ — J c(x) dx = N-— (A 1 c(X 1 ) + • • • + A„c(X n )) 

where c(Xj) denotes the "first return cycle" for the points x in the subinterval 
Xj, see Fig. 22. This gives an explicit formula for the asymptotic cycle 

c(x N) 1 / \ 
c = Jim = _ + • • • + Kc(X n )) (9) 

Note that the asymptotic cycle does not depend on the starting point x G X. 

Exercise. Show that the paths vx, V4 (see Fig. 21) represent a basis of 
cycles of the corresponding flat surface S of genus g = 2. Show that the first 
return cycles c{Xj) determine another basis of cycles and that one can pass 
from one basis to the other using the following relations (see Fig. 21): 
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c(Xi) = vi - v 4 c(X 2 ) = v x - v 3 - v 4 

c(X 3 ) = 2vi + v 2 - v 4 c(X 4 ) = vi + v 2 + v 3 - v 4 

Express the asymptotic cycle in the basis Vj in terms of the lengths and 
then in terms of the angle by which the regular octagon is turned with respect 
to the standard presentation. Locally the vertical foliation goes to the North. 
And globally? 

5.3 Time Acceleration Machine (Renormalization): Conceptual 
Description 

In the previous section we have seen why all trajectories of a typical directional 
flow wind around the surface following the same asymptotic cycle. We have 
also found an effective way to evaluate this asymptotic cycle: we have seen 
that it is sufficient to find an interval exchange transformation T : X — ► X 
induced on any transverse segment X as the first return map of the directional 
flow, and then to determine the "first return cycles" c(Xj), see Fig. 22. The 
linear combination of the cycles c(Xj) taken with weights proportional to the 
lengths Xj = \Xj\ of subintervals gives the asymptotic cycle, see (9). 

Let us proceed now with a more delicate question of deviation of a tra- 
jectory of the directional flow from the asymptotic direction. Without loss of 
generality we may assume that the directional flow under consideration is the 
vertical flow. We know that a very long cycle c(x, N) corresponding to a large 
number N of returns of the trajectory to the horizontal segment X stretches 
in the direction approaching the direction of the asymptotic cycle c. We want 
to describe how c(x, N) deviates from this direction (see Sec. 4.2 and Sec. 4.3). 

We have already seen that as soon as we have evaluated the "first return 
cycles" c(Xj), complete information about cycles representing long pieces of 
trajectories of the directional flow is encoded in the corresponding trajectory 
x,T(x), . . . ,T N ~ 1 (x) of the interval exchange transformation; see (8). 

An interval exchange transformations gives an example of a parabolic dy- 
namical system which is neither completely regular (like rotation of a circle) 
nor completely chaotic (like geodesic flow on a compact manifold of constant 
negative curvature, which is a typical example of a hyperbolic system). In some 
aspects interval exchange transformations are closer to rotations of a circle: 
say, as it was proved by A. Katok, an interval exchange transformation is 
never mixing [Kat2] (though, as it was very recently proved by A. Avila and 
G. Forni in [AvFor], gcncrically it is weakly mixing). However, the behavior of 
deviation from the ergodic mean resembles the behavior of a chaotic system. 

Our principal tool in the study of interval exchange transformations ex- 
ploits certain self-similarity of these maps. Choosing a shorter horizontal in- 
terval X' we make the vertical flow wind for a long time before the first 
return to X' . However, the new first return map in a sense would not be 
more complicated than the initial one: it would be again an interval exchange 
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transformation T 1 : X' — > X' of the same (or almost the same) number of 
subintervals. 

To check the latter statement let us study the nature of the points of 
discontinuity of the first return map T : X — > X. An interior point iglis 
a point of discontinuity either if the forward vertical trajectory of x hits one 
of the endpoints of X (as on Fig. 20) or if the forward vertical trajectory of x 
hits the conical point before coming back to X (see Fig. 21). A conical point 
having the cone angle 2ir(d + 1) has d+ 1 incoming vertical trajectories which 
land to this conical point. (Say, the flat surface represented on Fig. 21 has 
a single conical point with the cone angle 6ir; hence this conical point has 3 
incoming vertical trajectories.) Following them at the backward direction till 
the first intersection with X we find d + 1 points of discontinuity on X (see 
Fig. 21). Thus, all conical points taken together produce + 1) points of 

discontinuity on X. 

Generically two more points of discontinuity come from the backward tra- 
jectories of the endpoints of X. However, in order to get as small number of 
discontinuity points as possible we can choose X in such way that either back- 
ward or forward trajectory of each of the two endpoints hits some conical point 
before coming back to X. This eliminates these two additional discontinuity 
points. 

Convention 2. From now on we shall always choose any horizontal subinterval 
X in such way that the interval exchange transformation T : X — > X induced 
by the first return of the vertical flow to X has the minimal possible number 

n — ^^(rfj + 1) + 1 = 2g + (number of conical points) — 1 

3 

of subintervals under exchange. 

In the formula above we used the Gauss-Bonnet formula telling that 
J2j dj — 2g — 2, where g is the genus of the surface. 

Following Convention 2 we shall usually place the left endpoint of the 
horizontal interval X at the conical singularity. This leaves a discrete choice 
for the position of the right endpoint. 

Renormalization 

We apply the following strategy in our study of cycles c(x, N). Choose some 
horizontal segment X satisfying Convention 2. Consider vertical trajectories, 
which hit conical points. Follow them in backward direction till the first in- 
tersection with X. Consider the resulting decomposition X = X\ U • • • U X n , 
the corresponding interval exchange transformation T : X — > X and the "first 
return cycles" c(Xj). 

Consider a smaller subinterval X' C X satisfying Convention 2. Apply 
the above procedure to X'; let X ' = X[ U ■ • • U X' n be the corresponding 
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decomposition of X' . We get a new partition of our flat surface into a collection 
of n rectangles based over subintervals X[ U • • • U X' n . 

By construction the vertical trajectories of any two points x$,x G X' k 
follow the same high and narrow rectangle R' k of the new building up to 
their first return to X' . This implies that the corresponding new "first return 
cycles" c'(xo) = c'(x) are the same and equal to c'{X' k ). 

Both vertical trajectories of Xq,x £ X' k intersect the initial interval X 
many times before first return to X'. However, since these trajectories stay 
together, they visit the same intervals Xj k in the same order jo, Ji, ■ ■ ■ ,ji (the 
length I = l(k) of this trajectory depends on the subinterval X' k ). 

This means that we can construct an n x n-matrix Bjk indicating how 
many times a vertical trajectory emitted from a point x € X' k have visited 
subinterval Xj before the first return to X' . (By convention the starting point 
counts, while the first return point does not.) Here X = Xi U ■ ■ • U X n is the 
partition of the initial "long" horizontal interval X and X' = X[ U ■ ■ • U X' n 
is the partition of the new "short" subinterval X'. 

Having computed this integer matrix B we can represent new "first return 
cycles" c'(X' k ) in terms of the initial "first return cycles" c(Xj) as 

c'(X' k ) = Bi fc c(Xi) + • • • + B nk c(X n ) (10) 

Moreover, it is easy to see that the lengths X' k = \X' k \ of subintervals of the new 
partition are related to the lengths \Xj\ of subintervals of the initial partition 
by a similar relation 

\. /-..a; (ii) 

Note that to evaluate matrix B we, actually, do not need to use the vertical 
flow: the matrix B is completely determined by the initial interval exchange 
transformation T : X — > X and by the position of the subinterval X' C X. 

What we gain with this construction is the following. To consider a cycle 
c(x, N) representing a long piece of leaf of the vertical foliation we followed 
the trajectory x, T(x), . . . , T N (x) of the initial interval exchange transforma- 
tion T : X — > X and applied formula (8). Passing to a shorter horizontal 
interval X' C X we can follow the trajectory x,T'(x), . . . , (T') N (x) of the 
new interval exchange transformation T' : X' — > X' (provided x <E X 1 ). Since 
the subinterval X' is much shorter than X we cover the initial piece of tra- 
jectory of the vertical flow in a smaller number N' of steps. In other words, 
passing from T to T' we accelerate the time: it is easy to see that the tra- 
jectory x,T'(x),...,(T') N '(x) follows the trajectory x, T(x), . . . , T N (x) but 
jumps over many iterations of T at a time. 

Of course this approach would be non efficient if the new first return map 
T 1 : X' — * X' would be much more complicated than the initial one. But we 
know that passing from T to T' we stay within a family of interval exchange 
transformations of the fixed number n of subintervals, and, moreover, that the 
new "first return cycles" and the lengths of the new subintervals are expressed 
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in terms of the initial ones by means of the n x n-matrix B, which depends 
only on the choice of X' C X and which can be easily computed. 

Our strategy can be formalized as follows. In the next two sections we 
describe a simple explicit algorithm (generalizing Euclidean algorithm) called 
Rauzy-Veech induction which canonically associates to an interval exchange 
transformation T : X — > X some specific subinterval X' C X and, hence, a 
new interval exchange transformation T' : X' — > X' . This algorithm can be 
considered as a map from the space of all interval exchange transformations of 
a given number n of subintervals to itself. Applying recursively this algorithm 
we construct a sequence of subintervals X = X^ D X^ D X^ D . . . and 
a sequence of matrices B = B(X^), B(X^), . . . describing transitions form 
interval exchange transformation : X^ — > X^ to interval exchange 
transformation T^ 1 ) : X^+V -> X^ r+1 \ Rewriting equations (10) and (11) 
in a matrix form we get: 

H.uy.p| 

H.U-.H 

\\ n (x^))J V / W* (r) )/ 

Taking a product B^ = B(X^) ■ B(X^) B(X^ S '^) we can imme- 
diately express the "first return cycles" to a microscopic subinterval X^ in 
terms of the initial "first return cycles" to X by a linear expression analogous 
to (12). Note, however, that before coming back to this microscopic subin- 
terval X^ the vertical flow has to travel for enormously long time. The first 
return cycle to this very short subinterval X^ represents the cycle c(x, N) 
corresponding to very long trajectory x,T(x), ...,T N (x) of the initial interval 
exchange transformation with N ~ cxp(const ■ s). In other words, our renor- 
malization procedure plays a role of a time acceleration machine: instead of 
following patiently the trajectory x,T(x), ...,T N (x) of the initial interval ex- 
change transformation for the exponential time N ~ cxp(const ■ s) we obtain 
the cycle c(x, N) applying only s steps of renormalization! 

One can argue that in this way we can describe only very special parts of 
vertical trajectories: those which start and end at the same microscopically 
small subinterval X^ C X . This can be overdone by the following technique. 
Consider an enormously long trajectory x, T(x), . . . , T N (x) which starts and 
finishes at some generic points of X . One can choose s(N) in such way that the 
trajectory would get to X^ relatively soon (in comparison with its length TV); 
then would return back to X^ many times; and would reach the last point 
T N (x) relatively fast after the last visit to X 1 --^. That means that essentially 
(up to negligibly short starting part and ending part) one can assume that the 
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entire trajectory starts at and ends at X^ (returning to this subinterval 
many times). 

This simple idea can be developed and rigorously arranged (sec [Zo4] for 
details). To avoid overloading of this survey with technicalities we consider 
only a simplified problem giving a comprehensive description of the first return 
cycles to X^ s \ The nature of the asymptotic flag is especially transparent in 
this case. 

5.4 Euclidean Algorithm as a Renormalization Procedure in 
Genus One 

To illustrate the idea of renormalization we start with the "elementary" case, 
when the Riemann surface is a torus, the foliation is a standard irrational 
foliation, and the initial transversal X is a meridian. We have seen at Fig. 19 
that in this case the first return map T : X — > X is just a rotation of a circle. 

Consider rotation of a circle T : S 1 — > S 1 by an angle a. Let the length of 
the circle be normalized to one. Consider trajectory x, Tx, T 2 x, ... of a point 
x (see Fig. 23). Denote the length of the arc (x,Tx) by A = a/(2w). 




Fig. 23. Renormalization for rotation of a circle leads to Euclidean algorithm and 
to Gauss measure 

Cutting the circle at the point x we get an interval X; the rotation of 
the circle generates a map of the interval X to itself which we denote by 
the same symbol T : X — ► X. Unbend X isometrically to get a horizontal 
interval of unit length in such way that the counterclockwise orientation of 
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the circle gives the standard positive orientation of the horizontal interval. 
The map T acts on X as follows: it cuts the unit interval X into two pieces 
Xi U X 2 of lengths \Xi\ = 1 — A and \X 2 \ = A and interchanges the pieces 
preserving the orientation, see Fig. 23. In other words, the map T is an interval 
exchange transformation of two subintcrvals. (To avoid confusion we stress 
that X 2 = [T- l x,x] and not X = [x,Tx].) 

Suppose now that we are looking at X in the microscope which shows 
only the subinterval X' = [x,Tx[ (corresponding to the sector of angle a at 
Fig. 23). Consider the trajectory x, Tx, T 2 x, ... of the point x which is the left 
extremity of X\ . For the particular rotation represented at Fig. 23 the points 
Tx, T 2 x, T 3 x, T A x are outside of the sector of our vision; the next point of the 
trajectory which we see in X' is the point T 5 x. This is the first return point 
T'x = T 5 x to the subinterval X 1 . Following the trajectory T 5 x, T 6 x, ... further 
on we would not see several more points and then we shall see T 10 x = T'(T'x). 
This is the second return to the subinterval X' . 

Note that the distance between x and T'x = T 5 x is the same as the 
distance between T'x = T 5 x and T'(T'x) = T w x; it equals (1 - {1/A}) • A, 
where { } denotes the fractional part of a real number. It is easy to see that 
T' : X' — > X' is again an interval exchange transformation of two subintcrvals 
X[ U X' 2 . The lengths of subintcrvals are \X[\ = {1/A} • A and (1 - {1/A}) • A. 
After identification of the endpoints the segment X' becomes a circle and the 
map T" becomes a rotation of the circle T". Having started with a rotation T 
in a counterclockwise direction by the angle a = 2n ■ A we get a rotation T' 

in a clockwise direction by the angle a' = 2tt • ^ — j (please verify). 

One should not think that X" = T 5 identically. It is true for the points of 
first subinterval, T'\x> = T 5 . However, for the points of the second subinterval 
X 2 we have T'\x> = T 4 . In other words, for the points of the sector a, which 
are close to the extremity Tx, four iterations of T bring them back to the 
sector. Thus, the matrix B(X') of number of visits to subintervals has the 
form 

[B u B 12 \ [A 3\ 
\B 2l B 22 ) 

(please draw X\ U X 2 and X[ U X 2 and verify). We remind that Bjk indicates 
how many times a vertical trajectory emitted from a point x G X' k have visited 
subinterval Xj before the first return to X', where by convention the starting 
point counts, while the first return point does not. 

Thus we get a renormalization procedure as described in the previous 
section: confine the map T to a smaller subinterval X'; consider the resulting 
first return map T"; rescale X' to have unit length. Having started with an 
interval exchange transformation T of two intervals of lengths (1 — A, A), where 
A € (0, 1) we get (after rescaling) an interval exchange transformation T' of 
two intervals of lengths {1/A}, 1 — {1/A}. Or, in terms of rotations, having 
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started with a counterclockwise rotation by the angle a = 2ir\ we get a 

clockwise rotation by the angle a' = 2n ■ / — j> . 

One can recognize Euclidean algorithm in our renormalization procedure. 
Consider the "space of rotations" , where rotations are parametrized by the 
angle 2tt\, A G [0; 1[. The map 



9^^{t> (13) 




can be considered as a map from "the space of rotations" to itself, or what 
is the same, a map from "the space of interval exchange transformations of 
two subintervals" to itself. The map g is ergodic with respect to the invariant 
probability measure 

1 dX 

on the parameter space A e [0;1[ which is called the Gauss measure. This 
map is intimately related with the development of A into a continued fraction 

1 

A = 



1 

ni H 



1 

n 2 H — 

ri3 H 

We shall see another renormalization procedure related to map (13) in the 
next sections, in particular, in Sec. 5.9. 



5.5 Rauzy Veech Induction 

In the previous section we have seen an example of a renormalization pro- 
cedure for interval exchange transformations of two intervals. In this section 
we consider a similar renormalization procedure which now works for interval 
exchanges of any number of subintervals. As we have seen in the previous sec- 
tion, we do not need to keep information about the flat surface to describe the 
renormalization algorithm. Nevertheless, we prefer to keep track of zippercd 
rectangles decomposition of the surface corresponding to the sequence of the 
horizontal subintervals X = D X^- 1 ' D . . . in order to preserve geometric 
spirit of the algorithm. 

Consider a flat surface S; choose a horizontal interval X satisfying Conven- 
tion 2; consider the corresponding decomposition of the surface into zippercd 
rectangles as on Fig. 21. Let X\ U • • • UX n be the corresponding decomposition 
of the horizontal segment X in the base; let \j = \Xj\ denote the widths of 
subintervals. 
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Convention 3. We associate to a decomposition of a flat surface into rectangles 
a permutation tt in such way that the top horizontal segments of the rectangles 
arc glued to the bottom side of the interval X in the order 7r _1 (l), . . . , 7r _1 (n). 

Example. In the example presented at Fig. 21 the four rectangles R%, . . . , R4 
appear at the bottom side of X in the order i?3, R±, i?4, R2, so we associate 
to this way of gluing a permutation 

(3142) =(3,1,4,2)- 1 = (2,4,1,3) = tt 

Exercise. Show that intersection indices c{X{) o c(Xj) of the "first return 
cycles" (sec Sec. 5.2) are given by the following skew-symmetric matrix Q(tt) 
defined by the permutation tt: 

!1 if i < j and 7r _1 (i) > ir~ 1 j 
-1 if i > j and 7r _1 (z) < TT~ l j (15) 
otherwise 

Evaluate Q(tt) for the permutation in the Example above and compare the 
result with a direct calculation for the cycles c(Xj), j = 1, ... ,4, computed 
in the Exercise at the end of Sec. 5.2. 

Compare now the width A„ of the rightmost rectangle R n with the width 
A^-if,;,) of the rectangle which is glued to the rightmost position at the bottom 
of X. As a new subinterval X 1 C X consider the subinterval X' , which has the 
same left extremity as X, but which is shorter than X by min(A n , A^-i^)). 

The situation when A„ > A„— is represented at Fig. 24; the situation 
when A n < A w -i(„) is represented at Fig. 25. 

By construction the first return map T' : X — > X' has the same number n 
of subintervals in its decomposition. Observing Fig. 24 and 24 one can see that 
in the first case, when A„ > A w -i(„), the new decomposition X[ U ■ ■ • U X' n is 
obtained from the original decomposition X\ U • • • U X n by shortening the last 
interval by A^-i^) from the right. In the second case, when A rl < A^-i^), the 
new decomposition X[ U • • • U X' n is obtained from the original decomposition 
X\ U • • • U X n by eliminating the last subinterval X n and by partitioning the 
subinterval X^-ir n \ into two ones of lengths A^— i( n ) — A„ and A„ correspond- 
ingly. 

The order in which the rectangles of the new building are glued to the 
bottom of the interval X' changes. The new permutation tt' can be described 
as follows. Consider the initial permutation tt as a pair of ordcrings of a 
finite set: a "top" ordering 1, 2, ...,n (corresponding to the ordering of the 
rectangles along the top side of the base interval X) and a "bottom" ordering 
7r _1 (l), . . . , 7t _1 (to) (corresponding to the ordering of the rectangles along the 
bottom side of the base interval X). In the first case, when A n > A T -i(„), the 
new permutation tt' corresponds to the modification of the bottom ordering 
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Fig. 24. Type I modification: the rightmost rectangle 7?4 on top of X is wider than 



the rectangle i?i = _R„ 



-1(4) 



glued to the rightmost position at the bottom of X. 



by cyclically moving one step forward those letters occurring after the image 
of the last letter in the bottom line, i.e., after the letter n. In the second case, 
when A n < A 7r -i(„), the new permutation it' corresponds to the modification of 
the top ordering by cyclically moving one step forward those letters occurring 
after the image of the last letter in the top line, i.e., after the letter 7r _1 (n). 

Example. For the initial buildings at both Figures 24 and 25 the permutation 
7r corresponding to the initial interval exchange transformation T : X — * X is 
the same and equals 

'\ 2 3 4 N 
^43 2 1, 

Our modification produces permutation 



in the first case (when A„ > A w - 



12 3 4 



lj J \4 1 3 2 y 
l(„)) and permutation 
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Fig. 25. Type II modification: the rightmost rectangle R<t on top of X is narrower 
than the rectangle R\ — R^-iu) glued to the rightmost position at the bottom of 
X. 



(\ ^2^3^4^\/l 42 3\ A 2 3 4\ _ , 
^4 3 2 1 /l v 4 32iyl v 2 431/ 7r 

in the second case (when A rl < A T -i(„)). 

Note that in the second case (when A„ < A„—i( n )) passing to the new 
decomposition X { U • • ■ U X' n we have to change the initial enumeration of 
the subintervals though physically all subintervals but one stay unchanged. 
Another choice would be to assign "names" to subintervals once and forever. 
Under the first choice the permutations 

(\ 4 2 3\ , (\ 2 3 4\ 
(,432 lJ and (,2 4 3 1 J 

coincide; under the latter choice they become different permutations. The 
article [Y] in the current volume adopts the second convention. 

Similarly to the case of interval exchange transformations of two intervals 
the induction procedure described above can be described entirely in terms of 
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interval exchange transformations T : X — > X and T' : X' — > X'. Historically 
it was proposed by G. Razy [Ra] in these latter form and then was interpreted 
by W. Veech [Vc3] in terms of zippered rectangles. (Actually, the zippcrcd 
rectangles decomposition has appeared and was first studied in [Vc3].) 



5.6 Multiplicative Cocycle on the Space of Interval Exchanges 

The renormalization procedure constructed in Sec. 5.4 gives a map g from the 
space of rotations of a circle to itself, or, in other terms, from the space of 
interval exchange transformations of two subintervals to itself. The permuta- 
tion 7r corresponding to an interval exchange transformations of two intervals 
X\ UX2 is always equal to 7r = (2, 1), so such interval exchange transformation 
can be parametrized by a single real parameter A € (0,1), where A = \X\\. 
Here we assume that the total length \X\ \ + \X 2 \ = \X\ of the interval X is 
normalized as \X\ = 1. 

An interval exchange transformation of n subintervals X = X\ U ■ • • U X n 
is parametrized by a collection Ai, . . . , A„ of positive numbers representing 
the lengths of subintervals and by a permutation 7r £ & n . Assuming that the 
total length of the interval X is normalized as \X\ = 1 we see that the space 
of interval exchange transformations is parametrized by a finite collection of 
(n— l)-dimensional simplices A"^ 1 = {(Ai,...,A„) | Ai + - • -+A„ = l;Aj > 0}, 
where each simplex corresponds to some fixed permutation it. 

As a collection of permutations one can consider all permutations obtained 
from a given one by applying recursively the modifications described at the end 
of the previous section. Such collection of permutations is called a Rauzy class 
9t C 6 n . Figure 26 illustrates the Rauzy class of the permutation (4, 3, 2, 1), 
where the arrows indicated modifications of the first and of the second type. 



n(^(324 1 ! 




(2431) 



(2413) 



1 (4321) 11 > (4132) 11 




(3142) 



I II 

Fig. 26. Rauzy class of permutation (4, 3, 2, 1) 



The renormalization procedure described in the previous section (com- 
bined with rescaling of the resulting interval X' to the unit length) defines a 
map 

T : A n - X x m -> A 11 ' 1 x m 
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for each space of interval exchange transformations to itself. The following 
Theorem of W. A. Veech (see [Ve3]) is crucial in this story. 

Key Theorem (W. A. Veech). The map T is ergodic with respect to ab- 
solutely continuous invariant measure. 

Remark. It is easy to see that the Rauzy class presented at Fig. 26 does 
not depend on the starting permutation. Actually, the same is true for any 
Rauzy class. Moreover, for almost any flat surface S € H. comp (di, . . . , d m ) 
in any connected component of any stratum a finite set of all permutations 
realizable by hrst return maps to all possible horizontal segments X satisfying 
Convention 2 does not depend on the surface S. This set is a disjoint union 
of a finite collection of corresponding Rauzy classes 9\\ , . . . , ?Rj , where j is 
the number of distinct entries di 1 < di 2 < ■ ■ ■ < di j . The union Ul ex = £Hi U 
• • • U y\j is called the extended Rauzy class; it depends only on the connected 
component H comp (di, . . . , dm)- In particular, connected components of the 
strata arc characterized by the extended Rauzy classes, where the latter ones 
can be described in purely combinatorial terms. 

With the Theorem above we have almost accomplished our scheme for a 
renormalization procedure. There is only one trouble with the map T: the 
measure mentioned in the Theorem is infinite (the total measure of the space 
of interval exchange transformations is infinite). This technical problem can 
be fixed by the following trick. We shall modify the renormalization algorithm 
described in the previous section by making several modifications of the zip- 
pered rectangle at a time. At a single step of the new algorithm Q wc apply 
several steps of the previous one T. Namely, we keep going as soon as we ap- 
ply consecutive transformations T of the same type / or of the same type II. 
Conceptually it does not change the renormalization procedure, but now the 
renormalization develops faster than before. The following example illustrates 
the correspondence between the renormalization procedures T and Q : 



(A, tt) - T(A, tt) ^ T 2 (A, tt) ^ T 3 (A, tt) ^ T 4 (A, tt) ^ T 5 (A, tt) ^ 



The accelerated procedure Q was introduced in [Zo2] , where the following 
Theorem was proved. 

Theorem. The map Q is ergodic with respect to absolutely continuous invari- 
ant probability measure on each space of zippered rectangles. 

Now, when wc have elaborated almost all necessary tools we are ready to 
give an idea of the proof of the Theorem from Sec. 4.3 concerning asymptotic 
flag. The last element which is missing is some analysis of the matrices B(X, tt) 



II 

(A,7T) 



S(A,7T) 



II 



in (12). 
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Multiplicative Cocycle 

In our interpretation of a renormalization procedure as a map on the space of 
interval exchange transformations we have to consider matrix B as a matrix- 
valued function B(X, tt) on the space A 11 ^ 1 x 9^ of interval exchange transfor- 
mations. Our goal (as it was outlined in Sec. 5.3) is to describe the properties 
of the products B^ = B{X^)-B{X^) B(X^) of matrices correspond- 
ing to successive steps of renormalization. 

We shall keep the same notation B(X,w) for matrices corresponding to 
our fast renormalization procedure Q(X,tt). Consider the product B^(X, tt) 
of values of B(-) taken along the orbit (A, tt), G(X, tt), . . . ,Q s (X, tt) of length s 
of the map Q : 

B (s >(A,7r) = B(X,tt) ■ B(G(X,tt)) B(£ s (A,tt)) 

Such B^(X, 7r) is called a multiplicative cocycle over the map Q: 

B { P +q \X,Tr) = B ( p\X,tt) ■ B iq \g p {X,Tr)) 

Using ergodicity of Q we can apply multiplicative ergodic theorem (see Ap- 
pendix B) to describe properties of B^ s \ Morally, the multiplicative ergodic 
theorem tells that for large values of s the matrix B^(X,tt) should be con- 
sidered as a matrix conjugate to the s-th power of some constant matrix. (See 
Appendix B for a rigorous formulation.) The logarithms of eigenvalues of this 
constant matrix are called Lyapunov exponents of the multiplicative cocycle. 

Recall that matrix B T (X, tt) was defined as a matrix representing the new 
"first return cycles" in terms of the old ones, see (12). Actually, it can be also 
interpreted as a matrix representing a change of a basis in the first relative co- 
homology [Re a;] = (Ai, . . . , A„) £ H 1 (S, {conical singularities}; R). It is easy 
to check that it respects the (degenerate) symplectic form: the intersection 
form (15). Note that symplectic matrices have certain symmetry of eigenval- 
ues. In particular, it follows from the general theory that the corresponding 
Lyapunov exponents have the following symmetry: 

0i > 2 > 6> 3 > ■ • ■ > g > = Q = ■ ■ ■ = Q > -6 g > ■ ■ ■ > -0 2 > -6 l 

number of conical points —1 

Note that the first return cycles actually belong to the absolute homology 
group Hx(S;Z) c Hi(S;M.) ~ M. 2g . Passing to this 2g-dimensional space we 
get matrices which already preserve a nondegenerate symplectic form. They 
define the following subcollection 

Oi > 2 > 3 > ■ ■ ■ > 6 g > -6 g > ■ ■ ■ > -9 3 > -9 2 > -9 X 

of Lyapunov exponents. 

The rest is an elementary linear algebra. We want to describe how do the 
large powers s of a symplectic matrix with eigenvalues 
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exp(0i) > exp(6> 2 ) > • ■ ■ > exp(0 ff ) > exp(— S ) >■••> exp(-0 2 ) > cxp(-6' 1 ) 

act on a 2<7-dimensional symplcctic space. 

We know that the Lyapunov exponent 6 g > is nonnegativc. Assume 4 
that it is actually strictly positive: 8 g > 0. Then for half of dimensions our 
linear map is expanding and for half of dimensions it is contracting. In par- 
ticular, under the assumption that 9 g > we conclude that the linear map 
B( s > (A, 7r) projects all homology space to a Lagrangian subspace (spanned by 
eigenvectors corresponding to positive Lyapunov exponents). 

Assuming that the spectrum of Lyapunov exponents is simple, that is 
assuming that 

01 > 02 > 03 > • • • > 0g 

we get the entire picture of deviation. A generic vector in the homology space 
stretches along the principal eigenvector (the one corresponding to the eigen- 
value 9\) with a factor exp(s#i); it expands along the next eigenvector with 
a factor exp(s$2), etc, up to the order g; its deviation from the Lagrangian 
subspace spanned by the first g eigenvectors tends to zero. Hence, the norm I 
of the image of a generic vector under s-th power of our linear map is of the 
order / ~ exp(s0i); its deviation from the subspace Vi, which is spanned by 

the top eigenvector, is of the order exp(s02) = ; its deviation from the sub- 

83 

space V2 spanned by the two top eigenvectors is of the order exp(s03) = ZT , 
etc; there is no deviation from the Lagrangian subspace spanned by the top 
g eigenvectors. In particular, the exponent i/j responsible for deviation from 
the subspace Vj_i from the Theorem in Sec. 4.3 is obtained by normalization 
of the Lyapunov exponent 9j by the leading Lyapunov exponent 6\ : 



This completes the proof of the Theorem in the case when the vertical 
trajectory starts and ends at the same microscopic horizontal interval (in 
other words, when the piece of trajectory is "almost closed"). Applying some 
additional (relatively involved) ergodic machinery one can complete the proof 
of the Theorem for arbitrary long pieces of vertical trajectories; see [Zo4] for 
a complete proof. 

I would like to stress that the original Theorem proved in [Zo4] is condi- 
tional: the statement about Lagrangian subspace was proved modulo conjec- 
ture that 9 g > 0; the statement about a complete Lagrangian flag was proved 
modulo conjectural simplicity of the spectrum of Lyapunov exponents. 

Positivity 6 g > was proved by G. Forni in [Fori], and simplicity of the 
spectrum was recently proved by A. Avila and M. Viana [AvVi]; see more 
details in Sec. 5.8. As it was shown in [Zo2] the proof of the strict inequality 
6\ > 62 immediately follows from results of W. A. Veech. 

4 Actually, this assumption is a highly nontrivial Theorem of G. Forni [Fori]; see 
below, see also Sec. 5.8 
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Exercise. In analogy with what was done in Sec. 5.4 consider the Rauzy-Veech 
induction in the torus case applying it to interval exchange transformations 
of two subintervals. We have seen that in this case the space of interval ex- 
change transformations is just an interval (0, 1). Find an explicit formulae for 
the Rauzy-Veech renormalization map T : (0, 1) — ► (0, 1) and for the "fast" 
rcnormalization map Q : (0, 1) — > (0, 1). Explain why the invariant measure 
is infinite for the map T. Find a relation between Q and the Gauss map 

g : x I— > < — > . Let 

P 1 _ 1 

q s l 

rii H 



n 2 



be the s-th best rational approximation of a real number x £ (0, 1). In the 
torus case the spectrum of Lyapunov exponents reduces to a single pair B\ > 
—6\. Show that for almost all x £ (0, 1) the Lyapunov exponent d\ (called in 
number theory the Levy constant, see [Lv]) is responsible for the growth rate 
of the denominator of the continued fraction expansion of x: 



,. log<7s 
hm 



12 log 2 



5.7 Space of Zippered Rectangles and Teichmuller geodesic flow 

We have proved the Theorem about an asymptotic Lagrangian flag of sub- 
spaces responsible for deviation of the cycles c(x, N) from asymptotic direc- 
tion. We have also proved that the exponents Vj responsible for the quan- 
titative description of the deviation are expressed in terms of the Lyapunov 
exponents of the multiplicative cocycle corresponding to our rcnormalization 
procedure: Uj = dj/di- 

There remains a natural question why should we choose this particular 
renormalization procedure and not a different one. One more natural question 
is what is the relation between renormalization procedure and the flow induced 

by the action of the diagonal subgroup ^ cx pW ^ ^ \ on the space of flat 

surfaces which we agreed to call the Teichmuller geodesic flow; this relation 
was announced in Sec. 4.3. This section answers to these questions which are, 
actually, closely related. 

In our presentation we follow the fundamental paper [Ve3] of W. A. Veech; 
the material at the end of the section is based on the paper [Zo2] developing 
the initial paper [Ve3]. 
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Space of Zippered Rectangles 

We have seen that locally a flat surface S can be parametrized by a collection of 
relative periods of the holomorphic 1-form ui representing the flat surface, i.e. 
we can choose a small domain containing [uj] in the relative cohomology group 
H 1 (S, {Pi, . . . , P m };C) as a coordinate chart in the corresponding stratum 
H(di, . . . , d m ). 

Decomposition of a flat surface into zippered rectangles gives another sys- 
tem of local coordinates in the stratum. Namely, choose a horizontal segment 
X satisfying Convention 2 from Sec. 5.3 and consider the corresponding de- 
composition of S into zippered rectangles. Let Ai, . . . , A„ be the widths of the 
rectangles, hi, . . . , h n their heights, and ai, . . . , a m be the altitudes responsi- 
ble for the position of singularities (we zip the neighboring rectangles Rj and 
Rj+i from the bottom up to the altitude dj and then the rectangles split at 
the singularity, see Figures 21, 24, 25). There is one more parameter describ- 
ing a decomposition of a flat surface into zippered rectangles: a permutation 
7r G & n - This latter parameter is discrete. 

The vertical parameters hj, au and are not independent: they satisfy some 
linear equations and inequalities. Varying the continuous parameters A, h, a 
respecting the linear relations between parameters h and a we get a new 
set of coordinates in the stratum. These coordinates were introduced and 
studied by W. A. Veech in [Ve3]. In particular, it was proved that for any 
(A, 7r) G A 11 ' 1 x $K in the space of interval exchange transformations there 
is an n-dimcnsional open cone of solutions (h, a). In other words, having any 
interval exchange transformation T : X — > X one can always construct a flat 
surface S and a horizontal segment X C S inside it such that the first return 
of the vertical flow to X gives the initial interval exchange transformations. 
Moreover, there is a n-dimensional family of such flat surfaces - suspensions 
over the interval exchange transformation T : X — > X (see [Ma3], [Ve3]). 

Is there a canonical decomposition of a flat surface into zippered rectan- 
gles? A choice of horizontal segment X C S completely determines a decom- 
position of a generic surface S into zippered rectangles, so our question is 
equivalent to the problem of a canonical choice of a horizontal segment X 
satisfying Convention 2. The choice which we propose is almost canonical; it 
leaves an arbitrariness of finite order which is the same for almost all S in the 
stratum. Here is the choice. Let us place the left extremity of the horizontal 
segment X at one of the conical singularities, and let us choose the length 
\X\ of the segment in such way that X would be the shortest possible interval 
satisfying Convention 2 and condition \X\ > 1. 

In practice the interval X can be constructed as follows: start with a 
sufficiently long horizontal interval having its left extremity at a conical point 
and satisfying Convention 2. Apply the "slow" Rauzy- Veech algorithm as long 
as the resulting subinterval has length at least 1. For almost all flat surfaces 
after finite number of steps we obtain the desired interval X. 
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The surface S has finite number of conical singularities; each conical sin- 
gularity has finite number of horizontal prongs, so we get arbitrariness of 
finite order. Thus, the resulting space of zippered rectangles can be essentially 
viewed as a (ramified) covering over the corresponding connected component 
H comp (di , . . . , d m ) of the stratum. Passing to a codimcnsion one subspace fl 
defined by the condition A • h = 1 we get a space of zippered rectangles of area 
one covering the space y^° m% '(di, . . . , d m ) of flat surfaces of area one. Consider 
a codimension-one subspace T C fl of zippered rectangles which have unit 
area, and which have the base X of length one, \X\ = \\ + ■ ■ ■ + A„ = 1. The 
space T has a natural structure of a fiber bundle over the space A n ~ x x JH 
of interval exchange transformations: we associate to a zippered rectangle 
(A, h, a, 7r) € T the interval exchange transformation (A, tt). 

We would like to emphasize an interpretation of fl as a fundamental do- 
main in the space of all zippered rectangles of area one. As a fundamental 
domain fl is defined by the additional condition on the base: X is the shortest 
possible interval satisfying Convention 2 such that \X\ > 1. In this interpreta- 
tion Y is the boundary of the fundamental domain. Starting with an arbitrary 
zippered rectangle representation satisfying Convention 2 we can apply sev- 
eral steps of Rauzy-Veech algorithm (see Fig. 24 and Fig. 25), which does 
not change the surface S. After several iterations we get to the fundamental 
domain fl. 

We would use the same notation for fl considered as a fundamental domain 
and for fl considered as a quotient, when two boundary components of fl are 
identified by the modification of zippered rectangles as on Fig. 24 and Fig. 25. 

Teichmiiller Geodesic Flow and its First Return Map to a Cross-section 

Zippered rectangles coordinates are extremely convenient when working with 
the Teichmiiller geodesic flow, which we identify with the action of the di- 
agonal subgroup gt = ^ ex P^ ^ ^ \ . Namely, gt expands the horizontal 

parameters A by the factor exp(t) and contracts the vertical parameters h, a 
by the same factor. 

Consider a zippered rectangle S — (A, h, a, tt) G T with the base X of unit 
length. Applying g t to S with t continuously increasing from t — we shall 
eventually make the length of the base of the deformed zippered rectangle 
g t S = (exp(i)A, exp(— t)h, exp(— t)a; n) too long and thus we shall get outside 
of the fundamental domain fl. It is not difficult to determine an exact time 
to when it will happen. We get to the boundary of the fundamental domain 
fl at the time 

t (S) = -log(l-min(A„,A 7r -i (n) )). (17) 

The time to is chosen in such way that applying to the zippered rectangle 
g to S one step of the Rauzy-Veech induction (see Fig. 24 and Fig. 25) we 
get a new zippered rectangle with the base X' of unit length. To verify for- 
mula (17) for to it is sufficient to note that expansion-contraction commutes 
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with Rauzy-Veech induction. Thus, to evaluate t we can first apply one step 
of the Rauzy-Veech induction and then apply expansion-contraction for an 
appropriate time, which would bring us back to T, i.e. which would make the 
length of the base of the new building of zippered rectangles equal to one. 

In other words, starting at a point S £ T and following the flow for the 
time to(S) we get to the boundary of the fundamental domain in the space 
of zippered rectangles and we have to instantly jump back to the point of Y 
identified with gt S. One can recognize in this construction the first return 
map S : Y — > Y defined by the flow g t on the section Y: at the time to(S) 
the flow g t emitted from a point 5eT returns back to the codimension-one 
subspace Y transversal to the flow. 

Morally one should consider the map S as a map on some subspace of flat 
surfaces. Note, that S is not applicable to points of flat surfaces, it associates 
to a flat surface taken as a whole another flat surface taken as a whole. 

We see now that the Rauzy-Veech renormalization procedure S : Y — > Y 
performed on the level of zippered rectangles is nothing but discrete version 
of the Teichmullcr geodesic flow. Namely S is the first return map of the 
Teichmuller geodesic flow to a section Y. By construction the Rauzy-Veech 
induction T : A' 1-1 x EH — > Z\ n_1 x 5H on the space of interval exchange 
transformations is just a projection of S. In other words, the following diagram 

r(9t) — ^— » Y{m) 



A n ~ x x m — - — ► A n ~ x x eh 

is commutative, and the invariant measure on the space Z\™ _1 x EH of interval 
exchange transformations is a push forward of the natural invariant measure 
on the space Y of zippered rectangles. 

Choice of a Section 

Now we can return to the questions addressed at the beginning of this section. 
Ignoring an algorithmic aspect of the choice of renormalization procedure we 
see that conceptually, it is defined by a section of the Teichmuller geodesic 
flow. In particular, the "fast" renormalization procedure Q : A 71 ^ 1 x EH — ► 
A n ~ 1 x EH defined in the previous section corresponds to a choice of a sub- 
section Y' C Y. Luckily it has a simple algorithmic representation in terms of 
modification of the interval exchange transformation T(X, w), and, moreover, 
it has a simple description in terms of coordinates A, h, a, n in the space of 
zippered rectangles given by an extra condition for the parameter a n . 

Recall that parameters a,j are responsible for the position of singularities: 
we zip the neighboring rectangles Rj and Rj+i from the bottom up to the 
altitude a,j , see Fig. 21. In particular, by construction all a,j for j = 1, . . . , n— 1 
are positive. Parameter a n is, however, different from the others: the rectangle 
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R n is the rightmost rectangle in the collection. If there is a conical singular- 
ity located at the right side of this rightmost rectangle (sec, for example, 
the zippered rectangle decomposition of the flat surface on the top part of 
Fig. 24), then parameter a m is positive; it indicates as usual at what height 
is located the singularity. However, the right side of the rightmost rectangle 
might contain no singularity. This means that the singularity is located on 
the corresponding vertical trajectory below the zero level of the base X. The 
rectangle which is glued to X from below at the rightmost position is the 
rectangle R v -ir n y, the singularity is located on the right side of this rectangle 
(see, for example, the zippered rectangle decomposition of the flat surface on 
the bottom part of Fig. 24). In this case we let a n be negative indicating how 
low we have to descend along downward vertical trajectory emitted from the 
right endpoint of X to hit the singularity. 

The subsection Y' is defined by the following extra condition 

Y' ={ (A, h, a, it) G Y \ a n > when \ n > A ff -i(„)} U 
U { (A, h, a, it) £ T \ a n < when A„ < A„—i/ n )} 

Exercise. Check which zippered rectangles at Figures 21, 24, 25 satisfy the 
condition a n ■ (A„ — A^-i^) > and which do not. 

It can be verified (see [Zo2] ) that the section Y' is still a fiber bundle over 
the space of zippered rectangles and that the corresponding first return map 
S' : Y' — > Y' projects to the map Q: 

r'(m) r(st) 

I I (18) 

£n-l x jfl A n-\ x Sft 

Exercise. Verify that the definition of the renormalization procedure Q as 
a projection of the first return map of the Tcichmuller geodesic flow to Y' 
matches the intrinsic definition of Q given in Sec. 5.6. 

Different choices of the section also explain why the invariant measure 
on the space of interval exchange transformations A"" 1 x £H was infinite for 
the Rauzy-Veech induction T while is finite for the "fast" renormalization 
procedure Q. As a model case consider a directional flow on a torus and two 
different sections to this flow. Taking as a section the line Y represented on 
the left picture of Fig. 27 we get a section of infinite measure though the 
measure of the torus is finite and the flow is very nice. Taking as a section a 
finite piece Y' C Y as on the right side of Fig. 27 we get a section of finite 
measure. 

Similarly, the component 7i 1 omp (di, . . . , d m ) of the stratum has finite vol- 
ume and hence the space of zippered rectangles fi which is a finite covering 
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Fig. 27. The section Y on the left picture has infinite measure though the measure 
of the space is finite. The subsection Y' C Y on the right picture has finite measure. 
In both cases the first return map of the ergodic flow to the section is ergodic, but 
the mean return time to the left subsection is zero 

over Ti.1° mp (di, . . . ,d m ) also has finite volume. However, the initial section Y 
has infinite "hyperarea" while section Y' already has finite "hyperarea" . 

We complete this section with a several comments concerning Lyapunov 
exponents. Though these comments are too brief to give a comprehensive 
proof of the relation between exponents Uj responsible for the deviation and 
the Lyapunov exponents of the Teichmiiller geodesic flow, they present the 
key idea, which can be completed by an elementary calculation. 

It is clear that the Lyapunov exponents of the Teichmiiller geodesic flow 
g t on the stratum 7iJ omp (di, . . . , d m ) coincide with the Lyapunov exponents 
on its finite covering £1. The Lyapunov exponents of the flow g t differ from 
the Lyapunov exponents of the first return map S 1 to the section Y' only by a 
scaling factor representing the "hyperarea" of Y'. The map Q on the space of 
zippered rectangles coincides with the restriction of the map S' on the zippered 
rectangles restricted to horizontal parameters A and discrete parameter n. 
Thus, the Lyapunov exponents of Q form a subcollection of the Lyapunov 
exponents of S' corresponding to the subspace of horizontal parameters A. It 
remains to note that the Lyapunov exponents of the map Q are related to 
the Lyapunov exponents of the cocycle B(X, it) just by the scaling factor, and 
that we have already expressed the exponents Vj responsible for the deviation 
in terms of the Lyapunov exponents of the multiplicative cocycle B, see (16). 
Matching all the elements of this chain together we get a representation of the 
exponents Vj in terms of the Lyapunov exponents of the Teichmiiller geodesic 
flow given in Sec. 4.3. 

Zippered rectangles and Lyapunov exponents: more serious reading. More de- 
tails on Rauzy classes 9^, zippered rectangles, Lyapunov exponents of the 
Teichmiiller geodesic flow and their relation might be found in original papers 
of G. Rauzy [Ra], W. Veech [Ve3], [Ve6] and the author [Zo2], [Zo4]. 

5.8 Spectrum of Lyapunov Exponents (after M. Kontsevich, 
G. Forni, A. Avila and M. Viana) 

It should be mentioned that the statement that the subspace V g , such that 
\dist{c{x, N), V g )\ < const for any TV, has dimension exactly g was formulated 
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in the original paper [Zo4] as a conditional statement. It was based on the 
conjecture that the Lyapunov exponent v g is strictly positive. This conjecture 
was later proved by G. Forni in [Fori]. 

Theorem ( G. Forni). For any connected component Tt comp (di, . . . , d m ) of 
any stratum of Abelian differentials the first g Lyapunov exponents of the 
Teichmuller geodesic flow are strictly greater than 1: 

l + v g >\ 

As an indication why this positivity is not something which should be 
taken for granted we would like to give some precision about related results 
of G. Forni. 

The Lyapunov exponents of the Teichmuller geodesic flow play the role 
of (logarithms of) eigenvalues of a virtual "average monodromy of the tan- 
gent bundle along the flow". Instead of considering the tangent bundle to 
H(di , . . . , d m ) one can consider another vector bundle intimately related 
to the tangent bundle. This vector bundle has the space H l (S;M.) as a 
fiber. Since we know how to identify the lattices H 1 (S]X) C iJ 1 (S , ;M) and 
H X (S']Z) C H X {S'\^) in the fibers over two flat surfaces S and S' which 
are close to each other in H(di, . . . , d m ), we know how to transport the fiber 
if^SjK) over the "point" S in the base . . . , dm) to the fiber iJ^S^R) 

over the "point" S'. In other words, we have a canonical connection (called 
Gauss-Manin connection) in the vector bundle. Hence we can again study the 
"average monodromy of the fiber along the flow" . It is not difficult to show that 
the corresponding Lyapunov exponents are related to the Lyapunov exponents 
of the Teichmuller flow. Namely, the new collection of Lyapunov exponents 
has the form: 

1 > v-i > ■ ■ ■ > v g > —v g > ■ • ■ > —V2 > —I 

In particular, the collection of Lyapunov exponents of the Teichmuller geodesic 
flow can be obtained as follows: take two copies of the collection above; add 
+1 to all the entries in one copy; add —I to all entries in another copy; take 
the union of the resulting collections. The theorem of G. Forni tells that for 
any connected component of any stratum we have v g > 0. 

Consider now some SL(2, R)-invariant subvaricty Af C 7i(di, . . . , d m )- 
Consider the restriction of the vector bundle with the fiber H 1 (S; R) to N . We 
can compute the "average monodromy of the fiber along the Teichmuller flow" 
restricted to A/. It gives a new collection of Lyapunov exponents. Since the 
"holonomy" preserves the natural symplectic form in the fiber, the collection 
will be again symmetric: 

1 > ^ > ■ ■ ■ > v' g > -v' g > ■ ■ ■ > -v' % > -1 

G. Forni has showed [For2] that there are examples of invariant subvarieties Af 
such that all v'j, j = 2, ... ,g, are equal to zero! Moreover, G. Forni explicitly 



Flat Surfaces 65 



describes the locus where the monodromy does not change the fiber (or it 
exterior powers) too much, and where one may get multiplicities of Lyapunov 
exponents. 

Another conditional statement in the original paper [Zo4] concerns strict 
inclusions Vi C V2 C ■ • • C V g C H±(S; K). It was based on the other conjec- 
ture claiming that the Lyapunov exponents have simple spectrum. The first 
strict inequality v\ > V2 is an elementary corollary of general results of Veech; 
see [Zo3]. The other strict inequalities are much more difficult to prove. Very 
recently A. Avila and M. Viana [AvVi] have announced a proof of simplicity 
of the spectrum (19) for any connected component of any stratum proving 
the conjecture which was open for a decade. 

Theorem (A. Avila, M. Viana). For any connected component of any stra- 
tum H. cornp (di, . . . , d m ) of Abelian differentials the first g Lyapunov exponents 
are distinct: 

l + l/ 1 >l + l/ 2 >--->l+Vg (19) 

Sum of the Lyapunov exponents 

Currently there are no methods of calculation of Lyapunov exponents for 
general dynamical systems. The Teichmuller geodesic flow does not make an 
exception: there is some knowledge of approximate values of the numbers Vj 
obtained by computer simulations for numerous low-dimensional strata, but 
there is no approach leading to explicit evaluation of these numbers with 
exception for some very special cases. 

Nevertheless, for any connected component of any stratum (and, more 
generally, for any GL + (2; R)-invariant suborbifold) it is possible to evaluate 
the sum of the Lyapunov exponents v\ + . . . v g , where g is the genus. The 
formula for this sum was discovered by M. Kontsevich in [Kon]; it is given in 
terms of the following natural structures on the strata 7i(di, . . . , d m ). 

There is a natural action of C* on every stratum of the moduli space of 
holomorphic 1-forms: we can multiply a holomorphic form lj by a complex 
number. Let us denote by H(2) (di, ■ ■ ■ , d m ) the quotient of %(di, . . . , d m ) over 
C*. The space Tt(2){di, ■ ■ ■ , d m ) can be viewed as the space of flat surfaces of 
unit area without choice of distinguished direction. 

There are two natural holomorphic vector bundles over 7i( 2 )(<ii, . . . , d m )- 
The first one is the C*-bundle H(di, . . . , d m ) —* 7i(2)(di, . . . , dm). The second 
one is the C s -bundle, which fiber is composed of all holomorphic 1-forms in the 
complex structure corresponding to a flat surface S € W(2)(^i) ■ ■ ■ , d m ). Both 
bundles have natural curvatures; we denote by 71 and 72 the corresponding 
closed curvature 2-forms. 

Finally, there is a natural closed codimension two form (3 on every stra- 
tum 7i(2){d\, . . . , dm). To construct [3 consider the natural volume form Q 
on Ti.(di, . . . , dm). Four generators of the Lie algebra g[(2;R) define four 
distinguished vectors in the tangent space TsH{dx, . . . , d m ) at any "point" 
S G Ti(di, . . . , d m )- Plugging these four vectors in the first four arguments of 
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the volume form Q we get a closed codimension four form on Ti(di, . . . , d m ). It 
is easy to check that this form can be pushed forward along the C*-fibers of the 
bundle 7i{dx, . . . , d rn ) — > 7i(2)(di, ■ ■ ■ , d m ) resulting in the closed codimension 
two form on the base of this fiber bundle. 

Theorem (M. Kontsevich). For any connected component of any stratum 
the sum of the first g Lyapunov exponents can be expressed as 

J/3A 72 

V\ H V Vg = -j— , 

/ A 71 

where the integration is performed over the corresponding connected compo- 
nent ofH(2)(d\, . . . ,rf m ). 

As it was shown by G. Forni, this formula can be generalized for other 
GL + (2; R)-invariant submanifolds. 

The proof is based on two observations. The first one generalizes the fact 
that dynamics of the geodesic flow on the hyperbolic plane is in some sense 
equivalent to dynamics of random walk. One can replace Teichmiillcr geodesies 
by geodesic broken lines consisting of geodesic segments of unit length. Having 
a broken line containing n geodesic segments with the endpoint at the point 
S n we emit from S n a new geodesic in a random direction and stop at the 
distance one from S n at the new point S n +i. This generalization suggested 
by M. Kontsevich was formalized and justified by G. Forni. 

Consider the vector bundle over the moduli space of holomorphic 1-forms 
with the fiber H 1 (S; R) over the "point" S. We are interested in the sum i/i + 
• ■ ■ + v g of Lyapunov exponents representing mean monodromy of this vector 
bundle along random walk. It follows from standard arguments concerning 
Lyapunov exponents that this sum corresponds to the top Lyapunov exponent 
of the exterior power of order g of the initial vector bundle. In other words, 
we want to measure the average growth rate of the norm of a g-dimcnsional 
subspace in H 1 {S\M.) when we transport it along trajectories of the random 
walk using the Gauss-Manin connection. 

Fix a Lagrangian subspace L in the fiber H 1 (So;WL) over a "point" S n . 
Consider the set of points located at the Tcichmullcr distance 1 from S n . 
Transport L to each point S n +i of this "unit sphere" along the corresponding 
geodesic segment joining S n with S n +i) measure the logarithm of the change 
of the norm of L; take the average over the "unit sphere" . The key observation 
of M. Kontsevich in [Kon] is that for an appropriate choice of the norm this 
average growth rate is the same for all Lagrangian subspaces L in H 1 (Sq; R) 
and depends only on the point S n . A calculation based on this observation 
gives the formula above. 

Actually, formula above can be rewritten in a much more explicit form 
(which is a work in progress). The values of the sum given by this more 
explicit formula perfectly match numerical simulations. The table below gives 
the values of the sums of Lyapunov exponents for some low-dimensional strata; 
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this computation uses the results of A. Eskin and A. Okounkov [EOk] for the 
volumes of the strata. 



Conjectural values of v\ + • • • + v g for some strata 



W(2) 


H(l,l) 




W(4,M) 




1,1, 1,1,1) 


H(l, 1,1, 1,1, 1,1,1) 


4 
3 


3 
2 




1137 




839 
377 


235 761 


550 


93428 



In particular, since v\ = \ this information gives the exact value of the 
only nontrivial Lyapunov exponent V2 for the strata in genus two. Some extra 
arguments show that V2 = 1/3 for the stratum 7i(2) and for any GL + (2; R)- 
invariant submanifold in it; V2 = 1/2 for the stratum TC(1,1) and for any 
GL + (2; i?)-invariant submanifold in it. 



5.9 Encoding a Continued Fraction by a Cutting Sequence of a 
Geodesic 

We have seen that renormalization for a rotation of a circle (or equivalcntly 
for an interval exchange transformation of two subintervals) leads to the Eu- 
clidean algorithm which can be considered in this guise as a particular case 
of the fast Rauzy-Veech induction. 
The multiplicative cocycle 

^=G"')-a;)-"G"f)'U0- 

considered in section 5.6 corresponds to the decomposition of a real number 
x € (0, 1) into continued fraction, x = [0; n\, ri2, ■ ■ •], 

1 

x = 

1 



112 H 

A flat surface which realizes an interval exchange transformation of two 
subintervals is a fiat torus. The the moduli space of flat tori can be naturally 
identified with SL(2, R) / SL(2, Z) which in its turn can be naturally identified 
with the unit tangent bundle to the modular surface H 2 /SX(2, Z) see Sec. 3.2. 
Moreover, the Teichmuller metric on the space of tori coincides with the hy- 
perbolic metric on H 2 , and the Teichmuller geodesic flow on the moduli space 
of flat tori coincides with the geodesic flow on the modular surface. 



68 Anton Zorich 



Hence, the construction from the previous section suggests that the Eu- 
clidean algorithm corresponds to the following geometric procedure. There 
should be a section T in the (covering of) the unit tangent bundle to the 
modular surface and its subsection T' C Y such that the trajectory of the 
geodesic flow emitted from a point of Y' returns to T' after n\ intersections 
with Y, then after ni intersections with Y, etc. In other words there is a nat- 
ural way to code a continued fraction by a sequence of intersections (so called 
"cutting sequence" ) of the corresponding geodesic with some sections Y 1 C Y. 

Actually, a geometric coding of a continued fraction by a cutting sequence 
of a geodesies on a surface is known since the works of J. Nielsen and E. Art in 
in 20s and 30s. The study of the geometric coding was developed in the 80s 
and 90s by C. Series, R. Adler, L. Flatto and other authors. We refer to the 
expository paper [Ser] of C. Series for detailed description of the following 
geometric coding algorithm. 

Consider a tiling of the upper half plane with isometric hyperbolic triangles 
as at Fig. 29. A fundamental domain of the tiling is a triangle with vertices 
at 0,1 and oo; the corresponding quotient surface is a triple cover over the 
standard modular surface (see Fig. 47). This triangulation of H 2 by ideal 
triangles is also known as Farey tessellation. 

Consider a real number x £ (0, 1). Consider any geodesic 7 landing to the 
real axis at x such that 7 intersects with the imaginary axis; let iy be the 
point of intersection. Let us follow the geodesic 7 starting from iy in direction 
of x. Each time when we cross a triangle of our tiling let us note by the symbol 
L the situation when we have a single vertex on the left and two vertices on 
the right (see Fig 28) and by the symbol R the symmetric situation. 

single vertex single vertex 




Fig. 28. Coding rule: when we cross a triangle leaving one vertex on the left and 
two on the right we write symbol L; when there is one vertex on the right and two 
on the left we write symbol R 

Example. Following the geodesic 7 presented at Fig. 29 from some iy to x = 
( V85 - 5)/10) w 0.421954 we get a sequence R,R,L,L,R,L,L,R,R,L,... 
which we abbreviate as R 2 L 2 R 1 L 2 R 2 L 1 .... 

Theorem (C. Series). Let x € (0,1) be irrational. Let"/ be a geodesic emit- 
ted from some iy and landing at x; let R ni L n2 R n3 L" 4 . . . be the corresponding 
cutting sequence. Then x = [0; n\, ri2, «3, 714, . . . ] . 
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x 1 



Fig. 29. (After C. Series.) The cutting sequence defined by this geodesic starts with 
R, R, L, L, R, L, L, R,R,L,... which we abbreviate as R 2 L 2 R 1 L 2 R 2 L 1 .... The real 
number x £ (0, f ) at which lands the geodesic has continued fraction expansion 
a;=[0;2,2,l ) 2,2,l,...] 

Geometric symbolic coding: more serious reading. I can strongly recommend 
a paper of P. Arnonx [Arn] which clearly and rigorously explains the idea 
of suspension in the spirit of diagram (18) handling the particular case of 
Euclidean algorithm and of geodesic flow on the Poincare upper half-plane. As 
a survey on geometric coding we can recommend the survey of C. Series [Scr] 
(as well as other surveys in this collection on related subjects). 

6 Closed Geodesies and Saddle Connections on Flat 
Surfaces 

Emitting a geodesic in an irrational direction on a flat torus we get an irra- 
tional winding line; emitting it in a rational direction we get a closed geodesic. 
Similarly, for a flat surface of higher genus a countable dense set of directions 
corresponds to closed geodesies. 

In this section we study how many closed regular geodesies of bounded 
length live on a generic flat surface S. We consider also saddle connections 
(i.e. geodesic segments joining pairs of conical singularities) and count them. 

We explain a curious phenomenon concerning saddle connections and 
closed geodesies on flat surfaces: they often appear in pairs, triples, etc of 
parallel saddle connections (correspondingly closed geodesies) of equal length. 

When all saddle connections (closed geodesies) in such configuration be- 
come short the corresponding flat surface starts to degenerate and gets close 
to the boundary of the moduli space. Thus, a description of possible configu- 
rations of parallel saddle connections (closed geodesies) gives us a description 
of the multidimensional "cusps" of the strata. 
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6.1 Counting Closed Geodesies and Saddle Connections 

Closed Regular Geodesies Versus Irrational Winding Lines 

Consider a flat torus obtained by identifying pairs of opposite sides of a unit 
square. A geodesic emitted in an irrational direction (one with irrational slope) 
is an irrational winding line; it is dense in the torus. A geodesic emitted in a 
rational direction is closed; all parallel geodesic are also closed, so directional 
flow in a rational direction fills the torus with parallel periodic trajectories. 
The set of rational directions has measure zero in the set of all possible direc- 
tions. In this sense directions representing irrational winding lines are typical 
and directions representing closed geodesic - nontypical. 

The situation with flat surfaces of higher genera g > 2 is similar in many 
aspects, though more complicated in details. For example , for any flat surfaces 
S almost all directions are "irrational" ; any geodesic emitted in an irrational 
direction is dense in the surface. Actually, even stronger statement is true: 

Theorem (S. Kerckhoff, H. Masur, J. Smillie). For any flat surface 
directional flow in almost any direction is uniquely ergodic. 

For the torus the condition that directional flow is minimal (that is any 
trajectory going in this direction is dense in the torus) is equivalent to the 
condition that the flow is uniquely ergodic (the natural Lebesgue measure 
induced by the flat structure is the only finite measure invariant under direc- 
tional flow; see Appendix A for details). Surprisingly a directional flow on a 
surface of higher genus (already for g = 2) might be minimal but not uniquely 
ergodicl Namely, for some directions which give rise to a minimal directional 
flow it might be possible to divide the surface into two parts (of nonzero mea- 
sure) in such way that some trajectories would mostly stay in one part while 
other trajectories would mostly stay in the other. 

Closed geodesies on flat surfaces of higher genera also have some similar- 
ities with ones on the torus. Suppose that we have a regular closed geodesic 
passing through a point xq E S. Emitting a geodesic from a nearby point x in 
the same direction we obtain a parallel closed geodesic of the same length as 
the initial one. Thus, closed geodesies also appear in families of parallel closed 
geodesies. However, in the torus case every such family fills the entire torus 
while a family of parallel regular closed geodesies on a flat surfaces of higher 
genus fills only part of the surface. Namely, it fills a flat cylinder having a 
conical singularity on each of its boundaries. 

Exercise. Find several periodic directions on the flat surface from Fig. 12. 
Find corresponding families of parallel closed geodesies. Verify that each of 
the surfaces from Fig. 44 decomposes under the vertical flow into three cylin- 
ders (of different circumference) filled with periodic trajectories. Find these 
cylinders. 
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Counting Problem 

Take an arbitrary loop on a torus. Imagine that it is made from a stretched 
elastic cord. Letting it contract we get a closed regular geodesic (may be 
winding several times along itself). Now repeat the experiment with a more 
complicated flat surface. If the initial loop was very simple (or if we are ex- 
tremely lucky) we again obtain a regular closed geodesic. However, in general 
we obtain a closed broken line of geodesic segments with vertices at a collection 
of conical points. 

Similarly letting contract an elastic cord joining a pair of conical singular- 
ities we usually obtain a broken line composed from several geodesic segments 
joining conical singularities. In this sense torus is very different from a general 
flat surface. 

A geodesic segment joining two conical singularities and having no conical 
points in its interior is called saddle connection. The case when boundaries of 
a saddle connection coincide is not excluded: a saddle connection might join 
a conical point to itself. 

Convention 4- In this paper we consider only saddle connections and closed 
regular geodesies. We never consider broken lines formed by several geodesic 
segments. 

Now we are ready to formulate the Counting Problem. Everywhere in this 
section we normalize the area of flat surfaces to one. 

Counting Problem. Fix a flat surface S. Let N SC (S, L) be the number of 
saddle connections on S of length at most L. Let N cg (S,L) be the number of 
maximal cylinders filled with closed regular geodesies of length at most L on 
S. Find asymptotics of N SC (S, L) and N cg (S,L) as L — > oo. 

It was proved by H. Masur (see [Ma5] and [Ma6]) that for any flat surface 
S counting functions N(S, L) grow quadratically in L. Namely, there exist 
constants < consti(S) < const2(S) < oo such that 



for L sufficiently large. Recently Ya. Vorobets has obtained in [Vb2] uniform 
estimates for the constants consti(S) and const2(S) which depend only on 
the genus of S. 

Passing from all flat surfaces to almost all surfaces in a given connected 
component of a given stratum one gets a much more precise result; see [EMa]. 

Theorem (A. Eskin and H. Masur). For almost all flat surfaces S in a 
given connected component of a stratum 7i(di, . . . , d m ) the counting functions 
N SC (S,L) and N cg (S,L) have exact quadratic asymptotics 



consh(S) < N(S,L)/L 2 < const 2 (S) 



lim 



N SC {S, L) 
ttL 2 



const 



sc 



lim 



N cg (S,L) 
irL 2 



= const, 



(20) 
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where Siegel-Veech constants const sc and const cg are the same for almost all 
flat surfaces in the component Ti^ {d\, . . . , d m ). 

We multiply denominator by it to follow a conventional normalization. 

Phenomenon of Higher Multiplicities 

Let us discuss now the following problem. Suppose that we have a regular 
closed geodesic on a flat surface S. Memorize its direction, say, let it be the 
North- West direction. (Recall that by Convention 1 in Sec. 1.2 we can place a 
compass at any point of the surface and it will tell us what is the direction to 
the North.) Consider the maximal cylinder filled with closed regular geodesies 
parallel to ours. Take a point x outside this cylinder and emit a geodesic from 
x in the North- West direction. There are two questions. 

- How big is the chance to get a closed geodesic? 

- How big is the chance to get a closed geodesic of the same length as the 
initial one? 

Intuitively it is clear that the answer to the first question is: "the chances 
are low" and to the second one "the chances are even lower" . This makes the 
following Theorem (see [EMaZo]) somehow counterintuitive: 

Theorem (A. Eskin, H. Masur, A. Zorich). For almost all flat sur- 
faces S from any stratum different from H.i(2g — 2) or 7ii(di,d2) the func- 
tion N two cyi(S,L) counting the number of families of parallel regular closed 
geodesies filling two distinct maximal cylinders has exact quadratic asymp- 
totics 

.. Nfwo_cyl(S, L) 

lim — = consttwo cyl 

where Siegel-Veech constants const two C yl > depends only on the connected 
component of the stratum. 

For almost all flat surface S in any stratum one cannot find a single pair 
of parallel regular closed geodesies on S of different length. 

There is general formula for the Siegel Vecch constant const two cyl an( l f° r 
similar constants which gives explicit numerical answers for all strata in low 
genera. Recall that the principal stratum H(l, . . . , 1) is the only stratum of fill 
dimension in Tt g ; it is the stratum of holomorphic 1-forms with simple zeros 
(or, what is the same, of flat surfaces with conical angles 4ir at all cone points). 
Numerical values of the Siegel-Veech constants for the principal stratum are 
presented in Table 2. 

Comparing these values we see, that our intuition was not quite misleading. 
Morally, in genus g = 4 a closed regular geodesic belongs to a one-cylinder 
family with "probability" 97.1%, to a two-cylinder family with "probabil- 
ity" 2.8% and to a three-cylinder family with "probability" only 0.1% (where 
"probabilities" are calculated proportionally to Siegel-Veech constants). 
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9 — 1 


9 — 1 


3-3 




single 
cylinder 




5 1 

2 wr 1 - 52 


36 1 


3150 1 
377 'as)" 5 - 08 


two 
cylinders 






3 1 

uwr 0A3 


90 1 

377- CM* - 145 


three 
cylinders 











Table 2. Siegel-Veech constants const n _ cy i for the principal stratum Tti(l, •••,!) 



In theorem above we discussed closed regular geodesies. A similar phe- 
nomenon is true for saddle connections. Recall that the cone angle at a conical 
point on a flat surface is an integer multiple of 2ir. Thus, at a point with a cone 
angle 2im every direction is presented n times. Suppose that we have found 
a saddle connection of length going from conical point Pi to conical point 
P2. Memorize its direction (say, the North- West direction) and its length I. 
Then with a "nonzero probability" (understood in the same sense as above) 
emitting a geodesic from Pi in one of the remaining n — 1 North- West direc- 
tions we make it hit Pi at the distance I. More rigorously, the Siegel-Veech 
constant counting configurations of two parallel saddle connections of equal 
length joining Pi to P2 is nonzero. 

The explicit formula for any Siegel-Veech constant from [EMaZo] can be 
morally described as the follows. Up to some combinatorial factor responsible 
for dimensions, multiplicities of zeroes and possible symmetries any Siegel- 
Veech constant can be obtained as a limit 

1 Vol( "e-neighborhood of the cusp C ") 
c(C) = lim — - comp — (21) 

e^O -KE Z Volrt-L (di, . . . ,d m ) 

where C is a particular configuration of saddle connections or closed geodesies. 

Say, as a configuration C one can consider a configuration of two maximal 
cylinders filled with parallel closed regular geodesies of equal lengths. The 
e-neighborhood of the corresponding cusp is the subset of those flat surfaces 
S G Tii° mp (di 1 . . . ,d m ) which have at least one pair of cylinders filled with 
parallel closed geodesies of length shorter than e. 

As another example one can consider a configuration of three parallel 
saddle connections of equal lengths on S E 7ii(l,l,4, 8) joining zero Pi of 
degree 4 (having cone angle lOV) to zero P2 of degree 8 (having cone angle 
187r) separated by angles 2ir, 2ir, 6ir at Pi and by angles 6ir, lOir, 2ir at P2. The 
e-neighborhood H\(l, 1, 4, 8) C Hi(l, 1, 4, 8) of the corresponding cusp is the 
subset of those flat surfaces in Hi(l, 1, 4, 8) which have at least one triple of 
saddle connections as described above of length shorter than e. 
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We explain the origin of the key formula (21) in the next section. In sec- 
tion 6.4 we give an explanation of appearance of higher multiplicities. 

Other counting problems (after Ya. Vorobets) 

Having a flat surface S of unit area we have studied above the number of max- 
imal cylinders N cg (S, L) filled with closed regular geodesies of length at most 
L on S . (In this setting when we get in some direction several parallel maxi- 
mal cylinders of equal perimeter, we count each of them.) In the paper [Vb2] 
Ya. Vorobets considered other counting problems. 

In particular, among all maximal periodic cylinders of length at most L 
(as above) he counted the number N cg>a (S, L) of those ones, which have area 
greater than a. He also counted the total sum N area (S, L) of areas of all 
maximal cylinders of perimeter at most L and the number N X (S, L) of regular 
periodic geodesies of length at most L passing trough a given point i£S. 

Ya. Vorobets has also studied how the maximal cylinders filled with closed 
geodesies are distributed with respect to their direction and their area. He 
considered the induced families of probability measures on the circle S , on the 
unit interval [0, 1] and on their product S 1 x [0; 1]. Given a subset U C S\, V C 
[0; 1], W C S 1 x [0; 1] the corresponding measures dir^U), ar L (V),pair L (W) 
tell the proportion of cylinders of bounded perimeter having direction in U , 
area in V, or the pair (direction, area) in W correspondingly 

Using the general approach of A. Eskin and H. Masur, Ya. Vorobets has 
proved in [Vb2] existence of exact quadratic asymptotics for the counting 
functions introduced above. He has computed the corresponding Siegcl-Veech 
constants in terms of the Siegcl-Veech constant const cg in (20) and found the 
asymptotic distributions of directions and areas of the cylinders: 

Theorem (Ya. Vorobets). For almost any flat surface S of unit area in 
any connected component of any stratum W° (di, . . . , d rn ) and for almost 
any point x of S one has 

N cg , a (S,L) N area (S,L) N X (S,L) 
hm — = c cgi<J hm — = c area hm — — = c x , 

L->oo 7Tiv L^oo TrL z L^ao It 

\9 -vi- const cg 
where c ca a = (1 — a) 9 +m ■ const co and c area = c x — . 

cg.a \ } eg area x 2g - 2 + m 

For almost any flat surface S of unit area one has the following week 
convergence of measures: 

dirr, — > ip arL — > p pair^ — > ip x p, 

where ip is the uniform probability measure on the circle and p is the probability 
measure on the unit interval [0; 1] with the density (2g— 3+m)(l— x) 29 ~ +m dx. 



Flat Surfaces 75 



Directional flow: more serious reading. Theorem of S. Kerckhoff, H. Masur 
and J. Smillic is proved in [KMaS]. An example of minimal but not uniquely 
ergodic interval exchange transformations is constructed by W. Veech in [Vel] 
(using different terminology); an independent example (also using different 
terminology) was constructed at the same time by V. I. Oseledets. For flows 
such examples are constructed in the paper of A. Katok [Katl] and developed 
by E. Sataev in [Sat]. Another example was discovered by M. Keane [Kea2]. 
For alternative approach to the study of unique ergodicity of interval ex- 
change transformations see the paper of M. Boshernitzan [Ber2]. A very nice 
construction of minimal but not uniquely ergodic interval exchange transfor- 
mations (in a language which is very close to the language of this paper) can 
be found in the survey of H. Masur and S. Tabachnikov [MaT] or in the survey 
of H. Masur [Ma7]. 

6.2 Siegel Veech Formula 

We start from a slight formalization of our counting problem. As usual we 
start with a model case of the flat torus. As usual we assume that our flat 
torus is glued from a unit square. We count closed regular geodesies on T 2 of 
a bounded length. To mimic count of saddle connections we mark two points 
Pi P2 on T 2 and count geodesic segments of bounded length joining Pi and 
P2. 

Our formalization consists in the following construction. Consider an aux- 
iliary Euclidian plane M 2 . Having found a regular closed geodesic on T 2 we 
note its direction a and length / and draw a vector in R 2 in direction a hav- 
ing length I. We apply a similar construction to "saddle connections". The 
endpoints of corresponding vectors form two discrete subsets in IR 2 which we 
denote by V cg and V sc . 

It is easy to see that for the torus case a generic choice of Pi and P2 
generates a set V sc which is just a shifted square lattice, see Fig. 30. The 
set V cg is a subset of primitive elements of the square lattice, see Fig. 30. 
Since we count only regular closed geodesies which do not turn many times 
around themselves we cannot obtain elements of the form (kni,krii) with 
k, m, 712 £ Z. 

The corresponding counting functions N SC (T 2 ,L) and N SC (T 2 ,L) corre- 
spond to the number of element of V sc and V cg correspondingly which get to 
a disc of radius L centered in the origin. Both functions have exact quadratic 
asymptotics. Denoting by Xl{v), where dgI 2 the indicator function of such 
disc we get 

N SC (T 2 , L) = J2 veVsc Xl(v) ~ 1 • ttL 2 (22) 
N cg (T 2 ,L) = £ t , 6Vcg X l(v) -tL-kL 2 




The coefficients in quadratic asymptotics define the corresponding Siegel 

1 6 / 

Vccch constants const ,„ = 1 and const rn = —— = — . (Note that here we 

eg ^ ^2 v 

count every geodesic twice: once with one orientation and the other one with 
the opposite orientation. This explains why in this normalization we obtain 
the value of const cg twice as much as const cg for genus one in Table 2.) 

Consider now a more general flat surface S. Fix the geometric type of con- 
figuration C of saddle connections or closed geodesies. By definition all saddle 
connections (closed geodesies) in C are parallel and have equal length. Thus, 
similar to the torus case, every time we see a collection of saddle connections 
(closed geodesies) of geometric type C we can associate to such collection a 
vector in R 2 . We again obtain a discrete set V(S) C R 2 . 

Now fix C and apply this construction to every flat surface S in the stratum 
7i\(d\, . .. ,d m ). Consider the following operator / i— > / generalizing (22) from 
functions with compact support on R 2 to functions on H.i(di, . . . , d m ): 

f(S) = J2 /(«) 

Lemma (W. Veech). The functional 

/-> / f(S)dvi 

Jn c 1 omp (d 1 ,...,d m ) 

is SL(2,Wl)-invariant. 

Having proved convergence of the integral above the Lemma follows imme- 
diately from invariance of the measure dv\ under the action of SX(2, R) and 
from the fact that V(gS) = gV(s) for any flat surface S and any g G SL(2, R). 

Now note that there very few SL(2, R)-invariant functionals on functions 
with compact support in R 2 . Actually, there arc two such functionals, and the 
other ones are linear combinations of these two. These two functionals are the 
value of /(0) at the origin and the integral J R2 f{x,y)dxdy. It is possible to 
see that the value /(0) at the origin is irrelevant for the functional from the 
Lemma above. Hence it is proportional to J R2 f(x,y)dxdy. 
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Theorem (W. Veech). For any function f : M 2 — > R wrei/i compact support 
one has 



Here the constant C in (23) does not depend on /; it depends only on the 
connected component l-i c ° mv {d\^ . . . ,d m ) and on the geometric type C of the 
chosen configuration. 

Note that it is an exact equality, fn particular, choosing as / = \l the 
indicator function of a disc of radius L we see that for any given L G R, 
the average number of saddle connections not longer than L on flat surfaces 
S G Tii° mp (di, . . . , dm) is exactly C ■ irL 2 , where C does not depend on L. 

It would be convenient to introduce a special notation for such /. Let 



The Theorem of Eskin and Masur [EMa] cited above tells that for large values 
of L one gets /l(S') ~ c(C)ttL 2 for almost all individual flat surfaces S G 
Hl° mp (di, . . . ,d m ) and that the corresponding constant c(C) coincides with 
the constant C above. 

Formula (23) can be applied to /l for any value of L. In particular, in- 
stead of considering large L we can choose a very small value L = e. The 
corresponding function f e (S) counts how many collections of parallel e-short 
saddle connections (closed geodesies) of the type C we can find on the flat 
surface S. 

Usually there are no such saddle connections (closed geodesies), so for most 
flat surfaces f e (S) = 0. For some surfaces there is exactly one collection like 
this. We denote the corresponding subset by 7i\ ,thlck {C) C 7ii(di, . . . , d m ). 
Finally, for the surfaces from the remaining (very small) subset 7^'* m (C) 
one has several collections of short saddle connections (closed geodesies) of 
the type C. Thus, 




(23) 




for most of the surfaces S 



1 for S G Uf Uck (C) 
> 1 for S G Hf hin (C) 



and we can rewrite (23) for f e as 
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c(C)-ne =c(C)l Xe{x,y)dxdy 
1 



Vol Hi P (d U d m ) JHl om *(d u ...,d m ) 

1 



fe(S) dv x 
1 dv x + 



Vo\H c 1 omp (di, ...,d m ) J n ^,thia^ C) 

fe(S) dv x 



1 

VolHi° mp (di, . . . , d m ) Jn\' thin (c) 

It can be shown that though f £ (S) might be large on Hj 1 ' m (C) the mea- 
sure of this subset is so small (it is of the order e 4 that the last integral above 
is negligible in comparison with the previous one; namely it is o(e 2 ). (This is a 
highly nontrivial result of Eskin and Masur [EMa].) Taking into consideration 
that 

/ ldv! = VolH £ { thlck (C) = VolHf (C) + o(e 2 ) 

J H\' thzck {C) 

we can rewrite the chain of equalities above as 

2 VolWf(C) . ,. 

c(C)TC= VdCU...,4) + ^ 

which is equivalent to (21). 

Baby Case: Closed Geodesies on the Torus 

As an elementary application we can prove that proportion of primitive lattice 
points among all lattice points is 1/C(2). In other words, applying (21) we can 
prove asymptotic formula (22) for the number of primitive lattice points in a 
disc of large radius L. As we have seen at Fig. 30 this number equals to the 
number 7V cff (T 2 , L) of families of oriented closed geodesies of length bounded 
by L on the standard torus T 2 . 

We want to apply (21) to prove the following formula for the corresponding 
Siegel-Veech constant c+ (where superscript + indicates that we are counting 
oriented geodesies on T 2 ). 

A^ cg (T 2 ,L) 1 6 



l^L ttL 2 " C(2) 



r'2 



Note that the moduli space 7ii(0) of flat tori is a total space of a unit 
tangent bundle to the modular surface (see Sec. 3.2, Fig. 13; see also (30) in 
Sec. 9.1 for geometric details). Modular surface can be considered as a space 
of flat tori of unit area without choice of direction to the North. 

Measure on this circle bundle disintegrates to the product measure on 
the fiber and the hyperbolic measure on the modular curve. In particular, 
Vol(7ii(0)) = 7r-7r/3, where ir/3 is the hyperbolic area of the modular surface. 
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Similarly, Vol(Wf(0)) = ir ■ Arca(Cusp(e)), where Cusp(e) is a subset of 
the modular surface corresponding to those flat tori of unit area which have 
a geodesic shorter than e (see Fig. 13). 

Showing that Area(Cusp(e)) ~ e 2 we apply (21) to get 

1 Area(Cusp(e)) 1 e 2 + o(e 2 ) 1 

Cm = lim 



+ o 7T£ 2 Area(Modular surface) 7re 2 7r/3 2<"(2)" 

Note that the Siegel-Veech constant c cg corresponds to counting nonoriented 
closed geodesies on T 2 . Thus, finally we obtain the desired value c+ = 2c cg . 

In the next section we give an idea of how one can compute Vol TL\ (C) in the 
simplest case, In Sec. 6.4 we describe the phenomenon of higher multiplicities 
and discuss the structure of typical cusps of the moduli spaces 7i(di, . . . , d m ). 



6.3 Simplest Cusps of the Moduli Space 

In this section we consider the simplest "cusp" C on a stratum 7i(c?i, . . . , d m ) 
and evaluate VolWf (C) for this cusp. Namely, we assume that the flat surface 
has at least two distinct conical points Pi ^ P 2 ; let 27r(di + 1), 27r(di + 1) be 
corresponding cone angles. As a configuration C we consider a configuration 
when we have a single saddle connection p joining Pi to P2 and no other 
saddle connections on S parallel to p. In our calculation we assume that the 
conical points on every S € T~l{di, . . . , d m ) have names; we count only saddle 
connections joining Pi to Pi. 

Consider some S e (C) C H(di, d m ), that is a flat surface S 

having a single saddle connection joining Pi to P2 which is not longer than e 
and having no other short saddle connections or closed geodesies. 

We are going to show that there is a canonical way to shrink the saddle 
connection on S E TL £ ^ thlck [C) coalescing two conical points into one. We shall 
see that, morally, this provides us with an (almost) fiber bundle 

Hi'* %c (di, d2, ds, ■ ■ ■ , d m ) 

[dI (24) 

Hi(di +d 2 ,d 3 ,...,d m ) 

where D 2 is a ramified cover of order (di+efe+l) over a standard metric disc of 
radius e. Moreover, we shall see that the measure on lc (di , d% , d%, ... , d m ) 
disintegrates into a product of the standard measure on D 2 and the natural 
measure on 7ii(di + d%, d$, . . . , d m ). The latter would imply the following 
simple answer to our problem: 

Vol("e-neighborhood of the cusp C") = 

= Vo\H e 1 (C) = VolHf(d 1 ,d 2 ,d 3 ,...,d m ) w 

ss (di + d 2 + 1) ■ ne 2 ■ Vol + d 2 , d 3 ,..., d m ) (25) 
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Instead of contracting an isolated short saddle connection to a point we 
prefer to create it breaking a conical point P' £ S' of degree d\ + d 2 on a 
surface S' € H-i(di + d 2 , d^, . . . ,d m ) into two conical points of degrees di 
and d 2 joined by a short saddle connections. We shall sec that this surgery is 
invertiblc, and thus we shall get a coalescing construction. In the remaining 
part of this section we describe this surgery following [EMaZo] . 

Breaking up a Conical Point into Two 

We assume that the initial surface S' € 7ii(c?i + d 2 , d$, . . . , d m ) does not have 
any short saddle connections or short closed geodesies. 

Consider a metric disc of a small radius e centered at the point P' , i.e. 
the set of points Q' of the surface S' such that Euclidean distance from Q' 
to the point P 1 is less than or equal to e. We suppose that e > is chosen 
small enough, so that the e-disc does not contain any other conical points 
of the metric; we assume also, that the disc which we defined in the metric 
sense is homeomorphic to a topological disc. Then, metrically our disc has a 
structure of a regular cone with a cone angle 27r(<ii + d 2 + 1); here d\ + d 2 
is the multiplicity of the zero P' . Now cut the chosen disc (cone) out of the 
surface. We shall modify the flat metric inside it preserving the metric at the 
boundary, and then paste the modified disc (cone) back into the surface. 




Fig. 31. Breaking up a conical point into two (after [EMaZo]). 



Our cone can be glued from 2(fcj + 1) copies of standard metric half-discs 
of the radius e, see the picture at the top of Fig. 31. Choose some small 8, 
where < 6 < e and change the way of gluing the half-discs as indicated on 
the bottom picture of Fig. 31. As patterns we still use the standard metric 
half-discs, but we move slightly the marked points on their diameters. Now we 
use two special half-discs; they have two marked points on the diameter at the 
distance 8 from the center of the half disc. Each of the remaining 2ki half-discs 
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has a single marked point at the distance 5 from the center of the half-disc. 
We are alternating the half-discs with the marked point moved to the right 
and to the left from the center. The picture shows that all the lengths along 
identifications arc matching; gluing the half-discs in this latter way we obtain 
a topological disc with a flat metric; now the flat metric has two cone-type 
singularities with the cone angles 2ir(di + 1) and 2n(d2 + 1). 

Note that a small tubular neighborhood of the boundary of the initial cone 
is isometric to the corresponding tubular neighborhood of the boundary of the 
resulting object. Thus we can paste it back into the surface. Pasting it back 
we can turn it by any angle ip, where < ip < 2ir(di + c?2 + 1). 

We described how to break up a zero of multiplicity d\ + d2 of an Abelian 
differential into two zeroes of multiplicities d±,d2- The construction is local; 
it is parameterized by the two free real parameters (actually, by one complex 
parameter): by the small distance 26 between the newborn zeroes, and by 
the direction ip of the short geodesic segment joining the two newborn zeroes. 
In particular, as a parameter space for this construction one can choose a 
ramified covering of degree d\ + e?2 + 1 over a standard metric disc of radius 
e. 

6.4 Multiple Isometric Geodesies and Principal Boundary of the 
Moduli Space 

In this section we give an explanation of the phenomenon of higher multiplic- 
ities, we consider typical degenerations of flat surfaces and we discuss how 
can one use configurations of parallel saddle connections or closed geodesies 
to determine the orbit of a flat surface S. 

Multiple Isometric Saddle Connections 

Consider a collection of saddle connections and closed geodesies representing 
a basis of relative homology H 1 (S, {Pi, . . . , P m }; Z) on a flat surface S. Re- 
call, that any geodesic on S goes in some constant direction. Recall also that 
by Convention 1 any flat surface is endowed with a distinguished direction 
to the North, so we can place a compass at any point of S and determine 
in which direction goes our geodesic. Thus, every closed geodesic or saddle 
connection determines a vector v^- € M 2 ~ C which goes in the same direc- 
tion and have the same length as the corresponding geodesic element. Finally 
recall that collection of planar vectors {vi, . . . , V2 fl + m -i} considered as com- 
plex numbers provide us with a local coordinate system in 7i(di, . . . , d m ). In 
complex-analytic language these coordinates are the relative periods of holo- 
morphic 1-form representing the flat surface S, namely Vj = J p co, where pj 
is the corresponding geodesic element (saddle connection or closed geodesic). 

We say that two geodesic elements 71,72 (saddle connections or closed 
geodesies) are homologous on a flat surface S if they determine the same 
homology classes in H\(S, {Pi, ... , P m }; Z). In other words, 71 is homologous 



82 



Anton Zorich 



to 72 if cutting S by 71 and 72 we break the surface S into two pieces. For 
example, the saddle connections 71,72,73 on the right surface at the bottom 
of Fig. 32 are homologous. 

The following elementary observation is very important for the sequel. 
Since the holomorphic 1-form u> representing S is closed, homologous geodesic 
elements 71,72 define the same period: 



We intensively used the following result of H. Masur and J. Smillie [MaS] 
in our considerations in the previous section. 

Theorem (H. Masur, J. Smillie). There is a constant M such that for 
all e, k > the subset H\{di, . . . , d m ) of Hi(di, . . . , d rn ) consisting of those 
flat surfaces, which have a saddle connection of length at most e, has volume 
at most Me 2 . 

The volume of the set of flat surfaces with a saddle connection of length 
at most e and a nonhomologous saddle connection with length at most k is at 
most Me 2 n 2 . 

Morally, this Theorem says that a subset corresponding to one complex 
parameter with norm at most e has measure of order e 2 , and a subset cor- 
responding to two complex parameter with norm at most e has measure of 
order e 4 . In particular, this theorem implies that Yo\'H E 1 ' thm {d\, . . . , d m ) is of 
the order e 4 . 

In the previous section we considered the subset of flat surfaces S G 
Jt\ (e?i, c?2, <^3, ■ ■ • , d m ) having a single short saddle connection joining ze- 
roes of degrees d\ and d^- We associated to such surface S a new surface 
S' G Hf(di + d%, g?3, . . . , d m ) in the smaller stratum. Note that, morally, sur- 
faces S and S' have the same periods with a reservation that S has one more 
period than S': the extra small period represented by our short saddle con- 
nection. 

Metrically surfaces 5 and S' are almost the same: having a surface S we 
know how to contract our short saddle connection to a point; having a surface 
S' and an abstract short period v G C ~ M 2 we know how to break the corre- 
sponding zero on S' into two zeroes joined by a single short saddle connection 
realizing period v. (In the latter construction we have some additional discrete 
freedom: we can break the zeroes in direction v in d\ + efo + 1 different ways.) 

Our construction does not work when we have two nonhomologous short 
geodesic elements on the surface S. But we do not care since according to the 
Theorem of H, Masur and J. Smillie the subset Wf'* ,n (di, . . . , d m ) of such 
surfaces has very small measure (of the order e 4 ). 

Now consider a slightly more general surgery represented by Fig. 32. We 
take three distinct flat surfaces, we break a zero on each of them as it was 
done in the previous section. We do it coherently using the same direction 
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and the same distance 5 on each surface (left part of Fig. 32). Then we slit 
each surface along the newborn saddle connection and glue the surfaces in a 
close compound surface as indicated on the right part of Fig. 32. 




Fig. 32. Multiple homologous saddle connections; topological picture 



The resulting surface has three short parallel saddle connections of equal 
length. By construction they are homologous: cutting along any two of them 
we divide the surface into two parts. Thus, the resulting surface again has only 
one short period! Note that a complete collection of periods of the compound 
surface can be obtained as disjoint union 

periods of S[ U periods of S' 2 U periods of £3 U newborn short period 

Hence, any flat surface S with three short homologous saddle connections and 
no other short geodesic elements can be viewed as a nonconnected flat surface 
S[ U £2 LI £3 plus a memory about the short period vgC which we use when 
we break the zeroes (plus some combinatorial arbitrariness of finite order). 

The moduli space of disconnected flat surfaces of the same type as 
S[ U S' 2 LI S3 has one dimension less than the original stratum H(d\, . . . , d m ). 
Our considerations imply the following generalization of formula (25) for the 
volume of "e-neighborhood" of the corresponding cusp: 



Vol("e-neighborhood of the cusp C" ) = VolHf (C) = 
k- Vol Hi (stratum of S[)-Vo\Hi (stratum of S' 2 )- Vol Hi (stratum of S'^-Tie 2 
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where the factor k is an explicit constant depending on dimensions, possible 
symmetries, and combinatorics of multiplicities of zeroes. In particular, we 
get a subset of order e 2 . 

Now for any flat surface So € Ti(di, . . . , d m ) we can state a counting 
problem, where we shall count only those saddle connections which appear 
in configuration C of triples of homologous saddle connections breaking So 
into three surfaces of the same topological and geometric types as S[ , S' 2 , S' 3 . 
Applying literally same arguments as in Sec. 6.2 and 6.3 we can show that 
such number of triples of homologous saddle connections of length at most L 
has quadratic growth rate and that the corresponding Siegel-Veech constant 
c(C) can be expressed by the same formula as above: 

1 Vol("e-neighborhood of the cusp C ") 
C ( C ) = ]™ VolHT mp (d u ...,d m ) = 

IlLi VolHi(stratum of S' k ) 



(combinatorial factor) 



vo\nT mp (d 1 ,...,d m ) 



Principal Boundary of the Strata 

The results above give a description of typical degenerations of flat surfaces. 
A flat surface gets close to the boundary of the moduli space when some 
geodesic element (or a collection of geodesic elements) get short. Morally, we 
have described something like "faces" of the boundary, while there are still 
"edges" , etc, representing degenerations of higher codimension. Flat surfaces 
which are close to this "principal boundary" of a stratum 7i(c?i, . . . , d m ) have 
the following structure. 

If the short geodesic element is a saddle connection joining two distinct 
zeroes, then the surface looks like the one at Fig. 32. It can be decomposed to 
several flat surfaces with slits glued cyclically one to another. The boundaries 
of the slits form short saddle connections on the compound surface. All these 
saddle connections join the same pair of points; they have the same length and 
direction. They represent homologous cycles in the relative homology group 
H 1 {S,{P 1 ,...,P m };1). 

Figure 33 represents a flat surface So € 7i(l, 1) unfolded to a polygon. The 
two short bold segments represent two homologous saddle connections. The 
reader can easily check that on any surface S € Tt(l, 1) obtained from So by 
a small deformation one can find a pair of short parallel saddle connections 
of equal length. Cutting S by these saddle connections we get a pair of tori 
with slits. 

We did not discuss in the previous section the case when the short geodesic 
element is a regular closed geodesic (or a saddle connection joining a conical 
point to itself). Morally, it is similar to the case of saddle connections, but 
technically it slightly more difficult. A flat surface near the principal boundary 
of this type is presented on the right of Fig. 34. 
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Fig. 33. Flat surface with two short homologous saddle connections. Any small 
deformation of this surface also has a pair of short homologous saddle connections. 




Fig. 34. A typical (on the right) and nontypical (on the left) degenerations of a 
flat surface. Topological picture 

Similar to the case of saddle connections, the surface can be also decom- 
posed to a collection of well-proportional flat surfaces S[ , . . . , S' n of lower 
genera. Each surface S' k has a pair of holes. Each of these holes is realized by 
a saddle connection joining a zero to itself. The surfaces are cyclically glued 
to a "necklace" , where two neighboring surfaces might be glued directly or by 
a narrow cylinder. Since the waist curves of all these cylinders and all sad- 
dle connections representing boundaries of surfaces S' k are homologous, the 
corresponding closed geodesies on S are parallel and have equal length. A 
more artistic 5 image of a surface, which is located closed to the boundary of 
a stratum is represented on Fig. 35. 

The surface on the left of Fig. 34 is close to an "edge" of the moduli space 
in the sense that it represents a "nontypical" degeneration: a degeneration 
of codimension two. This surface has two nonhomologous closed geodesies 
shorter than e. Due to the Theorem of H. Masur and J. Smillic cited above, 
the subset Ti^f m (di, . . . , d m ) of such surfaces has measure of the order e 4 . 
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Fig. 35. Flat surface near the principal boundary 5 



Configurations of Saddle Connections and of Closed Geodesies as Invariants 
of Orbits 

Consider a flat surface So and consider its orbit under the action of SL(2, K). 
It is very easy to construct this orbit locally for those elements of the group 
GI/ + (2,R) which are close to identity. It is a fairly complicated problem to 
construct this orbit globally in 7i(di, . . . , d m ) and to find its closure. Ergodic 
theorem of H. Masur and W. Veech (see Sec. 3.5) tells that for almost ev- 
ery surface So the closure of the orbit of So coincides with the embodying 
connected component of the corresponding stratum. But for some surfaces 
the closures of the orbits are smaller. Sometimes it is possible to distinguish 
such surfaces looking at the configurations of parallel closed geodesies and 
saddle connections. Say, for Veech surfaces which will be discussed in Sec. 9.5 
the orbit of SL(2,M.) is already closed in the stratum, so Veech surfaces are 
very special. This property has an immediate reflection in behavior of parallel 
closed geodesies and saddle connections: as soon as we have a saddle connec- 
tion or a closed geodesic in some direction on a Veech surface, all geodesies 
in this direction are either closed or (finite number of them) produce saddle 
connections. 

Thus, it is useful to study configurations of parallel closed geodesies on 
a surface (which includes the study of proportions of corresponding maximal 
cylinders filled with parallel regular closed geodesies) to get information about 
the closure of corresponding orbit. 

One can also use configurations of parallel closed geodesies on a flat surface 
to determine those connected component of the stratum, to which belongs 
our surface So- Some configurations (say, g — 1 tori connected in a "necklace" 
by a chain of cylinders, compare to Fig. 35) are specific for some connected 

5 H. Matisse: La Danse. The State Hermitage Museum, St. Petersburg, (c) Suc- 
cession H. Matisse/VG Bild-Kunst, Bonn, 2005 
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components and never appear for surfaces from other connected components. 
We return to this discussion in the very end of Sec. 9.4 where we use this 
idea to distinguish connected components of the strata in the moduli space of 
quadratic differentials. 

6.5 Application: Billiards in Rectangular Polygons 

Consider now a problem of counting generalized diagonals of bounded length 
or a problem of counting closed billiard trajectories of bounded length in a 
billiard in a rational polygon 77. Apply Katok-Zcmliakov construction (see 
Sec. 2.1) and glue a very flat surface S from the billiard table 77. Every 
generalized diagonal (trajectory joining two corners of the billiard, possibly, 
after reflections from the sides) unfolds to a saddle connection and every 
periodic trajectory unfolds to a closed regular geodesic. 

It is very tempting to use the results described above for the counting 
problem for the billiard. Unfortunately, the technique elaborated above is not 
applicable to billiards directly The problem is that flat surfaces coming from 
billiards form a subspace of large codimension in any stratum of flat surfaces; 
in particular, this subspace has measure zero. Our "almost all" technique does 
not see this subspace. 

However, the problems are related in some special cases; see [EMaScm] 
treating billiard in a rectangle with a barrier. As another illustration we 
consider billiards in "rectangular polygons". These results represent the 
work [AthEZo] which is in progress. We warn the reader that we are extremely 
informal in the remaining part of this section. 

Rectangular Polygons 

Figure 36 suggests several examples of rectangular polygons. The "polygons" 
are allowed to have ramification points at the boundary, with restriction that 
the angles at ramification points are integer multiples of tt/2. Note that we 
do not identify the side P$Pq with a part of the side P4P5 in the right poly- 
gon. This polygon should be considered with a cut along the side PsTV The 
corresponding billiard has a "barrier" along the side Pr^Pe- 




PC > 2 

Fig. 36. Rectangular polygons 



88 Anton Zorich 



Formally speaking, by a rectangular polygon II we call a topological disc 
endowed with a flat metric, such that the boundary dll is presented by a 
finite broken line of geodesic segments and such that the angle between any 
two consecutive sides equals kir/2, where k G N. 

Consider now our problem for a standard rectangular billiard table (the 
proportions of the sides do not matter). We emit a trajectory from some corner 
of the table and want it arrive to another corner after several reflections from 
the sides. 

When our trajectory reflects from a side it is convenient to prolong it as 
a straight line by making a reflection of the rectangle with respect to the 
corresponding side (see Katok-Zemliakov construction in Sec. 2.1). Unfolding 
our rectangular table we tile the plane M. 2 with a rectangular lattice. Our 
problem can be reformulated as a problem of counting primitive lattice points 
(see the right part of Fig. 30). 

We are emitting our initial trajectory from some fixed corner of the billiard. 
It means that in the model with a lattice in the plane we are emitting a straight 
line from the origin inside one of the four quadrants. Thus, we are counting 
the asymptotics for the number of primitive lattice points in the intersection 
of a coordinate quadrant with a disc of large radius L centered at the origin. 
This gives us 1/4 of the number of all primitive lattice points. Note that 
in our count we have fixed the initial corner, but we let our trajectory hit 
any of the remaining three corners. Thus, if we count only those generalized 
diagonals which are launched from some prescribed corner Pj and arrive to 
a prescribed corner Pj (different from initial one) we get 1/3 of the previous 
number. Hence, the number Nij(L) of generalized diagonals joining p with Pj 
is 1/12 of the number N cg (T 2 ,L) ~ (l/£(2)) • irL 2 of primitive lattice points, 
see 22. 

In our calculation we assumed that the billiard table has area one. It is 
clear that asymptotics for our counting function is homogeneous with respect 
to the area of the table. Adjusting our formula for a rectangular billiard table 
of the area different from 1 we get the following answer for the number of 
generalized diagonals of length at most L joining prescribed corner Pj to a 
prescribed corner Pj different from the first one: 

1 L 2 
Nij (L) ~ — ■ Area of the billiard table ^ 

Now, having studied a model case, we announce two examples of results 
from [AthEZo] concerning rectangular polygons. 

Consider a family of rectangular polygons having exactly k > angles 
37r/2 and all other angles ir/2 (see Fig. 37). Consider a generic billiard table 
in this family (in the measure-theoretical sense). Fix any two corners Pj =/= Pj 
having angles 7r/2. The number iVy(L) of generalized diagonals of length at 
most L joining Pj to Pj is approximately the same as for a rectangle: 

r ,n 1 L 2 
Nij(L)~—- 



2tt Area of the billiard table 
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Pi I T P . 

Fig. 37. Family of rectangular polygons of the same geometry and the same area. 
The shape of these polygons is quite different. Nevertheless for both billiard tables 
the number of trajectories of length at most L joining the right-angle corner Pi to 
the right-angle corner Pj is approximately the same as the number of trajectories 
of length at most L joining two right-angle corners of a rectangle of the same area. 



We have to admit that we are slightly cheating here: the equivalence "~" 
which we can currently prove is weaker than in (26); nevertheless, we do 
not want to go into technical details. 

Note that the shape of the polygon within the family may vary quite con- 
siderably, see Fig. 37, and this does not affect the asymptotic formula. How- 
ever, the answer changes drastically when we change the family. For rectan- 
gular polygons having several angles of the form mr the constant in quadratic 
asymptotics is more complicated. This is why we do not expect any elemen- 
tary proof of this formula (our proof involves evaluation of corresponding 
Siegel-Veech constant). 

Actually, naive intuition does not help in counting problems of this type. 
Consider, for example, an L-shaped billiard table as on Fig. 38. 




©Moon Duchin 

Fig. 38. L-shaped billiard table 
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The angle at the vertex P is Stt/2 which is three times larger than the 
angle n/2 at the other five vertices Pi, . . . ,P 5 . However, the number 

^ 2 I? 

3 it Area of the billiard table 

of generalized diagonals of length at most L joining P$ to Pj, where 1 < j < 5, 
is four times bigger than the number Nij(L) of generalized diagonals joining 
two corners with the angles tt/2. 



7 Volume of Moduli Space 

In Sec. 3.4 we defined a volume clement in the stratum H(di, . . . , d m ). We used 
linear volume element in cohomological coordinates i? 1 (S', {Pi, . . . , P m }', C) 
normalized in such way that a fundamental domain of the lattice 

H\S,{P 1 ,...,P m };Z®V^lZ) c H 1 (S,{P 1 ,...,P m };C) 

has unit volume. The unit lattice does not depend on the choice of cohomo- 
logical coordinates, its vertices play the role of integer points in the moduli 
space Tt(di, . . . , d m ). In Sec. 7.1 we suggest a geometric interpretation of flat 
surfaces representing integer points of the strata. 

Using this interpretation we give an idea for counting the volume ( "hyper- 
area") of the hypersurfacc 7i-i{d\, . . . , d m ) C 7i{d\, . . . , d m ) of flat surfaces of 
area one. We apply the strategy which can be illustrated in a model example 
of evaluation of the area of a unit sphere. We first count the asymptotics for 
the number N(R) of integer points inside a ball of huge radius R. Clearly 
N(R) corresponds to the volume of the ball, so if we know the asymptotics 
for N(R) we know the formula for the volume Vol(i?) of the ball of radius R. 
d , . 

The derivative -^yol(R)\ R=1 gives us the area of the unit sphere. 

Similarly, to evaluate the "hyperarea of a unit hyperboloid" Hi (d± , . . . , d m ) 
it is sufficient to count the asymptotics for the number of integer points inside 
a "hyperboloid" Hfi(di, . . . , d m ) of huge "radius" R. The role of the "radius" 
R is played by the positive homogeneous real function R = area(S) defined 
on H(di, . . . , d m ). 

Note that the volume v(l~l<R(di, . . . , d m )) of a domain bounded by the 
"hyperboloid" 7iR{d\, . . . , d m ) is a homogeneous function of R of the weight 
diniR TC(di, . . . , d m )/2 while the volume of a ball of radius R is a homogeneous 
function of R of the weight which equals the dimension of the space. This 
difference in weights explains the factor 2 in the formula below: 

Vo\(H 1 (d 1 ,...,d m )) = 2-^y{H< R {d 1 ,...,d m )) 

= dim K (Hi(di, . . .,d m )) ■ v(H<i{d\, . . .,d m )) (27) 
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This approach to computation of the volumes was suggested by A. Es- 
kin and A. Okounkov and by M. Kontsevich and the author. However, the 
straightforward application of this approach, described in Sec. 7.1, gives the 
values of the volumes only for several low-dimensional strata. The general so- 
lution of the problem was found by A. Eskin and A. Okounkov who used in 
addition powerful methods of representation theory. We give an idea of their 
method in Sec. 7.2. 

7.1 Square-tiled Surfaces 

Let us study the geometric properties of the flat surfaces S represented by 
"integer points" S £ H 1 ^, {Pi, . . . , P m }', Z © y^lZ) in cohomological co- 
ordinates. Let oj be the holomorphic one-form representing such flat surface 
S. Since [w] £ H 1 ^, {Pi, . . . , P m }; Z © \/ r lZ) all periods of uj (including 
relative periods) belong to Z © ^/— 1 Z. Hence the following map / w from the 
fiat surface S to the standard torus T 2 = C/(Z © \J — 1 Z) is well-defined: 



where Pi is one of the conical points. It is easy to check that f u is a ramified 
covering, moreover, it has exactly m ramification points, and the ramification 
points are exactly the zeros Pi,.. . , P m of u>. Consider the flat torus T 2 as 
a unit square with the identified opposite sides. The covering f u : S — ► T 2 
induces a tiling of the flat surface S by unit squares. Note that all unit squares 
are endowed with the following additional structure: we know exactly which 
edge is top, bottom, right, and left; adjacency of the squares respects this 
structure in a natural way: we glue vertices to vertices and edges to edges, 
moreover, the right edge of a square is always identified to the left edge of some 
square and top edge is always identified to the bottom edge of some square. 
We shall call a flat surface with such tiling a square-tiled surface. We see that 
the problem of counting the volume of 7ii(di, . . . ,d m ) is equivalent to the 
following problem: how many square-tiled surfaces of a given geometric type 
(determined by number and types of conical singularities) can we construct 
using at most N unit squares. Say, Fig. 46 gives the list of all square-tiled 
surfaces of genus g > 1 glued from at most 3 squares; they all belong to TC(2). 

In terms of the coverings our Problem can be formulated as follows. Con- 
sider the ramified coverings p : S — > T 2 over the standard torus T 2 . Fix the 
number m of branching points, and ramification degrees d%, . . . , d m ) at these 
points. Assume that all ramification points Pi,... , P m project to the same 
point of the torus T 2 . Enumerate ramified coverings of any given ramification 
type having at most N 3> 1 sheets. Here pairs of coverings forming commu- 
tative diagrams as below are identified: 




S 



S 



\ y 

T 2 



(28) 
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Counting of Square-tiled Tori 

Let us count the number of square-tiled tori tiled by at most N 3> 1 squares. 
In this case our square-tiled surface has no singularities at all: we have a flat 
torus tiled with the unit squares in a regular way. Cutting our flat torus by a 
horizontal waist curve we get a cylinder with a waist curve of length w S N 
and a height del, see Fig. 39. The number of squares in the tiling equals 
w ■ h. We can reglue the torus from the cylinder with some integer twist t, 
see Fig. 39. Making an appropriate Dehn twist along the waist curve we can 
reduce the value of the twist t to one of the values 0, 1, . . . , w — 1. In other 
words, fixing the integer perimeter w and height h of a cylinder we get w 
different square-tiled tori. 




Fig. 39. A square-tiled surface is decomposed into several cylinders. Each cylinder is 
parametrized by its width (perimeter) w and height h. Gluing the cylinders together 
we get also a twist parameter t, where < t < w, for each cylinder 



Thus the number of square tiled tori constructed by using at most N 
squares is represented as 

v-1 ( NY N 2 N 2 tt 2 

wh<N W< N 

— h 

Actually, some of the tori presented by the first sum are equivalent by an 
affinc diffcomorphism, so we are counting them twice, or even several times. 
Say, the tori w = 2;h = l;t = and w = 1; h = 2; t = are equivalent. 
However, this happens relatively rarely, and this correction does not affect 
the leading term, so we simply neglect it. 

d 

Applying the derivative 2^^- (see (27)) we finally get the following 

AT=1 

value for the volume of the space of flat tori 



Vol(Wi(0)) 
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Decomposition of a Square- Tiled Surface into Cylinders 

Let us study the geometry of square-tiles surfaces. Note that all leaves of 
both horizontal and vertical foliation on every square-tiled surface are closed. 
In particular the union of all horizontal critical leaves (the ones adjacent to 
conical points) forms a finite graph r. The collection Pi, ... , P m of conical 
points forms the set of the vertices of this graph; the edges of the graph are 
formed by horizontal saddle connections. The complement S — r is a union 
of flat cylinders. 

For example, for the square-tiled surfaces from Fig. 46 we get the following 
decompositions into horizontal cylinders. We have one surface composed from 
a single cylinder filled with closed horizontal trajectories; this cylinder has 
width (perimeter) w = 3 and hight h = 1. Two other surfaces are composed 
from two cylinders. The heights of the cylinders are hi = h% = 1, the widths 
are wi = 1 and u>2 = 2 correspondingly. Observing the two-cylinder surfaces 
at Fig. 46 we see that they differ by the twist parameter £2 (see Fig. 39) of the 
wider cylinder: in one case £2 = and in the other case £2 = 1. By construction 
the width iUj and height hi of any cylinder are strictly positive integers; the 
value of the twist U is a nonncgative integer bounded by the width of the 
cylinder: < U < Wi- 

Separatrix Diagrams 

Let us study in more details the graphs r of horizontal saddle connections. 

We start with an informal explanation. Consider the union of all saddle 
connections for the horizontal foliation, and add all critical points (zeroes 
of oj). We obtain a finite oriented graph r . Orientation on the edges comes 
from the canonical orientation of the horizontal foliation. Moreover, graph _T 
is drawn on an oriented surface, therefore it carries so called ribbon structure 
(even if we forget about the orientation of edges), i.e. on the star of each vertex 
P a cyclic order is given, namely the counterclockwise order in which edges 
are attached to P. The direction of edges attached to P alternates (between 
directions toward P and from P) as we follow the counterclockwise order. 

It is well known that any finite ribbon graph r defines canonically (up to 
an isotopy) an oriented surface S(r) with boundary. To obtain this surface 
we replace each edge of J 1 by a thin oriented strip (rectangle) and glue these 
strips together using the cyclic order in each vertex of r. In our case surface 
S(r) can be realized as a tubular e-ncighborhood (in the sense of transversal 
measure) of the union of all saddle connections for sufficiently small e > 0. 

The orientation of edges of r gives rise to the orientation of the bound- 
ary of S(r). Notice that this orientation is not the same as the canonical 
orientation of the boundary of an oriented surface. Thus, connected compo- 
nents of the boundary of S(r) are decomposed into two classes: positively 
and negatively oriented (positively when two orientations of the boundary 
components coincide and negatively, when they are different). The comple- 
ment to the tubular e-neighborhood of r is a finite disjoint union of open 
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cylinders foliated by oriented circles. It gives a decomposition of the set of 
boundary circles ir (d(S(r))) into pairs of components having opposite signs 
of the orientation. 

Now we are ready to give a formal definition: 

A separatrix diagram is a finite oriented ribbon graph r , and a decompo- 
sition of the set of boundary components of S(r) into pairs, such that 

the orientation of edges at any vertex is alternated with respect to the 
cyclic order of edges at this vertex; 

there is one positively oriented and one negatively oriented boundary com- 
ponent in each pair. 

Notice that ribbon graphs which appear as a part of the structure of a 
separatrix diagram are very special. Any vertex of such a graph has even 
degree, and the number of boundary components of the associated surface 
with boundary is even. Notice also, that in general the graph of a separatrix 
diagram is not planar. 

Any separatrix diagram (repairing) defines a closed oriented surface to- 
gether with an embedding of r (up to a homcomorphism) into this surface. 
Namely, we glue to the surface with boundary S(r) standard oriented cylin- 
ders using the given pairing. 

In pictures representing diagrams we encode the pairing on the set of 
boundary components painting corresponding domains in the picture by some 
colors (textures in the black-and-white text) in such a way that every color 
appears exactly twice. We will say also that paired components have the same 
color. 




Fig. 40. An example of a separatrix diagram. A detailed picture on the left can be 
encoded by a schematic picture on the right. 
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Example. The ribbon graph presented at Figure 40 corresponds to the hor- 
izontal foliation of an Abclian differential on a surface of genus g = 2. The 
Abclian differential has a single zero of degree 2. The ribbon graph has two 
pairs of boundary components. 

Any separatrix diagram represents an orientable measured foliation with 
only closed leaves on a compact oriented surface without boundary. We say 
that a diagram is realizable if, moreover, this measured foliation can be chosen 
as the horizontal foliation of some Abelian differential. Lemma below gives a 
criterion of realizability of a diagram. 

Assign to each saddle connection a real variable standing for its "length" . 
Now any boundary component is also endowed with a "length" obtained as 
sum of the "lengths" of all those saddle connections which belong to this 
component, ff we want to glue flat cylinders to the boundary components, the 
lengths of the components in every pair should match each other. Thus for 
every two boundary components paired together (i.e. having the same color) 
we get a linear equation: "the length of the positively oriented component 
equals the length of the negatively oriented one" . 

Lemma. A diagram is realizable if and only if the corresponding system of 
linear equations on "lengths" of saddle connections admits strictly positive 
solution. 

The proof is obvious. 

Example. The diagram presented at Fig. 40 has three saddle connections, all 
of them are loops. Let Pi6,P52,P34 be their "lengths". There are two pairs 
of boundary components. The corresponding system of linear equations is as 
follows: 

f P34 = Pie 

\ Pl6 + P52 = P34 + P52 

Exercise. Check that two separatrix diagrams at Fig. 4f are realizable, and 
one - not. Check that there are no other realizable separatrix diagrams for 
the surfaces from the stratum 7i(2). Find all realizable separatrix diagrams 
for the stratum H(l, 1). 

Counting of Square-tiled Surfaces in Ti(2) 

To consider one more example we count square-tiled surfaces in the stratum 
H(2). We have seen that in this stratum there are only two realizable separa- 
trix diagrams; they are presented on the left and in the center of Fig. 41. 

Consider those square tiled surfaces from H(2) which correspond to the 
left diagram from Fig. 41. In this case the ribbon graph corresponding to the 
separatrix diagram has single "top" and single "bottom" boundary compo- 
nent, so our surface is glued from a single cylinder. The waist curve of the 
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Fig. 41. The separatrix diagrams from left to right a) represent a square-tiled 
surface glued from one cylinder of width w = p\ + P2 + P3', b) represent a square- 
tiled surface glued from two cylinder of widths wi = pi and W2 — pi +P2; c) do not 
represent any square-tiled surface 

cylinder is of the length w = p\ + P2 + P3, where p\,p2,pz are the lengths of 
the separatrix loops. As usual, denote the height of the cylinder by h. The 
twist t of the cylinder has an integer value in the interval [0, w — 1}. Thus the 
number of surfaces of this type of the area bounded by N is asymptotically 
equivalent to the sum 

3 ^ (Pi+f2+P3)~^-C(4) = ^-^ 

(pi+P2+P3)h<N 

(see more detailed computation in [Zo5]). The coefficient 1/3 compensates 
the arbitrariness of the choice of enumeration of pi,p2,p% preserving the cyclic 
ordering. Similar to the torus case we have neglected a small correction coming 
from counting equivalent surfaces several times. 

Exercise. Check that for pi = P2 = P3 = 1 all possible values of the twist 
t = 0, 1, 2 give equivalent flat surfaces; see also Fig. 46 

Consider now a ribbon graph corresponding to the middle diagram from 
Fig. 41. It has two "top" and two "bottom" boundary components. Thus, 
topologically, we can glue in a pair of cylinders in two different ways. However, 
to have a flat structure on the resulting surface we need to have equal lengths 
of "top" and "bottom" boundary components. These lengths are determined 
by the lengths of the corresponding separatrix loops. It is easy to check that 
one of the two possible gluings of cylinders is forbidden: it implies that one of 
the separatrix loops has zero length, and hence the surface is degenerate. 

The other gluing is realizable. In this case there is a pair of separatrix 
loops of equal lengths p±, see Fig. 41. The square-tiled surface is glued from 
two cylinders: one having a waist curve wi = Pi, and the other one having 
waist curve W2 = Pi +P2- Denote the heights and twists of the corresponding 
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cylinders by hi , h 2 and ti , i 2 . The twist of the first cylinder has the value in 
the interval [0, W\ — 1]; the twist of the second cylinder has the value in the 
interval [0, w 2 — 1]- Thus the number of surfaces of 2-cylinder type of the area 
bounded by N is asymptotically equivalent to the sum 



E 

Pl,P2,hi,h 2 
Plh 1 + (p 1 +p 2 )h 2 <N 



Plipi +P2) 



E 



p\ 



P1P2 



Pi,p 2 ,hi,h 2 
p 1 (h 1 +h 2 )+p 2 h 2 <N 



It is not difficult to represent these two sums in terms of the multiple zeta 
values C(l>3) and C(2, 2) (see detailed computation in [Zo5]). Applying the 
relations <;(1, 3) = C(4)/4 and C(2,2) = 3C(4)/4 we get the following asymp- 
totic formula for our sum: 



E 



P1+P1P2 



N 4 N 4 5 

2<(1,3)+C(2,2)) = — - T C(4) 



Pi,p 2 ,hi,h 2 
pi(h 1 +h 2 )+p 2 h 2 <N 



24 



24 4 



Joining the impacts of the two diagrams and applying 2 ■ 



dN 



N 4 5 7T 4 

24 '4' 90 



(see (27)) 



we get the following value for the volume of the stratum H(2): 

Vol(Hi(2)) 



N=l 



7T 

120 



Separatrix diagrams: more serious reading. Technique of separatrix diagrams 
is extensively used by K. Strebel in [Str] and in some articles like [KonZo]. 



7.2 Approach of A. Eskin and A. Okounkov 

It is time to confess that evaluation of the volumes of the strata by means 
of naive counting square-tiled surfaces suggested in the previous section is 
not efficient in general case. The number of realizable separatrix diagrams 
grows and it is difficult to express the asymptotics of the sums for individual 
separatrix diagrams in reasonable terms (say, in terms of multiple zeta values). 
In general case the problem was solved using the following approach suggested 
by A. Eskin and A. Okounkov in [EOk]. 

Consider a square-tiled surface S G 7i(e?i, . . . , d m ). Enumerate the squares 
in some way. For the square number j let TT r (j) be the number of its neighbor 
to the right and let it u (j) be the number of the square atop the square number 
j. Consider the commutator tt 1 = 7r r 7r tl 7r J T 1 7r~ 1 of the resulting permutations. 
When the total number of squares is big enough, for most of the squares 

Geometrically the resulting permutation tt' corresponds to the following 
path: we start from a square number j, then we move one step right, one 
step up, one step left, one step down, and we arrive to n'(j). When the total 
number of squares is large, then for majority of the squares such path brings 
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us back to the initial square; for these values of index j we get 7r'(j) = j. 
However, we may have more than 4 squares adjacent to a conical singularity 
Pk £ S. For the squares adjacent to a singularity the path right-up-left-down 
docs not bring us back to the initial square. It is easy to check that the 
commutator ir' = 7r r 7r ll 7r r T 1 7r~ 1 is represented as a product of m cycles of 
lengths (di + 1), . . . , (d m + 1) correspondingly 

We conclude that a square-tiled surface S € H(di, . . . ,d m ) can be defined 
by a pair of permutations 7r r ,7r u , such that the commutator 7r r 7r u 7r ) T 1 7r, [ 7 1 
decomposes into given number m of cycles of given lengths (di + 1), . . . , (d m + 
1). Choosing another enumeration of the squares of the same square-tiled 
surface S we get two new permutations ir r ,ir u . Clearly the permutations in this 
new pair are conjugate to the initial ones by means of the same permutation 
p responsible for the change of enumeration of the squares: w r = pn r p~ x , ir u = 

piTup- 1 . 

We see that the problem of enumeration of square-tiled surfaces can be 
reformulated as a problem of enumeration of pairs of permutations of at most 
N elements such that their commutator decomposes into a given number 
of cycles of given lengths. Here the pairs of permutations are considered up 
to a simultaneous conjugation. This problem was solved by S. Bloch and 
A. Okounkov by using methods of representation theory. However, it is not 
directly applicable to our problem. Describing the square-tiled surfaces in 
terms of pairs of permutations one has to add an additional explicit constraint 
that the resulting square-tiled surface is connectedl Taking a random pair of 
permutations of very large number N 3> 1 of elements realizing some fixed 
combinatorics of the commutator we usually get a disconnected surface! 

The necessary correction is quite nontrivial. It was performed by A. Eskin 
and A. Okounkov in [EOk]. In the further paper A. Eskin, A. Okounkov 
and R. Pandharipande [EOkPnd] give the volumes of all individual connected 
components of the strata; see also very nice and accessible survey [E]. 

For a given square-tiled surface S denote by Aut(S) its automorphism 
group. Here we count only those automorphisms which isometrically send 
each square of the tiling to another square. For most of the square-tiled sur- 
faces Aut(S) is trivial; even for those rare square-tiled surfaces, which have 
nontrivial inner symmetries the group Aut(S) is obviously finite. We com- 
plete this section with the following arithmetic Theorem which confirms two 
conjectures of M. Kontsevich. 

Theorem (A. Eskin, A. Okounkov, R. Pandharipande). For every 
connected component of every stratum the generating function 



oo 



1 






N=l 



N-square-tiled 
surfaces S 



is a quasimodular form, i.e. a polynomial in Eisenstein series G2{q), G&(q), 
G 6 (q). 
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Volume Vol(Hj (di, . . . ,d m ) of every connected component of every 
stratum is a rational multiple of ir 2g , where g is the genus, d\ + • ■ ■ + d m = 
2g - 2. 

8 Crash Course in Teichmiiller Theory 

In this section we present the Teichmiiller theorem about extremal quasicon- 
formal map and define Teichmiiller metric. This enables us to explain finally 
in what sense the "Teichmiiller geodesic flow" (which we initially defined in 
terms of the action of the subgroup of diagonal matrices in SL(2,IR) on flat 
surfaces) is a geodesic flow. 

8.1 Extremal Quasiconformal Map 

Coefficient of Quasiconformality 

Consider a closed topological surface of genus g and two complex structures 
on it. Let So and Si be the corresponding Riemann surfaces. When the com- 
plex structures are different there are no conformal maps from So to Si. A 
smooth map f : So —> Si (or, being more precise, its derivative Df) sends an 
infinitesimal circle at x £ So to an infinitesimal ellipse, see Fig. 42. 




Fig. 42. Quasiconformal map 

Coefficient of quasiconformality of / at x G Sq is the ratio 

KM) = ~ 

of demi-axis of this ellipse. Coefficient of quasiconformality of / is 

K(f) = sup K x {f) 

Though So is a compact Riemann surface we use sup and not max since the 
smooth map / is allowed to have several isolated critical points where Df = 
(and not only det(Df) = 0) and where K x (f) is not defined. 
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Half-translation Structure 

There is a class of flat metrics which is slightly more general than the very flat 
metrics which we consider in this paper. Namely, we can allow to a flat metric 
to have the most simple nontrivial linear holonomy which is only possible: 
we can allow to a tangent vector to change its sign after a parallel transport 
along some closed loops (see the discussion of linear holonomy in Sec. 1.2). 

Surfaces endowed with such flat structures are called half-translation sur- 
faces. A holomorphic one- form (also called a holomorphic differential or an 
Abelian differential) is an analytic object representing a translation surface 
(in our terminology, a very flat surface). A holomorphic quadratic differential 
is an analytic object representing a half-translation surface. 

In local coordinate w a quadratic differential has the form q — q(w)(dw) 2 . 
In other words, the tensor rule for q is 



under a change of coordinate w = w(w'). 

One should not confuse (dw) 2 with a wedge product dw Adw which equals 
to zero! It is just a tensor of the type described by the tensor rule (29). In 
particular, any holomorphic one- form defined in local coordinates as uj = 
u>(w)dw canonically defines a quadratic differential uj 2 = uj 2 (w)(dw) 2 . 

Reciprocally, a holomorphic quadratic differential q = q(w)(dw) 2 locally 
defines a pair of holomorphic one forms ± \J q(w) dw in any simply-connected 
domain where q(w) ^ 0. However, for a generic holomorphic quadratic differ- 
ential neither of these 1-forms is globally defined: trying to extend the local 
form lo+ = q{uj) dw along a closed path we may get back with the form 
LO- = —y/q(w) dw. 

Recall that there is a bijection between very flat (^translation) surfaces 
and holomorphic 1-forms. There is a similar bijection between half-translation 
surfaces and holomorphic quadratic differentials, where similar to the "very 
flat" case a flat surface corresponding to a quadratic differential is polarized: 
it is endowed with canonical vertical and horizontal directions. (They can be 
defined locally using the holomorphic one-forms uj± = ±y // q(w) dw.) Note, 
however, that we cannot distinguish anymore between direction to the North 
and to the South, or between direction to the East and to the West unless 
the quadratic differential q is a global square of a holomorphic 1-form uj. In 
particular, the vertical and horizontal foliations are nonorientable for generic 
quadratic differentials. 

Teichmuller Theorem 

Choose any two complex structures on a topological surface of genus g > 1 ; 
let So and Si be the corresponding Riemann surfaces. Developing ideas of 
Grotzsch Teichmuller has proved a Theorem which we adapt to our language. 




(29) 
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Note that flat structure used in the formulation of the Theorem below is 
slightly more general than one considered in Sec. 1.2 and in Convention 1: it 
corresponds to a half-translation structure and to a holomorphic quadratic dif- 
ferential (see above in this section). In particular, speaking about a flat metric 
compatible with a given complex structure we mean a flat metric correspond- 
ing to a quadratic differential holomorphic in the given complex structure. 

Theorem (Teichmuller). For any pair So, S\ of Riemann surfaces of genus 
g > 1 there exist an extremal map fo : Sq — * S± which minimizes the coefficient 
of quasiconformality K(f). 

For this extremal map fo the coefficient of quasiconformality is constant 
on Sq ( outside of a finite collection of singular points ) 



One can choose a flat metric (half-translation structure) compatible with 
the complex structure in which foliation along big (correspondingly small) 
demi-axis of ellipses is the horizontal (correspondingly vertical) foliation in 
the flat metric. 

In flat coordinates the extremal map fo is just expansion- contraction with 
coefficient \f~K . 



Fig. 43. In flat coordinates the extremal map fo is just an expansion-contraction 
linear map 



8.2 Teichmuller Metric and Teichmuller Geodesic Flow 

Now everything is ready to define the Teichmuller metric. In this metric we 
measure the distance between two complex structures as 



where fo : Sq — > S± is the extremal map. 

It means that a holomorphic quadratic differential defines a direction of 
deformation of the complex structure and a geodesic in the Teichmuller metric. 



K x (fo) = K(fo) Vx £ S -{Px,..., P m } 




dist(S ,Si) = -logK(f ), 
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Namely, a holomorphic quadratic differential defines a flat metric. A one- 
parameter family of maps which in the flat coordinates are defined by diagonal 
matrices 



is a one-parameter family of extremal maps, so it forms a Tcichmuller geodesic. 
According to the definition above we have 



Remark. Note that the Teichmiiller metric is not a Ricmannian metric but 
a Finsler metric: it does not correspond to a quadratic form in the tangent 
space, but just to a norm which depends continuously on the point of the 
space of complex structures. 

It is known that the space of complex structures on a surface of genus 
g > 2 has complex dimension 3<7 — 3. We have seen, that the space of pairs 
(complex structure, holomorphic quadratic differential) can be identified with 
a tangent space to the space of complex structures, in particular, it has com- 
plex dimension 617 — 6. Taking into consideration the functorial behavior of the 
space of pairs (complex structure, holomorphic quadratic differential) one can 
check, that it should be identified with a total space of a cotangent bundle. 

9 Hope for a Magic Wand and Recent Results 

This section is devoted to one of the most challenging problems in the theory 
of flat surfaces: to the problem of complete classification of the closures of 
all orbits of Gi + (2,M) on the moduli spaces of Abelian (and quadratic) dif- 
ferentials. This problem was very recently solved for genus two in the works 
of K. Calta and of C. McMullen; we give a short survey of their results in 
Sec. 9.7 and 9.8. 

9.1 Complex Geodesies 

In this section we are following the geometric approach of C. McMullen de- 
veloped in [McM2] and [McM3]. 

Fix the genus g of the surfaces. We have seen in the previous section 
that we can identify the space Q of pairs (complex structure, holomorphic 
quadratic differential) with the total space of the (co)tangent bundle to the 
moduli space M. of complex structures. Space H of pairs (complex structure, 
holomorphic quadratic differential) can be identified with a subspace in Q 
of those quadratic differentials, which can be represented as global squares 
of holomorphic 1-forms. This subspace forms a vector subbundle of special 
directions in the (co)tangent space which we denote by the same symbol TL. 




dist(S ,g t So) = t. 
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The "unit hypcrboloid" 7t l c TL of holomorphic 1-forms corresponding to flat 
surfaces of unit area can be considered as a subbundle of unit vectors in "hi- 
lt is invariant under the Teichmiiller geodesic flow - the geodesic flow for the 
Teichmiiller metric. 

One can check that an SL(2; M)-orbit in Hi descends to a commutative 
diagram 

SX(2;R) > Hi 

(30) 



SL(2;R)/SO(2;R) 
which we interpret as 



Unit tangent 
bundle to 
.hyperbolic plane , 



Hyperbolic 
plane 



M, 



(Unit tangent ^ 
subbundle to 
moduli spacey 



/ Moduli 
I space 



Moreover, it can be checked that the map H 2 — > Ai in this diagram is 
an isometry with respect to the standard hyperbolic metric on H 2 and Te- 
ichmiiller metric on M.. Thus, following C. McMullcn it is natural to call the 
projections of SL(2; R)-orbits to M. (which coincide with the images of H 2 ) 
complex geodesies. Another name for these projections is Teichmiiller discs. 



9.2 Geometric Counterparts of Ratner's Theorem 

Though it is proved that the moduli space of complex structures is not a 
hyperbolic manifold (see [Mai]) there is a strong hope that with respect to 
SL(2, R)-action on TL and on Q the moduli space behaves as if it is. 

In this section we present several facts about group actions on homoge- 
neous spaces and several related facts about geodesic submanifolds. We warn 
the reader that our selection is not representative; it opens only a narrow slit 
to the fascinating world of interactions of group actions, rigidity, hyperbolic 
geometry, dynamics and number theory. 

We start with an informal formulation of part of Ratner's Theorem (see 
much better exposition adopted to our subject in the survey of A. Eskin [E]). 

A discrete subgroup T of a Lie group G is called a lattice if a homogeneous 
space G/r has finite volume. 

Theorem (M. Rattier). Let G be a connected Lie group and U a connected 
subgroup generated by unipotent elements. Then, for any lattice T C G and 
any x € G/T the closure of the orbit Ux in G/T is an orbit of some closed 
algebraic subgroup of G. 
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We would like to point out why this theorem is so remarkably powerful. 
Considering a dynamical system, even an ergodic one, it is possible to get a 
lot of information about a generic (in measure-theoretical sense) trajectory. 
However, usually there arc plenty of trajectories having rather particular be- 
havior. It is sufficient to consider geodesic flow on a surface with cusps to find 
trajectories with closures producing fairly wild sets. Ratner's theorem proves, 
that the closure of any orbit of a unipotcnt group acting on a homogeneous 
space is a nice homogeneous space. 

Ratner's theorem has numerous important relations with geometry of ho- 
mogeneous spaces. As an illustration we have chosen a result of N. Shah [Sh] 
and a generalization of his result for noncompact hyperbolic manifolds ob- 
tained by T. Payne [Pa]. 

Theorem (N. Shah). In a compact manifold of constant negative curva- 
ture, the closure of a totally geodesic, complete (immersed) submanifold of 
dimension at least 2 is a totally geodesic immersed submanifold. 

Theorem (T. Payne). Let f : M\ — ► M2 be a totally geodesic immersion 
between locally symmetric spaces of noncompact type, with Mi of finite volume. 
Then the closure of the image of f is an immersed submanifold. Moreover, 
when Mi and M2 have the same rank, the closure of the image is a totally 
geodesic submanifold. 



9.3 Main Hope 

Main Conjecture and its Possible Applications 

If only the moduli space of complex structures M. would be a homogeneous 
space we would immediately apply the Theorem above to diagram (30) and 
would solve considerable part of our problems. But it is not. Nevertheless, 
there is a strong hope for an analogous Theorem. 

Problem. Classify the closures ofGL + (2,M.)-orbits in7i g and in Q g . Classify 
the orbit closures of the unipotent subgroup ^ on 7i g and on Q g . 

The following Conjecture is one of the key conjectures in this area. 

Conjecture. The closure C(S) of a GL + (2,M.)-orbit of any flat surface S £ TL 
(or S € Q) is a complex-algebraic suborbifold. 

Remark. Wc do not discuss here the problems related with possible compact- 
ifications of the moduli spaces Tl g and Q g . A complex-analytic description 
of a compactification of Q g can be found in the papers of J. Fay [Fay] and 
H. Masur [Ma2]. 
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Recall that according to Theorem of M. Kontsevich (see Sec. 3.6) any 
GL + (2, R)-invariant complex suborbifold in 7i is represented by an affine sub- 
space in cohomological coordinates. Thus, if the Conjecture above is true, the 
structure of orbit closures of the action of GL + (2, R) on Tt and on Q (and of 
SL(2, R) on Hi and on Qi) would be as simple as in the case of homogeneous 
spaces. 

We have not discussed the aspects of Ratncr's Theorem concerning the 
measures; it states more than we cited above. Actually, not only orbit clo- 
sures have a nice form, but also invariant ergodic measures; namely, all such 
measures are just the natural measures supported on orbits of closed sub- 
groups. Trying to make a parallel with Ratner's Theorem one should extend 
the Conjecture above to invariant measures. 

In the most optimistic hopes the study of an individual flat surface 
S G H(di, . . . ,d m ) would look as follows. (Frankly speaking, here we slightly 
exaggerate in our scenario, but after all we are describing the dreams.) To 
describe all geometric properties of a flat surface S we first find the orbit clo- 
sure C(S) = GL+(2,R) S C H(di, . . . , d m ) (our optimistic hope assumes that 
there is an efficient way to do this). Then we consult a (conjectural) classifica- 
tion list and find C (S) in some magic table which gives all information about 
C(S) (like volume, description of cusps, Siegel-Veech constants, Lyapunov ex- 
ponents, adjacency to other invariant subspaces, etc). Using this information 
we get answers to all possible questions which one can ask about the initial 
flat surface S. 

Billiards in rational polygons give an example of possible implementa- 
tion of this optimistic scenario. Fixing angles of the polygon which defines 
a billiard table we can change the lengths of its sides. We get a family B 
of polygons which induces a family B of flat surfaces obtained by Katok- 
Zcmlyakov construction (see Sec. 2.1). This family B belongs to some fixed 
stratum B C Ti(di, . . . 7 d m ). However, it has a nontrivial codimension in the 
stratum, so B has measure zero and one cannot use ergodic theorem naively 
to get any information about billiards in corresponding polygons. Having a 
version of ergodic theorem which treats all orbits (like in Ratner's Theorem) 
presumably it would be possible to get a powerful tool for the study of rational 
billiards. 

Exercise. Consider the family B of billiard tables as on Fig. 37. Determine 
the stratum H.(di, . . . , d m ) to which belong the corresponding flat surfaces 
and compute the codimension of the resulting family B C TL{d\ 1 . . . , d m ). 

Content of Remaining Sections 

The Conjecture above is trivial for genus one, since in this case the "Te- 
ichmullcr space of Riemann surfaces of genus one" coincides with an upper 
half-plane, and the entire space coincides with a single Tcichmuller disc (image 
of H 2 in diagram (30)). 
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Very recently C. McMullcn proved the Conjecture in genus two [McM2], 
[McM3], and this is a highly nontrivial result. We give a short report of rev- 
olutionary results of K. Calta [Clt] and of C. McMullcn [McM2]-[McM6] in 
Sec. 9.7 below. However, their techniques, do not allow any straightforward 
generalizations to higher genera: Ricmann surfaces of genus two are rather 
special, in particular, every such surface is hyperelliptic. 

In the next two sections 9.4 and 9.5 we try to give an idea of what is known 
about invariant submanifolds in higher genera (which is an easy task since, 
unfortunately, little is known). 

Having an invariant submanifold K, C TL g (or K, C Q g ) in genus g one can 
construct a new invariant submanifold K, C TLg (correspondingly /C C Qj) in 
higher genus g > g replacing every S € K. by an appropriate ramified covering 
S over S of some fixed combinatorial type. We do not want to specify here 
what does a "fixed combinatorial type" mean; what we claim is that having 
an invariant manifold /C there is some procedure which allows to construct a 
whole bunch of new invariant submanifolds K, for higher genera g > g. 

In some cases all quadratic differentials in the invariant submanifold JC 
obtained by a ramified covering construction from some K, C Q g might become 
global squares of Abelian differentials. Hence, using special ramified coverings 
one can construct GL + (2, R)-invariant submanifolds K, C TLg from invariant 
submanifolds /C C Q g . 

What is really interesting to understand is what invariant manifolds form 
the "roots" of such constructions. Such invariant manifolds are often called 
the primitive ones. 

In the following two sections 9.4 and 9.5 we consider the two extremal 
classes of primitive invariant submanifolds: the largest ones and the small- 
est ones. Namely, in Sec. 9.4 we present a classification of connected com- 
ponents of the strata TC(di, . . . , d m ). It follows from crgodicity of SL(2, En- 
action on TLi(d±, . . . , d m ) that the orbit closure of almost any flat surface 
in TL(d±, . . . , d m ) coincides with the embodying connected component of 
H(di, . . . , d m ). 

In section 9.5 we consider the smallest possible GL + (2, R)-invariant sub- 
manifolds: those which correspond to closed orbits. Teichmuller discs obtained 
as projections of such orbits to the moduli space M. of complex structures form 
the "closed complex geodesies" — Riemann surfaces with cusps. 

9.4 Classification of Connected Components of the Strata 

In order to formulate the classification theorem for connected components of 
the strata TL(di, . . . , d m ) we need to describe the classifying invariants. There 
are two of them: spin structure and hyperellipticity. Both notions are appli- 
cable only to part of the strata: flat surfaces from the strata TL(2di, . . . , 2d m ) 
have even or odd spin structure. The strata TL(2g — 2) and TL{g — 1, g — 1) have 
special hyperelliptic connected component. 
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Spin Structure 

Consider a flat surface S from a stratum TC(2di, . . . ,2d m ). Let p : S 1 — > S 
be a smooth closed path on S; here S 1 is a standard circle. Note that at 
any point of the surfaces S we know where is the "direction to the North" . 
Hence, at any point p(t) — x G S we can apply a compass and measure the 
direction of the tangent vector x. Moving along our path p(t) we make the 
tangent vector turn in the compass. Thus we get a map G{p) : S 1 — ► S 1 from 
the parameter circle to the circumference of the compass. This map is called 
the Gauss map. We define the index ind(p) of the path p as a degree of the 
corresponding Gauss map (or, in other words as the algebraic number of turns 
of the tangent vector around the compass) taken modulo 2. 

ind(p) = deg G(p) mod 2 

It is easy to see that ind(p) does not depend on parameterization. More- 
over, it does not change under small deformations of the path. Deforming 
the path more drastically we may change its position with respect to conical 
singularities of the flat metric. Say, the initial path might go on the left of Pk 
and its deformation might pass on the right of This deformation changes 
the degG(p). However, if the cone angle at Pk is of the type 2ir(2dk + 1), 
then degG(p) mod 2 does not change! This observation explains why ind(p) 
is well-defined for a free homotopy class [p] when S S H(2di, . . . , 2d m ) (and 
hence, when all cone angles arc odd multiples of 2ir) . 

Consider a collection of closed smooth paths ai, b\, . . . , a g , b g representing 
a symplectic basis of homology Hi(S, Z/2Z). We define the parity of the spin- 
structure of a flat surface S £ H(2di, . . . , 2d m ) as 

9 

<j)(S) = (indK) + 1) (indOi) + 1) mod 2 (31) 

i=l 

Lemma. The value (friS) does not depend on the choice of symplectic basis 
of cycles {ai,bi}. It does not change under continuous deformations of S in 
W(2di,...,2d m ). 

Lemma above shows that the parity of the spin structure is an invariant 
of connected components of strata of those Abelian differentials, which have 
zeroes of even degrees. 

Exercise. Consider two flat surfaces presented at Fig. 44. They are obtained 
by a surgery which attaches a handle to a flat surface obtained from a regular 
octagon. Note, however, that the handles are attached in two different ways 
(see the identifications of vertical sides). Check that both surfaces belong to 
the same stratum W(4). 

Consider a symplectic basis of cycles of the initial surface (corresponding 
to the regular octagon) realized by paths which do not pass through the 
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Fig. 44. Attaching a handle to the fiat surface So £ in two different ways 

we get two flat surfaces in H(4) with different parities of spin structure. Hence the 
resulting flat surfaces live in different connected components of 7i(4) 

conical singularity. Show that a symplectic basis of cycles for each of two new 
surfaces can be obtained by completion of the initial basis with a pair of cycles 
a, b representing the attached handle, where the cycle a is the waist curve of 
the handle. Calculate ind(a) and ind(&) for each of two surfaces. Check that 
ind(b) are different, and thus our two flat surfaces have different parities of 
spin structure and hence belong to different connected components of 7i(4). 

Spin structure: more serious reading. We have hidden under the carpet geom- 
etry of the "spin structure" defining the "parity-of-the-spin-structure" . The 
reader can find details in [KonZo] and in original papers of M.Atiyah [At], 
J. Milnor [Mil], D.Mumford [Mum] and D. Johnson [J]. Recent paper of C. Mc- 
Mullen [McM4] contains further applications of spin structures to fiat surfaces. 

Hyperellipticity 

A flat surface S may have a symmetry; one specific family of such fiat sur- 
faces, which are "more symmetric than others" is of a special interest for 
us. Recall that there is a one-to-one correspondence between flat surfaces 
and pairs (Ricmann surface AI , holomorphic 1-form uj), see Sec. 3.3. When 
the corresponding Riemann surface is hyperelliptic the hyperelliptic involution 
t : M — > M acts on any holomorphic 1-form ui as t*lj = —uj. 

We say that a flat surface S is a hyperelliptic flat surface if there is an 
isomctry r : S — > S such that r is an involution, r o r = id, and the quotient 
surface S/t is a topological sphere. In flat coordinates differential of such 
involution obviously satisfies Dt = — Id. 

Exercise. Check that the flat surface S from Fig. 12 is hyperelliptic, and that 
the central symmetry of the polygon induces the hyperelliptic involution of S. 

In a general stratum H(di, . . . , d m ) hyperelliptic surfaces form a small 
subspace of nontrivial codimension. However, there are two special strata, 
namely H(2g — 2) and Ti(g — l,g — 1), for which hyperelliptic surfaces form 
entire hyperelliptic connected components H hvp (2g — 2) and TL hyp (g— l,g— 1) 
correspondingly. 
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Note that in the stratum TC(g — l,g—l) there are hypcrelliptic flat surfaces 
of two different types. A hypcrelliptic involution tS — > S may fix the conical 
points or might interchange them. It is not difficult to show that for surfaces 
from the connected component Tt hyp (g — 1, g — 1) the hypcrelliptic involution 
interchanges the conical singularities. 

The remaining family of those hyperelliptic flat surfaces in H(g — 1, g — 1), 
for which the hypcrelliptic involution keeps the saddle points fixed, forms a 
subspace of nontrivial codimension in the complement H(g—1, g—l)—H vp (g— 
1, g— 1). Thus, the hyperelliptic connected component Ti hyp (g — 1, g— 1) does 
not coincide with the space of all hyperelliptic flat surfaces. 

Classification Theorem for Abelian Differentials 

Now, having introduced the classifying invariants we can present the classifi- 
cation of connected components of strata of Abelian differentials. 

Theorem (M. Kontsevich and A. Zorich). All connected components of 
any stratum of Abelian differentials on a curve of genus g > 4 are described 
by the following list: 

The stratum H.{2g — 2) has three connected components: the hyperelliptic 
one, Tt hyp (2g — 2), and two nonhyperelliptic components: 7i even (2g — 2) and 
H. odd (2g — 2) corresponding to even and odd spin structures. 

The stratum H.(2d, 2d), d > 2 has three connected components: the hyperel- 
liptic one, TL hyp (2d, 2d), and two nonhyperelliptic components: H. even (2d,2d) 
and H odd (2d, 2d). 

All the other strata of the form TL(2di, . . . , 2d m ) have two connected com- 
ponents: 7i even (2di, . . . , 2d m ) and TL odd {2di, . . . , 2d„), corresponding to even 
and odd spin structures. 

The stratum TL{2d — \,2d — 1), d > 2, has two connected components; one 
of them: H hyp (2d- 1, 2d- 1) is hyperelliptic; the other H nonhyp (2d-l, 2d- 1) 
is not. 

All the other strata of Abelian differentials on the curves of genera g > 4 
are nonempty and connected. 

In the case of small genera 1 < g < 3 some components are missing in 
comparison with the general case. 

Theorem. The moduli space of Abelian differentials on a curve of genus g = 2 
contains two strata: 7i(l, 1) andTL{2). Each of them is connected and coincides 
with its hyperelliptic component. 

Each of the strata Ti(2, 2), 7i(4) of the moduli space of Abelian differentials 
on a curve of genus g = 3 has two connected components: the hyperelliptic one, 
and one having odd spin structure. The other strata are connected for genus 
.9 = 3. 

Since there is a one-to-one correspondence between connected components 
of the strata and extended Rauzy classes (see Sec. 5.6 and paper [Y] in this 
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collection) the Classification Theorem above classifies also the extended Rauzy 
classes. 

Classification Theorem for Quadratic Differentials 

Note that for any partition d± + ■ • • + d m = 2g — 2 of a positive even integer 
2g — 2 the stratum TC(di, . . . ,d m ) of Abelian differentials is nonempty For 
meromorphic quadratic differentials with at most simple poles there are four 
empty strata! Namely. 

Theorem (H. Masur and J. Smillie). Consider a partition of the number 
Ag — A, where g > into integers di + • • • + d m with all dj S N U {— 1}. 
The corresponding stratum Q(d±, . . . , d m ) is non-empty with the following four 
exceptions: 

Q(0), 2(1,-1) {in genus g = 1) and 2(4), 2(1, 3) (m ^enus 5 = 2) 

Classification of connected components of the strata of meromorphic 
quadratic differentials with at most simple poles was recently obtained by 
E. Lanneau [Lan]. 

Theorem (E. Lanneau). Four exceptional strata Q(— 1,9), Q(— 1,3, 6), 
Q(— 1,3,3, 3) and 2(12) of meromorphic quadratic differentials contain ex- 
actly two connected components; none of them hyper elliptic. 
Three series of strata 

Q(2(g - k) - 3, 2{g - k) - 3, 2k + 1, 2k + 1) k > -1, 5 > 1, g - k > 2 
Q(2(s - fc) - 3,2(5 - fc) - 3,4fc + 2) fc>0, 5>1 and g - k > 1 

Q(4(5 - fe) - 6, 4fc + 2) k>0, g>2andg-k>2 

contain hyperelliptic connected components. The strata from these series in 
genera g > 3 and the strata Q(— 1, — 1, 3, 3), Q(— 1,-1,6) in genus g = 2 
contain exactly two connected components; one of them - hyperelliptic, the 
other one - not. 

The remaining strata from these series, namely, 2(1,1,1,1), 2(1,1,2), 
2(2,2) in genus g = 2 and 2(1,1,-1,-1), Q(— 1,-1,2) in genus g = 1 
coincide with their hyperelliptic connected component. 

All other strata of meromorphic quadratic differentials with at most simple 
poles are connected. 

Recall that having a meromorphic quadratic differential with at most sim- 
ple poles one can associate to it a surface with a flat metric which is slightly 
more general that our usual very flat metric (see Sec. 8.1 for a discussion of 
half-translation structures). 

It is easy to verify whether a half-translation surface belongs to a hyper- 
elliptic component or not. However, currently there is no simple and efficient 
way to distinguish half-translation surfaces from the four exceptional compo- 
nents. 
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Problem. Find an invariant of the half-translation structure which would 
be easy to evaluate and which would distinguish half-translation surfaces 
from different connected components of the four exceptional strata Q{— 1,9), 
e(-l,3,6), e(-l,3,3,3) and Q(12). 

Currently there are two ways to determine to which of the two connected 
components of an exceptional stratum belongs a flat surface S. 

The first approach suggests to find a "generalized permutation" for an 
analog of the first return map of the vertical flow to a horizontal segment 
and then to find it in one of the two extended Rauzy classes (see Sec. 5.6) 
corresponding to two connected components. Note, however, that already for 
the stratum Q(— 1,9) the corresponding Rauzy classes contain 97 544, and 
12 978 generalized permutations; the Rauzy classes for the components of 
Q(12) contain already 894 117 and 150 457 elements. 

In the second approach one studies configurations of saddle connections 
(see Sec. 6.4) on the surface S and tries to find a configuration which is forbid- 
den for one of the two connected components of the corresponding stratum. 

For example, for surfaces from one of the two connected components of 
Q(— 1, 9) as soon as we have a saddle connection joining the simple pole with 
the zero we necessarily have a closed geodesic going in the same direction. 
Thus, if we manage to find on a surface S £ Q(— 1,9) a saddle connection 
which is not accompanied by a parallel closed geodesic, S belong to the other 
connected component of Q(— 1,9). 

9.5 Veech Surfaces 

For almost every fiat surface S in any stratum . . . , d m ) the orbit 

5L(2,R) • S is dense in the stratum and for any g\ ^ g 2 € SX(2,R) we 
have giS ^ g 2 S. However, some flat surfaces have extra symmetries. When a 
flat surface So has an affine automorphism, i.e. when for some go £ SX(2, R) 
we get goS = S the orbit of So is smaller than usual. 

The stabilizer Stab(S) £ <SL(2;R), that is a subgroup of those g £ 
SX(2,R) for which g S = S 7 is called the Veech group of the flat surface S 
and is denoted SL(S). In representation of the flat surface S in terms of a 
pair (Riemann surface X, holomorphic 1-form ui on it) the Veech group is 
denoted as SL(X,u>) following the notation of C. McMullen [McM2]. 

Some exceptional flat surfaces S possess very large group of symmetry 
and their orbits are very small. The flat surfaces having the largest possible 
symmetry group are called Veech surfaces. More precisely, a flat surface is 
called a Veech surface if its Veech group SL(S) is a lattice in SX(2, R) (that 
is the quotient SL(2,M.)/SL(S) has finite volume). 

Theorem (J. Smillie). An SL(2,M.)-orbit of a flat surface S is closed in 
Tii(di, . . . , d m ) if and only if S is a Veech surface. 
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Forgetting polarization (direction to the "North" ) of a flat surface we get 
a Teichmiiller disc of S (see (30) and the comments below it) 

50(2, R) ^ SL ( 2 ' R)/ SL(S) = H/ SL(S) 

A flat surface S is a Veech surface if its Teichmiiller disc M 2 /SL(S) has 
finite volume. However, even for a Veech surface the Teichmiiller disc is never 
compact: it necessarily contains at least one cusp. The Teichmiiller discs of 
Veech surfaces can be considered as closed complex geodesies (see discussion 
at the end of Sec. 9.1). 

Consider an elementary example. As a flat surface take a flat torus ob- 
tained from a unit square. Fig. 45 shows why the clement g + = ^ 
SX(2,Z) belongs to a stabilizer of S. 




Similarly the clement g_ = ^ G SL(2, R) also belongs to the Veech 

group SL(T 2 ) of T 2 . Since the group SL(2, Z) is generated by g+ and g- 
we conclude that SL(2,Z) C SL(T 2 ). It is easy to check that, actually, 
SL(2, 1) = SX(T 2 ). As the Teichmiiller disc of T 2 we get the modular curve 
H 2 /SX(2,Z) (see Fig. 13) which, actually, coincides with the moduli space of 
complex structures on the torus. 
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Fig. 46. There are three 3-square-tiled surfaces in TL(2). Our picture shows that 
they all belong to the same 51/(2; Z)-orbit 
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Consider a slightly more complicated example. 

Exercise. Verify that square-tiled surfaces presented at Fig. 46 belong to the 
stratum TL{2). Show that there are no other 3-square-tiled surfaces. Verify 
that the linear transformations indicated at Fig. 46 act as it is described on 
the Figure; check that the surfaces belong to the same SL(2, R)-orbit. Find 
Veech groups of these three surfaces. Show that these fiat surfaces are Vcech 
surfaces. Verify that the corresponding Tcichmuller disc is a triple cover over 
the modular curve (see Fig. 47). 




/ 



\ 



_L 



Fig. 47. Teichmiiller discs of a 3-square-tiled surface is 
modular curve 



triple cover over the 



Primitive Veech Surfaces 

It is not difficult to generalize the Exercise above and to show that any square- 
tiled surface (see Sec. 7.1) is necessarily a Veech surface. 

A square-tiled surface is a ramified covering over a fiat torus, such that 
all ramification points project to the same point on the fiat torus, which is a 
Veech surface. One can generalize this observation. Having a Veech surface S 
one can construct a ramified covering S — > S such that all ramification points 
on S project to conical singularities on S. One can check that any such S is a 
Vcech surface. Thus, having a Veech surface we can construct a whole bunch 
of Veech surfaces in higher genera. 

Veech surfaces which cannot be obtained from simpler Veech surfaces by 
the covering construction are called primitive. For a long time (and till recent 
revolution in genus two, see Sec. 9.7, the list of known primitive Veech surfaces 
was very short. Very recently C. McMullen has found infinitely many Veech 
surfaces in genera 3 and 4 as well, see [McM7] . All other known primitive Veech 
surfaces of genus g > 2 can be obtained by Katok-Zemlyakov construction 
(see Sec. 2.1) from triangular billiards of the first three types in the list below: 
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for n > 6 (discovered by W. Veech) 

for n > 7 (discovered by W. Veech) 

for n > 4 (discovered by Ya. Vorobets) 

(discovered by W. Veech)) 

(discovered by Ya. Vorobets) 

(discovered by R. Kenyon 

and J. Smillie) 

(discovered by W. P. Hooper) 

- The fiat surface corresponding to the isosceles triangle with the an- 
gles 7r/n, (n— l)7r/(2n), (n— l)ir/{2n) belongs to the hypcrclliptic component 
n hyp (2g-2) when n = 2g and to the hyperelliptic component H hyp (g-1, g-1) 
when n = 2g + 1. The surface can be unwrapped to the regular 2n-gon with 
opposite sides identified by parallel translations, see Fig. 21. 

- The flat surface corresponding to the isosceles triangle with the angles 
n/n, n/n, (n — 2)n/(2n) belongs to the hyperelliptic component H hyp (2g — 2) 
when n = 2g + 1 and to the hyperelliptic component Ti hyp (g — 1, g — 1) when 
n = 2g + 2. The surface can be unwrapped to a pair of regular n-gons glued 
by one side. Each side of one polygon is identified by a parallel translation 
with the corresponding side of the other polygon, see Fig. 7. 

- The flat surface corresponding to the obtuse triangle with the angles 
7r/n, 7r/(2n), (2n— 3)n/(2n) belongs to one of two nonhypcrelliptic components 
of the stratum TC(2g — 2) where n = g + 1. 

- The flat surface corresponding to the acute triangle 7r/3, 7r/4, 57r/12 
belongs to the nonhypcrelliptic component H odd (4); here g = 3. 

- The flat surface corresponding to the acute triangle 7r/3, 7r/5, 77r/15 
belongs to the nonhyperelliptic component H. even (6); here 5 = 4. 

- The flat surface corresponding to the acute triangle 2ir/9, 3ir/9, Att/9 
belongs to the stratum H(3, 1); here .9 = 3. 

- The flat surface corresponding to the obtuse triangle tt/3, 7r/12, 77r/12 
belongs to the stratum 7i(6); here g = 4 (the information that this is a Veech 
surface is taken from [McM7]). 

The details on unwrapping of these surfaces and on cylinder decomposi- 
tions of some of them can be found in the paper of Ya. Vorobets [Vbl]. 
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It is proved that unwrapping triangular billiards in other acute, rectan- 
gular or isosceles triangles does not give new Veech surfaces in genera g > 2 
(see [KcnS], [Pu], [Vbl] and further references in these papers). For obtuse 
triangles the question is open. 

We discuss genus g — 2 separately in the next section: very recently 
K. Calta [Clt] and C. McMullen [McM2] have found a countable family of 
primitive Veech surfaces in the stratum 7i(2) and proved that the list is com- 
plete. However, even in genus g = 2 the situation with the stratum H(l, 1) is 
drastically different: using the results of M. Moeller [Mol]-[Mo3] very recently 
C. McMullen has proved [McM6] the following result. 

Theorem (C. McMullen). The only primitive Veech surface in the stra- 
tum H(l, 1) is the surface represented by the regular decagon with identified 
opposite sides. 

Thus, it is not clear, what one should expect as a solution of the following 
general problem. 

Problem. Find all primitive Veech surfaces. 

An algebro-geometric approach to Veech surfaces suggested by M. Moller 
in [Mol] and [Mo2] might help to shed some light on this Problem. 

Veech surfaces: more serious reading. I recommend the survey paper [HuSdt5] 
of P. Hubert and T. Schmidt as an introduction to Veech surfaces. A canonical 
reference for square-tiled surfaces (also called arithmetic Veech surfaces) is the 
paper of E. Gutkin and C. Judge [GuJg]. More information about Veech sur- 
faces can be found in the pioneering paper of W. Veech [Ve7] and in the paper 
of Ya. Vorobets [Vbl]. For the most recent results concerning Veech groups 
and geometry of the Tcichmullcr discs see the original papers of P. Hubert and 
T. Schmidt [HuSdtl], [HuSdt2], [HuSdt3], [HuSdt4], of C. McMullen [McMl], 
[McM7] and of P. Hubert and S. Lelievre [HuLel], [HuLc2] . 

9.6 Kernel Foliation 

In this section we describe some natural holomorphic foliation on the moduli 
space of Abelian differentials. In higher genera little is known about this foli- 
ation (though it seems to be worth of study) . We use this foliation in the next 
section to describe GL(2, R)-invariant submanifolds of "intermediate type" 
discovered by K. Calta and by C. McMullen in genus two. 

We have seen that any stratum Ti(d\, . . . , d m ) can be locally parameterized 
by a collection of basic relative periods of the holomorphic one-form w, or, in 
other words, that a neighborhood U([ui}) C H 1 (S; {Pi, . . . , P m }', C) gives a 
local chart in 7i(di, . . . , d m ). 

Let S € 7i(l,l). Let closed paths 01,02,61,62 represent a basis of cycles 
in Hi(S]7*). Any path c joining conical singularities P\ and Pi represents a 
relative cycle in Hi(S, {Pi, P 2 }; Z). Let Ai,A 2 ,Bi 7 B 2 ,C € C be the periods 
of cu: the integrals of oj over a\, a 2 , 61, 62, c correspondingly. 
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Fig. 48. A deformation of a fiat surface inside the kernel foliation keeps the absolute 
periods unchanged 



Example. The collection of cycles a,, &», i = 1, 2, on the surfaces from Fig. 48 
represent a basis of cycles in Hi(S; Z). All horizontal geodesies on these sur- 
faces are closed; each surface can be decomposed into three cylinders filled 
with horizontal geodesies. Let wi 7 w 2 ,W3 = w\ + w 2 be the widths (perime- 
ters) of these cylinders; h\,h2,hs be their heights; tx,t-z,tz their twists, see 
Fig. 39. It is easy to check that 

A\ = I to = —w\ 

J ai 



A 2 = oj = -{wi + w 2 ) 

J a-2 




The kernel foliation in 1) is the foliation defined in local coordinates by 
equations 

{Ai = constn J B\ = const\2 

A 2 = const 2 i )B 2 = const 22 

In other words, this is a foliation which is obtained by fixing all absolute peri- 
ods and changing the relative period C = jp^ lo. Similarly, the kernel foliation 
in arbitrary stratum TL{d\, . . . , d m ) is a foliation which in cohomological coor- 
dinates is represented by parallel complex (to— l)-dimcnsional affinc subspaces 
obtained by changing all relative periods while fixing the absolute ones. 

Passing to a finite cover over H(di, . . . , d m ) we can assume that all zeroes 
Pi, ... , P m are named (i.e. having two zeroes Pj, of same degrees, we know 
exactly which of the two is Pj and which is Pk). Now we can fix an arbitrary 
subcollection of zeroes and define a kernel "subfoliation" along relative periods 
corresponding to chosen subcollection. 

Recall that the area of a flat surface is expressed in terms of the absolute 
periods (sec Riemann bilinear relation in Table 1 in Sec. 3.3). Thus, moving 
along leaves of kernel foliation we do not change the area of the surface. In 



Flat Surfaces 117 



particular, we can consider the kernel foliation as a foliation of the "unit 
hyperboloid" Hi(di, . . . , d m ). 

Exercises on Kernel Foliation 

Exercise. To deform the flat surface on the left of Fig. 48 along the kernel 
foliation we have to keep all Ai, A 2 , B\, P>2 unchanged. Hence, we cannot 
change the widths (perimeters) of the cylinders, since they are expressed in 
terms of A\ and A 2 . Increasing the height of the second cylinder by e we 
have to increase the height of the first cylinder by the same amount e to 
keep B\ = (ti — t 2 ) + V—l((h 2 + e ) — {hi + s)) unchanged; we also have to 
decrease the height /13 of the third cylinder by e to preserve the value of B 2 . 
Similarly, increasing the twist t\ by 6 we have to increase the twist t 2 by the 
same amount 5 and to decrease t% by 5. 

Exercise. It is convenient to consider the kernel foliation in the total moduli 
space Ti. g of all holomorphic 1-forms without subdivision of Ti. g into strata 
Ti.(di, . . . , d m ), where J2j dj = 2.9 — 2. In particular, to deform a surface S £ 
~H(2) C TL 2 along the kernel foliation we have to break the double zero into 
two simple zeroes preserving the absolute periods. The corresponding surgery 
is presented at Fig. 31. 

The leaves of the kernel foliation are naturally endowed with a flat struc- 
ture, which has conical singularities at the points of intersection of the leaf 
with the smaller strata and with degenerate strata. 

Assuming that the zeroes Pi, P 2 of a surface S <E 7i(l, 1) are named show 
that the intersection of the kernel foliation with the stratum H(2) corresponds 
to a conical point with the cone angle 6ir, while the intersections with the two 
strata of degenerate flat surfaces (determine which ones) are just the regular 
points of the flat structure. 

In the exercise below we use a polygonal representation of a flat surface 
(compare to Fig. 8 in the paper [Clt] of K. Calta). 

Exercise. Consider a regular decagon. Imagine that there are springs inside 
its sides so that we can shrink or expand the sides keeping them straight 
segments. Imagine that we hammer a nail in the center of each side. Though 
the centers of the sides are now fixed our decagon is still flexible: we can pull a 
vertex and the whole frame will follow, see Fig. 49b. We assume that under any 
such deformation each nail stays exactly in the middle of the corresponding 
side. 

Prove that the deformed polygon is again centrally symmetric with the 
same center of symmetry. Prove that the opposite sides of the deformed poly- 
gon are parallel and have equal length. Prove that the resulting flat surface 
lives in the same leaf of the kernel foliation. Show that the "nails" (the cen- 
ters of the sides) and the center of symmetry are the Weierstrass points of the 
corresponding Ricmann surface (fixed points of the hypcrclliptic involution). 
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Fig. 49. This cartoon movie represents a path living inside a leaf of the kernel 
foliation. When at the stage c it lands to a surface from the stratum TC(2) we 
remove the nails from one pair of vertices and hammer them in the other symmetric 
pair of vertices (stage d). Then we continue the trip inside the kernel foliation 

Move a vertex towards the center of the side which this vertex bounds. 
The side becomes short (see the upper side on Fig. 49b) and finally contracts 
to a point (Fig. 49c). What flat surface do we get? 

Show that the deformation presented at Fig. 49c brings us to a surface 
So in the stratum Ti(2). Remove the pair of clues, which is hammered at the 
vertices of the resulting octagon. Hammer them to another pair of symmetric 
vertices (see Fig. 49d). We can declare that we have a new centrally-symmetric 
decagon with a pair of sides of zero lengths. Stretching this pair of sides and 
making them have positive length (see Fig. 49e) we continue our trip inside 
the kernel foliation. Show that for a given small value v € C of the saddle 
connection joining two conical singularities, there are exactly three different 
surfaces obtained as a small deformation of the surface So € W(2) as on 
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Fig. 49c along the leaf of the kernel foliation and having a saddle connection 
v. (Use previous Exercise and Fig. 31.) 

Following our path in the kernel foliation we get to the surface Si as 
on Fig. 49f. Is this surface singular? To what stratum (singular stratum) it 
belongs? 

Compact Leaves of the Kernel Foliation 

Let us show that the leaf JC of the kernel foliation passing through a square- 
tiled surface S(u>q) of genus two is a compact square-tiled surface. To simplify 
notations suppose that the absolute periods of the holomorphic 1-form ujq, 
representing the flat surface S{uiq), generate the entire integer lattice Z + 



Consider the "relative period" map p : JC — > T 2 from the corresponding leaf 
JC of the kernel foliation containing S(u>o) to the torus. The map p associates 
to a flat surface S(u>) £ JC the relative period C taken modulo integers, 



K,3S(uj)Ac= lu modZ + V^TZ £ C/(Z + \^-f Z) = T 2 



Since the flat surface S(uj) belongs to the same leaf JC as the square-tiled 
surface S(u>o) € JC, the absolute periods of u> are the same as the ones of 
wo, an d hence the integral above taken modulo integers does not depend on 
the path on S(u>) joining Pi and P%. It is easy to check that the map p is a 
finite ramified covering over the torus T 2 , and thus the leaf JC is a square-tiled 
surface. 

Those flat surfaces S(u>) £ JC, which have integer relative period C £ 
Z + yj—ll,, have all periods in Z + \J — 1 Z. Hence, these flat surfaces are 
square-tiled. Since S(ui) and S(u)q) have the same area, the number N of 
squares tiling S(u>) and S(uq) is the same. Thus, JC has a structure of a 
square-tiled surface such that the vertices of the tiling are represented by 
jV-square-tiled surfaces S(u) £ JC. 

To discuss the geometry of JC we need to agree about enumeration of zeroes 
P\,Pi of a surface S £ 7~t(l, 1). We choose the convention where the zeroes 
are named. That is, given two zeroes of order 1 we know which of them is Pi 
and which of them is P^. Under this convention the square-tiled surface JC is 
a translation surface; it is represented by a holomorphic one-form. (Accepting 
the other convention we would obtain the quotient of JC over the natural 
involution exchanging the names of the zeroes. In this latter case the zeroes 
of S £ H(l, 1) arc not distinguishable; the leaf of the kernel foliation gives a 
flat surface represented by a quadratic differential.) 

The lattice points of the square-tiled surface JC are represented by N- 
squarc-tiled surfaces S £ JC of several types. We have the lattice points repre- 
sented by iV-square tiled surfaces from 7i(l, 1). These points are the regular 
points of the flat metric on JC. 




Z. 




Pi 
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There arc points of intersection of JC with Tt(2). Such point S(u>) € K,fXH{2) 
is always represented by a square-tiled surface, and hence gives a vertex of 
the tiling of JC. We have seen (see Fig. 31) that given a surface S{w) € H(2) 
and a small complex period C one can construct three different flat surfaces 
S(u>x), S(u)2), S(u>3) e H(l, 1) with the same absolute periods as to and with 
the relative period C . Thus, the points S(u>) 6 /CflW(2) correspond to conical 
points of JC with the cone angles 3 ■ 2ir when the zeroes are named (and with 
the cone angle 37r, when they are not named). The total number of such points 
was computed in the paper of A. Eskin, H. Masur and M. Schmol [EMaScm]; 
it equals 



3 / 2 t — r ( 1 N 

number of conical points on JC = — (N — 2)N II - 

p\N \ P / 



(32) 



There remain vertices of the tiling of JC represented by degenerate square- 
tiled surfaces S(ui). It is not difficult to show that these points are regular 
for the flat metric on JC(co) when the zeroes are named. (They correspond to 
conical singularities with the cone angle 7r, when the zeroes are not named). 
The degenerate iV-square-tiled surfaces S(u>) are of two types. It might be an 
-ZV-square-tilcd torus with two points of the tiling identified. It might be a pair 
of square-tiled tori with a vertex of the tiling on one torus identified with a 
vertex of the tiling on the other torus. Here the total number of squares used 
to tile these two tori is N. The total number of the vertices of the tiling of JC 
of this type is computed in the paper [Schll] ; it equals 

number of special points on JC = ^(5^ + 6)iV 2 J^J |l J 

p\N \ P ) 

Summarizing we conclude that the translation surface JC lives in the stra- 
tum H(2, . . . , 2), where the number k of conical points is given by formula (32). 

k 

We complete this section with an interpretation of a compact leaf K. 
as a space of torus coverings; this interpretation was introduced by A. Es- 
kin, H. Masur and M. Schmol in [EMaScm] and developed by M. Schmoll 
in [Schll], [Schl2]. 

We have seen that a nondegencrate flat surface S(u>q) £ 1) repre- 
senting a vertex of the square tiling of K, is an A^-square-tiled surface. Hence, 
S(u>) is a ramified covering over the standard torus T 2 of degree N having two 
simple ramification points, which project to the same point of the torus. A 
non vertex point S(u>) G K, is also a ramified covering over the standard torus 
T 2 . To see this consider once more the period map, but this time applied to 



■ r p — 

S{lo)3P p ^ uj modZ + v^lZ G C/(Z + v /r lZ) 



T 2 . 
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Since all absolute periods of ui live in Z + Z the integral taken modulo 
integers does not depend on the path joining the marked point (conical singu- 
larity) Pi with a point P of the flat surface S(u>). The map prof is a ramified 
covering. 

The degree of the covering can be computed as the ratio of areas of S(uj) 
and of the torus T 2 , which gives N. The covering has precisely two simple 
ramification points, which are the conical points P\,Pi of S(u>). This time 
they project to two different points of the torus. 

Recall that by convention we assume that the absolute periods of us gen- 
erate the entire lattice Z + Z. They corresponds to primitive covers: the 
ones which do not quotient through a larger torus. 

Proposition (M. Schmoll). Consider primitive branched covers over the 
standard torus T 2 . Fix the degree N of the cover. Let the cover have exactly 
two simple branch points. 

The space of such covers is connected; its natural compactification coin- 
cides with the corresponding leaf JC of the kernel foliation. The Veech group 
of the square-tiled surface JC coincides with SL(2,Z). 

Connectedness of the space of covers is not quite obvious (actually, it was 
proved earlier in other terms by W. Fulton [Ful]). An observation above shows, 
that the leaf JC coincides with a connected component of the space of covers. 
Thus, connectedness of the space of covers implies that this space coincides 
with JC. The group SL(2, Z) acts naturally on the space of covers; in particular 
it maps the space of covers to itself. This implies that SL(2, Z) belongs to the 
Veech group of the square-tiled surface JC. It is easy to show, that it actually 
coincides with SX(2,Z). 

Corollary (M. Schmoll). Consider square-tiled surfaces S (to) of genus two 
such that the absolute periods of us span the entire integer lattice Z + \/ — 1 Z. 
For any given N > 3 all such N-square tiled surfaces belong to the same 
compact connected leaf JC(N) of the kernel foliation. 

For more information on kernel foliation of square-tiled surfaces in genus 
two see the papers of A. Eskin, H. Masur and M. Schmol [EMaScm] and of 
M. Schmoll [Schll], [Schl2]. In particular, the latter papers propose a beautiful 
formula for Siegel- Veech constants of any flat surface S € JC(N) in terms of 
geometry of the cylinder decomposition of the square-tiled surface JC(N). 

9.7 Revolution in Genus Two (after K. Calta and C. McMullen) 

In this section we give an informal survey of recent revolutionary results in 
genus .9 = 2 due to K. Calta [Clt] and to C. McMullen [McM2]. 

Using different methods they found a countable collection of primitive 
Veech surfaces in the stratum H(2), proved that this collection describes 
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all Veech surfaces, and gave efficient algorithms which recognize and classify 
Veech surfaces in W(2). 

This result is in a sharp contrast with the Theorem of C. McMullen [McM6] 
cited above, which tells that in the other stratum Tt(l, 1) in genus g = 2 there 
is only one primitive Veech surface. 

This discovery of an infinite family of primitive Veech surfaces in the stra- 
tum 7i(2) is also in a sharp contrast with our poor knowledge of primitive 
Veech surfaces in higher genera: as we have seen in the previous section, 
primitive Veech surfaces in higher genera g > 3 are currently known only in 
some special strata (mostly hyperelliptic) , and even in these special strata we 
know only finite number of primitive Veech surfaces (basically, only one). 

Another remarkable result is a discovery by K. Calta and by C. McMullen 
of nontrivial examples of invariant submanifolds of intermediate dimension: 
larger than closed orbits and smaller than the entire stratum. 

One more revolutionary result in genus two is a Classification Theorem 
due to C. McMullen [McM3] which proves that a closure of any GL + (2,R)- 
orbit is a nice complex-analytic variety which is either an entire stratum, or 
which has one of the types mentioned above. 

Algebro- geometric Approach To avoid overloading of this survey I had to sacri- 
fice beautiful algebro-geometric part of this story developed by C. McMullen; 
the reader is addressed to original papers [McM2]-[McM6] and to a short 
overview presented in [HuSdt5] . 

Periods of Veech Surfaces in Genus g = 2 

If S is a Veech surface then the fiat surface gS is also a Veech surface for any 
g € GL + (2,R). Thus, speaking about a finite or about a countable collection 
of Veech surfaces we, actually, choose some family of representatives {Sk} of 
the orbits GL+(2,R) • S of Veech surfaces. 

The question, which elements of our collection {Sk} belong to the same 
GL + (2, IR)-orbit and which ones belong to different orbits is a matter of a 
separate nontrivial study. A solution was found by C. McMullen in [McM4]; 
it is briefly presented in the next Sec. 9.8. In this section we present effective 
algorithm due to K. Calta and to C. McMullen which enables to determine 
whether a given flat surface in H(2) is a Veech surface. 

Following K. Calta we say that a flat surface S can be rescaled to a flat 
surface S" if S and S' belong to the same GL + (2, R)-orbit. We say that a flat 
surface S is quadratic if for any homology cycle c G Hi(S; Z) we have 



In other words, we say that a flat surface S defined by a holomorphic 1-form 



We can considerably restrict the area of our search using the following 
Lemma of W. Thurston. 




eQ 
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Lemma (W. Thurston). Any Veech surface in genus g = 2 (no matter 
primitive or not, in the stratumTi(2) or Ti(l, 1)) can be rescaled to a quadratic 
surface. 

Using this Lemma K. Calta suggest the following algorithm deciding 
whether a given flat surface S <E H(2) is a Veech surface or not. 

Algorithm of Calta 

Recall, that if S is a Veech surface (of arbitrary genus), then by Veech alter- 
native (see Sec. 3.7) a directional flow in any direction is either minimal or 
completely periodic, that is a presence of a closed geodesic going in some direc- 
tion implies that all geodesies going in this direction are periodic. Moreover, it 
was proved by Veech that as soon as there a saddle connection going in some 
direction, this direction is also completely periodic. In both cases the surface 
decomposes into a finite collection of cylinders; each boundary component of 
each cylinder contains a conical singularity (see Sec. 7.1). 

The algorithm works as follows. Having a flat surface S € H(2) it is easy to 
find some closed geodesic on S (which is allowed to be a closed geodesic saddle 
connection). Since "rescaling" the surface S (i.e. applying a linear transfor- 
mation from GL + (2,R)) we replace a Veech surface by a Veech surface, we 
can turn S in such way that the direction of the closed geodesic will become 
horizontal. We denote the resulting surface by the same symbol S. 

Since S lives in H.(2) it has a single conical point P with the cone angle 
6ir. In particular, there are exactly three geodesies leaving the conical point 
in the positive horizontal direction (to the East). If at least one of these three 
horizontal geodesies does not come back to P the surface 5" is not a Veech 
surface. Otherwise our test continues. 

As we have seen in Sec. 7.1 there are two possible ways in which three 
horizontal geodesies emitted from P to the East can return to P. Either 
all three geodesies return at the angle 37r, or one of them returns at the 
angle 3ir and two others return at the angle 7r, see Fig. 41. In both cases all 
horizontal geodesies are closed. In the first case the surface decomposes into 
a single cylinder; in the second case the surface is glued from two cylinders, 
see Sec. 7.1. 

If the surface is decomposed into a single cylinder, it is sufficient to compare 
the lengths pi,p2,P3 of three horizontal saddle connections, see Fig. 41. The 
flat surface S is a Veech surface if and only if pi,P2,P3 are commensurable. 
Moreover, if pi,P2,P3 are commensurable, we can rescale S to a square-tiled 
surface. It can be done in several elementary steps. First we rescale S in 
the horizontal direction making p\,p2,P3 rational and then integer. Then we 
rescale 5* in the vertical direction making the height h of the cylinder integer. 



Finally, we apply an appropriate parabolic linear transformation . It 

does not change neither p\ 1 p2 1 P3 nor h, but when s € R varies the twist t 
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also varies continuously (sec Fig. 39) taking all values in R; in particular, we 
can achieve t — 0. We get a square-tiled surface. 

Consider the second case, when the surface decomposes into two cylinders. 
Let wx,W2,hx,h%,ti, ti be the widths (perimeters), heights and twists of this 
cylinders correspondingly (compare to Sec. 7.1). In a complete analogy with 
the one-cylinder case we can rescale the surface horizontally, then vertically, 
and finally apply an appropriate parabolic linear transformation in order to 
make the width (perimeter) w\ and height hi of the first cylinder equal to 
one, hi = W\ = 1, and the twist t\ = equal to zero. Applying an appropriate 
Dehn twist to the second cylinder we can assure < t 2 < w-2- 

If after our rescaling all parameters 102,^2,^2 characterizing the second 
cylinder do not get to the same quadratic field Q(Vd) for some d £ N, the 
surface S is not a Veech surface. 

If u>2,/i2,i2 are rational (i.e. if d is a complete square), the surface S can 
be rescaled to a square-tiled surface. 

The remaining case is treated by one of the key Theorems in the paper of 
K. Calta [Clt]. 

Theorem (K. Calta). Let all parameters Wj,hj,tj, j = 1,2 of a two- 
cylinder decomposition of a flat surface S £ H(2) belong to the same quadratic 
field <Q(\/d) with d £ N not a complete square. Then S is a Veech surface if 
and only if the parameters satisfy the following system of equations: 



Actually, we kept the system of equation above as it is written in the origi- 
nal paper [Clt] . In this form it can be adopted to a more general normalization 
of parameters: it is sufficient to rescale surface 5" to bring all Wj, hj, tj, j = 1, 2 
to a quadratic field. 

Remark. Similar necessary conditions for Veech surfaces in Ti(2) were ob- 
tained by D. Panov independently of K. Calta and of C. McMullcn. 

Since by Lemma of Thurston any Veech surface in Ti (2) can be rescaled to a 
quadratic surface, taking a collection of all quadratic surfaces decomposed into 
two horizontal cylinders satisfying the condition above, we get representatives 
of the GL+(2, R)-orbits of all flat surfaces in H(2). 

Exercise. Show that the Katok-Zcmlyakov construction applied to an L- 
shaped billiard as on Fig. 50 (see also Fig. 38) generates a surface S £ H(2). 
Show that this surface is decomposed into two cylinders filled by closed hor- 
izontal geodesies and that these cylinders have parameters w± = 2, hi = 




(33) 
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2a — 2, t\ = for the first cylinder and u>2 = 2b, h 2 = 2, t<i = for the 
second cylinder. Using the condition above prove the following Theorem of 
C. McMullen [McM2]: 



Fig. 50. L-shaped billiard table P(a, b) and its unfolding into a flat surface S £ H(2) 
(after C. McMullen [McM2]) 



Theorem ( C. McMullen). The L-shaped billiard table P(a, b) as on Fig. 50 
generates a Veech surface if and only if a and b are rational or 

a = x + zVd and b = y + z\fd 

for some x,y,z€Q with x + y = 1 and d > in Z. 

Kernel Foliation in Genus 2 

The following elementary observation explains our interest to kernel foliation 
in the content of our study of GL + (2, R)-invariant subvarieties. Let J\f C 
H(di, . . . ,d n ) be a GL + (2, R)-invariant submanifold. The "germ" of the kernel 
foliation at M is equivariant with respect to the action of GL + (2,R). 

In other words this statement can be described as follows. Denote by t(S) 
a translation by S along kernel foliation defined in a neighborhood of S G 
W(l,l). Here i5 £ C is a small parameter. Let g £ GL + (2,R) be close to 
identity. Then 

g o t(S) ■ S = t(gS) o gS 

where g acts on a complex number 6 as on a vector in R 2 . Moving along the 
kernel foliation, and then applying an element g of the group is the same 
as applying first the same element of the group and then moving along an 
appropriate translation along the kernel foliation. A similar construction works 
in a general stratum. 
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Wc get the following very tempting picture. Suppose, that we have a closed 
GI/ + (2, K)-orbit Af in H(l, 1) or in TC(2). For example, suppose that Af is an 
orbit of a Veech surface. Consider a union of leaves of the kernel foliation 
passing through Af. Due to the remark above this is a GL + (2, K)-invariant 
subset! 

The weakness of this optimistic picture is that there is no a priori reason 
to hope that the resulting invariant subset in H(l, 1) would be closed. And 
here the magic comes. The following statement was proved independently by 
K. Calta [Clt] and by C. McMullen [McM2]. 

Theorem (K. Calta; C. McMullen). For any Veech surface So £ H(2) the 
union of leaves of the kernel foliation passing through the GL + (2, M.)- orbit of 
Sq is a closed GL + (2, R)- invariant complex orbifold Af of complex dimension 
3. 

In particular, the complex dimension of the resulting orbifold is the sum 
of the complex dimension of the GL + (2, K)-orbit of So and of the complex 
dimension of the kernel foliation 3 = 2 + 1. 

K. Calta has found the following beautiful geometric characterization of 
flat surfaces living in an invariant subvaricty Af as above. Surfaces in any 
such Af are completely periodic: as soon as there is a single closed trajectory 
in some direction, all geodesies going in this direction are closed. This condi- 
tion is necessary and sufficient condition for a surface to live in an invariant 
subvaricty Af as above. In particular, this shows that not only Veech surfaces 
have this property. 

An algorithm analogous to the algorithm determining Veech surfaces in 
H(2) (see above) allows to K. Calta to determine whether a given surface 
S G H(l,l) is completely periodic or not (and hence, whether it belongs to 
an invariant subvariety Af as above or not). As before one starts with finding 
some closed geodesic or some saddle connection. If the surface is completely 
periodic in the corresponding direction, then all other geodesies going in this 
direction are periodic and the surface decomposes into cylinders. After an 
appropriate rotation this periodic direction becomes horizontal. Without loss 
of generality, we may assume that S decomposes into three cylinders. By un, 
hi and t, with 1 < i < 3, wc denote the widths, heights and twists. After 
renumbering, we may assume that u>3 = W\ + w-i- Define s\ = h\ + /13, 
S2 = hi + /13, Ti = t\ + £3, T2 = <2 + £3- If the surface is completely periodic its 
absolute periods can be rescaled to get to Q(Vd) + iQ(Vd) (compare to the 
algorithm for Veech surfaces). Leaving the elementary case when d £ N is a 
complete square, the following characteristic equations obtained by K. Calta 
(analogous to equations (33) above) tell whether our flat surface is completely 
periodic or not: 

WlSl = —W2S2, 

win + W2T2 = wif\ + W2T2, < n < Wi + W3. 
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Note that we consider any invariant subvariety TV as above as a subvariety 
in H 2 = H(l, 1) U H(2). Clearly the intersection TV n H(l, 1) and TV n H(2) 
results in close GL + (2, R)-invariant subvaricties in the corresponding strata. 
Note that dimcAf = 3. Since TV is a union of leaves of the kernel foliation this 
implies that dimcAf n H(2) = 2. This means that any surface S G TVn H(2) 
is a Veech surface! 

Classification Theorem of McMullen 

We complete this section with a description of the wonderful result of C. Mc- 
Mullen [McM3] realizing a dream of a complete classification of closures of 
GL + (2, R)-orbits in genus g = 2. The classification is astonishingly simple. 
We slightly reformulate the original Theorem using the notions of kernel foli- 
ation and of completely periodic surface (and, hence, using implicitly results 
of K. Calta [Clt]). 

Theorem (C. McMullen). 

• If a surface S € 7~t(2) is a Veech surface, its GL + (2,M)-orbit is a closed 
complex 2-dimensional subvariety; 

• Closure of GL + (2,WL)-orbit of any surface S G H(2) which is not a Veech 
surface is the entire stratum T~t(2); 

• If a surface S G H(l, 1) is a Veech surface, its GL + (2,M.)-orbit is a closed 
complex 2-dimensional subvariety; 

• If a surface S G 7i(l, 1) is not a Veech surface but is a completely peri- 
odic surface, then the closure of its GL + (2,WL)-orbit is a closed complex 
3- dimensional subvariety TV foliated by leaves of the kernel foliation as 
described above; 

• If a surface S £ W(l, 1) is not completely periodic, then the closure of its 
GL + (2,WL)-orbit is the entire stratum H. (1,1). 

Actually, the Theorem above is even stronger: connected components of 
these invariant submanifolds are basically also classified. We have seen that 
the GL + (2, R)-orbit of any Veech surfaces in TC(2) has a representative with 
all periods in a quadratic field. The discriminant D = b 2 — 4c > is a positive 
integer: a discriminant of the corresponding quadratic equation x 2 + bx + c = 
with integer coefficients. The discriminant is an invariant of an GL + (2,M)- 
orbit. Since for any integer b the number b 2 mod 4 can be either or 1, the 
discriminant D mod 4 = 0, 1. The values D = 1,4 are not realizable, so the 
possible values of D are 5, 8, 9, 12, 13, ... . 

We postpone the description of results of P. Hubert and S. Lelievre [HuLel] 
and of C. McMullen [McM4] on classification of the orbits of Veech surfaces 
in Tt(2) to the next section. Here we state the following result of C. Mc- 
Mullen [McM3]. By Af(D) denote the 3-dimcnsional invariant submanifold 
obtained as a union of leaves of the kernel foliation passing through the orbits 
of all Veech surfaces corresponding to the given discriminant D. 
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Theorem (C. McMullen). The invariant subvariety J\f(D) is nonempty 
and connected for any D = 0, 1 mod 4, D G N, D > 5. 

Ergodic Measures 

Actually, the Classification Theorem of C. McMullen is even stronger: it also 
classifies the invariant measures. Consider the "unit hyperboloids" Hi (2) and 
7ii(l, 1): the subvarieties of real codimcnsion one representing flat surfaces of 
area 1. The group SL(2, R) acts on 7ii(2) and 1) preserving the measure 

induced on these "unit hyperboloids", see Sec. 3.4. Let us discuss what other 
SL(2, R)-invariant measures do we know. 

When we have an SL ( 2, R) -invariant subvariety, we can get an invari- 
ant measure concentrated on this subvariety. For example, as we know an 
SL(2, R)-orbit of a Veech surface S is closed; it is isomorphic to the quotient 
SL(2, R)/r , (5 l ), where, by definition of a Veech surface, this quotient has finite 
volume. Thus, Haar measure on iSX(2,R) induces a finite invariant measure 
on the SL(2, R)-orbit of a Veech surface S. 

Consider now a "unit hyperboloid" J\f\ C N in the manifold N obtained 
as a union of leafs of the kernel foliation passing through SL(2, R)-orbit of a 
Veech surface S. Note that by Ricmann bilinear relations the area of a flat 
surface can be expressed in terms of absolute periods. Thus, moving along the 
kernel foliation we do not change the area of the surface. We have seen that 
every leaf of the kernel foliation is flat. Consider the corresponding Euclidean 
volume element in each leaf. The group 5L(2,R) maps leaves of the kernel 
foliation to leaves and respects this volume element. Thus we get an invariant 
measure on Afi C Af; near an SL(2, R) -orbit of a Veech surface it disintegrates 
to a product measure. 

We have associated to any connected invariant subvariety of each of four 
types as above a natural SL ( 2, R) -measure supported on it. One more result 
of C. McMullen in [McM3] tells that there are no other ergodic measures. 

Other Properties 

The invariant subvarieties have numerous wonderful geometric properties. In 
particular, their projections to the moduli space M of complex structures on 
a surface of genus two are also nice subvarieties. C. McMullen has showed that 
the GL + (2, R)-orbit of a Veech surface projects to an isometrically immersed 
algebraic curve and J\f(D) projects to a complex surface. Such surfaces (of 
complex dimension two) are called Hilbert modular surfaces. 

One more surprising phenomenon proved by C. McMullen in [McM3] con- 
cerns Veech groups of flat surfaces in genus g = 2. 

Theorem (C. McMullen) . If the Veech group T{S) of a flat surface S G 
TL{2) contains a hyperbolic element, the flat surface S is a Veech surface; in 
particular, its Veech group T(S) is a lattice in SX(2,R). 
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If the Veech group r(S) of a flat surface S € 7i(l, 1) contains a hyperbolic 
element, and the flat surface S is not a Veech surface, then S is completely 
periodic. In this case the Veech group r(S) is infinitely generated. 

The Veech group of a flat surface S contains a hyperbolic element if and 
only if S admits an affine pseudoanosov diffeomorphism. 

We complete this section with the following natural problem for higher 
genera g > 3. 

Problem. Let JC be a GL + (2, M.) -invariant subvariety in some stratum of 
holomorphic one-forms Ti(d\, ■ ■ ■ , d m ) C Ti g . Consider the unionU of leaves of 
the kernel foliation passing through K,. IsU a closed subvariety in 7i(l, ■ ■ ■ , 1) ? 
Similar question for other strata. 

9.8 Classification of Teichmiiller Discs of Veech Surfaces in T~C(2) 

It is easy to check that any square-tiled surface (see Sec. 7.1) is a Veech 
surface, and thus an SL(2, R)-orbit of any square-tiled surface is closed. Such 
orbit contains other square-tiled surfaces. Since the SL{2, R)-action does not 
change the area of a surface these other square-tiled surfaces are tiled with 
the same number of squares, see Fig. 46 in Sec. 9.5 and the Exercise related 
to this Figure. 

For a fixed integer n the number of n-square-tiled surfaces is finite. It would 
be interesting to know (and this is certainly a part of general Problem from 
Sec. 9.3) how the square-tiled surfaces are arranged into orbits of SL{2, K). 
Say, we have seen in the previous section that there are exactly three 3- 
squarc-tiled surfaces in Ti(2) (see Fig. 46) and that they belong to the same 
orbit. For n = 4 there are already nine 4-square-tiled surfaces in H(2) and 
they still belong to the same SL{2, M)-orbit. The corresponding Teichmiiller 
disc is a 9-fold cover over the modular curve. For n = 5 there are twenty 
seven 5-square-tiled surfaces in Ti(2) and they split into two different orbits 
of SL(2,R). 

Generalizing a result of P. Hubert and S. Lclicvre [HuLel] obtained for 
prime number n C. McMullen has recently proved the following conjecture of 
P. Hubert and S. Lelievre. 

Theorem (C. McMullen). All n- square-tiled surfaces inH(2), which can- 
not be tiled with p x q-rectangles with p or q greater than 1, get to the same 
SL(2,M.)-orbit when n > 4 is even and get to exactly two distinct orbits when 
n > 5 is odd. 

Actually, C. McMullen has classified in [McM4] the orbits of all Veech 
surfaces in H(2). As we have seen in the previous section Veech surfaces in 
Tt(2) are characterized by an integer parameter, called the discriminant D. 
For n-square-tiled surfaces the discriminant equals D = n 2 . 
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Veech surfaces which cannot be rescaled to a square-tiled surface are called 
nonarithmetic Veech surfaces. Any nonarithmctic Veech surface in Tl(2) can 
be rescaled to a flat surface having all periods in a quadratic field (see the 
Lemma of W. Thurston in the previous section). The discriminant correspond- 
ing to a nonarithmetic Veech surface is the discriminant of this quadratic field. 
Of course a GL + (2, R)-orbit of a nonarithmetic Veech surface might have dif- 
ferent representatives S, such that all periods of S belong a quadratic field. 
Nevertheless, for Veech surfaces in genus g — 2 the discriminant is well-defined: 
it is an invariant of a GL + (2, K)-orbit. The discriminant is a positive integer 
D = 0, 1 mod 4, D > 5. 

C. McMullen has proved the following classification Theorem [McM4]: 

Theorem ( C. McMullen) . For D = 1 mod 8, D > 9, all Veech surfaces in 
Tl(2) corresponding to discriminant D get to exactly two distinct GL + (2,IR)- 
orbits. For other values D = 0, 1 mod 4, D > 5, they belong to the same 
GL+(2,R) -orbit. 



Fig. 51. The L-shaped billiard table L(b, e) generating a "canonical" Veech surface 



Moreover, C. McMullen proposed the following canonical representative for 
any such GL + (2, R)-orbit, see [McM4]. Consider an L-shaped billiard L(b,e) 
as on Fig. 51, where 



'6,eeZ; A= (e + Ve 2 + 46)/2 
e = -l,0orl 

e+Kb v ' 

if e = 1 then b is even 

The billiard table L(b, e) generates a flat surface S(b, e) in H(2), see Fig. 50. 

Theorem (C. McMullen). The flat surface 5(6, e) generated by the L- 
shaped billiard table L(b, e) with parameters 6, e satisfying (34) is a Veech 
surface. Any closed GL + ' (2, M.)- orbit in TL{2) is represented by one of such 
S(b, e) and this representation is unique. The discriminant D of 5(6, e) equals 
D = e 2 + 46. 
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Exercise. Using linear transformations and scissors rescale the fiat surface 
obtained from the "double pentagon" (sec Fig. 7) to the flat surface obtained 

/ \ 1 + 

from a "golden cross" P[a,b), with a = b = — - — ; see Fig. 50 for the 

definition of P(a,b). Using linear transformations and scissors rescale any of 
these flat surfaces to a surface obtained from the billiard table L(l, — 1), see 
Fig. 51, proving that these Veech surfaces have discriminant D = 5. (For 
solutions see Fig. 4 in [McM2]; see also and [McM5]). 

Exercise (T. Schmidt). Using linear transformations and scissors rescale the 
flat surface obtained from the regular octagon to a surface obtained from 
the billiard table £(2,0), see Fig. 51, proving that these Veech surfaces have 
discriminant D = 8. (See [McM5]). 

Stratum H(l,l) 

We have discussed in details Veech surfaces in Tl(2). We did not discussed 
the Veech surface for the stratum H(l,l) because for square-tiled surfaces 
in W(l, 1) (also called arithmetic Veech surfaces) the classification of the 
GL + (2, K)-orbits is knot known yet... 

The classification of nonarithmetic Veech surfaces in 7i(l, 1) is, however, 
known (since very recently), and is quite surprising. Using the results of 
M. Moeller [Mol]-[Mo3] C. McMullcn has proved in [McM6] the following 
Theorem. 

Theorem (C. McMullen). Up to a reseating the only primitive nonarith- 
metic Veech surface in 1) is the one obtained from the regular decagon by 
identification of opposite sides. In other words, any primitive nonarithmetic 
Veech surface inH(l, 1) belongs to the GL + (2,M.)-orbit of the surface obtained 
from the regular decagon. 

For higher genera nothing is known neither about the number of SL(2, R)- 
orbits of 7i-square-tiled surfaces, nor about their geometry. 

Problem. Classify orbits of square-tiled surfaces in any stratum, in particular 
inW(l,l). 

Square-tiled Surfaces: more serious reading. An elementary introduction can 
be found in [Zo5]. Paper [HuLel] of P. Hubert and S.Lelievre and [McM4] of 
C. McMullcn classify orbits of square-tiled surfaces in TC(2). See also the paper 
of G. Schmithiiscn [Schn] for an algorithm of evaluation of the Veech group of a 
square-tiled surface and for examples of square-tiled surfaces having SL(2, R) 
as a Veech group. Another such example due to M. Moller is presented in the 
survey [HuSdt5] of P. Hubert and T. Schmidt. 
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10 Open Problems 

Flat surfaces with nontrivial holonomy and billiards in general 
polygons 

Problem 1 (Geodesies on general flat surfaces; Sec. 1.1). 

Describe behavior of geodesies on general flat surfaces with nontrivial 
holonomy. Prove (or disprove) that geodesic flow is ergodic on a typical (in 
any reasonable sense) flat surface. 

Does any (almost any) flat surface has at least one closed geodesic which 
does not pass through singular points? 

If yes, are there many regular closed geodesies? Namely, find the asymp- 
totics for the number of closed geodesies of bounded length as a function of 
the bound. 

Problem 2 (Billiards in general polygons; Sec. 2.1). 

Describe the behavior of a generic regular billiard trajectory in a generic 
triangle, in particular, prove (or disprove) that the billiard flow is ergodic. 

Does any (almost any) billiard table has at least one regular periodic 
trajectory? If the answer is affirmative, does this trajectory survive under 
deformations of the billiard table? 

If a periodic trajectory exists, are there many periodic trajectories like 
that? Namely, find the asymptotics for the number of periodic trajectories of 
bounded length as a function of the bound. 

More problems on billiards can be found in the survey of E. Gutkin [Gu2]. 

Problem 3 (Renormalization of billiards in polygons; Sec. 2.1 and Sec. 5). 

Is there a natural dynamical system (renormalization procedure) acting 
on the space of billiards in polygons? 

Classification of orbit closures in H g and Q g 

Problem 4 (Orbit closures for moduli spaces; Sec. 9.3). 

Is it true that the closures of GL + (2, K)-orbits in Jig and Q g are always 
complex-analytic (complex-algebraic?) orbifolds? Classify these closures. Clas- 
sify ergodic measures for the action of SX(2,R) on "unit hyperboloids" . 

Suppose that these orbit closures are described by an explicit list. Find 
natural intrinsic invariants of a flat surface S which would allow to determine 
the closure of the orbit of S in the list. 

To be honest, even having obtained a conjectural classification above, one 
would need to develop a serious further machinery to get full variety of inter- 
esting applications. The situation with the problem below is quite different: 
a reasonable solution of this problem would immediately give a burst of ap- 
plications since such a machinery already exists. For the experts interested 
in ergodic aspects, counting problems, etc, the measure-theoretic analogue of 
Ratner's Theorem discussed below is the biggest open problem in the area. 
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Problem 5 (A. Eskin: "Ratner's theorem" for moduli spaces; Sec. 9.3). 

Does the unipotent subgroup of 51,(2, R) act nicely on Tt g and Q g or not? 
Are the orbit closures always nice (for example, real-analytic) orbifolds or one 
can get complicated closures (say, Kantor sets)? 

If the action is "nice" , classify the closures of orbits of the unipotent sub- 
group ( ,1 !) in TL g and Q g . Classify crgodic measures for the action of the 

unipotent subgroup on "unit hyperboloids" . 
Solve the problem for genus g = 2. 

The problem above is solved in the particular case when the unipotent 
flow acts on a SL(2, R)-invariant submanifold in TL g obtained by a ramified 
covering construction from a Veech surface; see the papers of A. Eskin, H. Ma- 
sur and M. Schmoll [EMaScm] and of A. Eskin, J. Marklof, D. Witte Mor- 
ris [EMkWt]. 

The next problem concerns possibility of construction of Gi + (2,R)- 
invariant submanifolds in higher genera using kernel foliation. 

It follows from the results of K. Calta and C. McMullen that for any Veech 
surface So £ 7~l{2) the union of complex one-dimensional leaves of the kernel 
foliation passing through the complex two-dimensional GL + {2, R)-orbit O(Sq) 
of So is a closed complex orbifold Af of complex dimension 3, see Sec. 9.7. By 
construction it is GL + (2, R)-invariant. 

Problem 6 (Kernel foliation; Sec. 9.7). 

Let O C H(di, . . . , d m ) C TL g be a GL(2,R)-invariant submanifold (sub- 
orbifold). Let TL(d' 1: . . . ,d' n ) C Ti g be a bigger stratum adjacent to the first 
one, dj = d' k + ■ ■ ■ + d' k . , i = 1, . . . , m. 

Consider the closure of the union of leaves of the kernel foliation in 
H(d[,...,d' n ) (or in H g ) passing through O. We get a closed Gi+(2,R)- 
invariant subset Af C H(d[, . . . ,d' n ) (correspondingly Af C Ti. g ). 

Is Af a complex-analytic (complex-algebraic) orbifold? When Af does not 
coincide with the entire connected component of the stratum W(di, . . . , d' n ) 
(correspondingly Tl g )7 When dime Af = dime O + (n — m) (correspondingly 
dime Af = dime O + (2g — 2 — m)l Here (n — m) and (2g — 2 — m) is the 
complex dimension of leaves of the kernel foliation in Ti.(d[, . . . , d' n ) and in TL g 
correspondingly. 

Particular cases of classification problem 

Problem 7 (Exceptional strata of quadratic differentials; Sec. 9.4). 

Find an invariant which would be easy to evaluate and which would dis- 
tinguish half-translation surfaces from different connected components of the 
four exceptional strata Q(-l,9), Q(-l,3,6), Q(-l,3,3,3) and Q(12). 

Problem 8 (Veech surfaces; Sec. 9.5). 
Classify primitive Veech surfaces. 
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Problem 9 (Orbits of square-tiled surfaces; Sec 9.8). 

Classify orbits of square-tiled surfaces in any stratum. Describe their Te- 
ichmuller discs. 

Same problem for the particular case 7i(l, 1). 

Geometry of individual flat surfaces 

Problem 10 (Quadratic asymptotics for any surface; Sec. 6.1). 

Is it true that any very flat surface has exact quadratic asymptotics for the 
number of saddle connections and for the number of regular closed geodesies? 

Problem 11 (Error term for counting functions; Sec. 6.1). 

What can be said about the error term in quadratic asymptotics for count- 
ing functions N(S, L) ~ c ■ L? on a generic flat surface SI In particular, is it 
true that the limit 

log|iV(5,i)- C -L 2 | 7 
Jim sup - < 2 

L~,oo logi 

is strictly less than two? Is it the same for almost all flat surfaces in a given 
connected component of a stratum? 

One of the key properties used by C. McMullcn for the classification of the 
closures of orbits of GL(2,M.) in Tt(l, 1) was the knowledge that on any flat 
surface in this stratum one can find a pair of homologous saddle connections. 
Cutting the surface along these saddle connections one decomposes the surface 
into two tori and applies machinery of Ratncr theorem. 

Problem 12 (A. Eskin; C. McMullen: Decomposition of surfaces; Sec. 6.4). 

Given a connected component of the stratum Tl(d\, . . . ,d m ) of Abelian 
differentials (or of quadratic differentials Q(d\, . . . , d m )) find those configu- 
rations of homologous saddle connections (or homologous closed geodesies), 
which are presented at any very flat surface S in the stratum. 

For quadratic differentials the notion of "homologous" saddle connections 
(closed geodesies) should be understood in terms of homology with local co- 
efficients, see [MaZo]. 

Topological, geometric, and dynamical properties of the strata 

Problem 13 (M. Kontsevich: Topology of strata; Sec. 3). 

Is it true that strata H(d\, . . . , d m ) and Q(qi, . . . , q n ) are K(n, l)-spaces 
(i.e. their universal covers are contractible)? 

Problem 14 (Compactification of moduli spaces; Sec. 6.3 and 6.4). 

Describe natural compactifications of the moduli spaces of Abelian differ- 
entials Tig and of the moduli spaces of meromorphic quadratic differentials 
with at most simple poles Q g . Describe natural compactifications of corre- 
sponding strata TL(d\, . . . , d m ) and Q(qi, ■ . ■ , q n ). 
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Problem 15 (Dynamical Hodge decomposition; Sec. 4-3, 5.7 and Appendix B). 

Study properties of distributions of Lagrangian subspaces in H 1 (S; R) de- 
fined by the Teichmuller geodesic flow, in particular, their continuity. Is there 
any topological or geometric way to define them? 

Problem 16 (Lyapunov exponents; Sec. 5.7 and Appendix B). 

Study individual Lyapunov exponents of the Teichmuller geodesic flow 

- for all known SL(2; R)-invariant subvarieties; 

- for strata in large genera. 

Are they related to characteristic numbers of some natural bundles over 
appropriate compactifications of the strata? 

Some other open problems can be found in [HuMSdtZ]. 
A Ergodic Theorem 

We closely follow the presentation in Chapter 1 of [CFSin]. However, for the 
sake of brevity we do not consider flows; for the flows the theory is absolutely 
parallel. 

Ergodic Theorem 

Consider a manifold M" with a measure p,. We shall assume that the measure 
comes from a volume form on M n , and that the total volume (total measure) 
of M n is finite. We shall consider only measurable subsets of M™. 

Let T : M™ — * M n be a smooth map. We do not assume that T is a 
bijection unless it is explicitly specified. We say that T preserves measure 
/j, is for any subset U C M n measure /i(T _1 C/) of the prcimage coincides 
with measure p(U) of the set. For example the double cover T : S 1 — > S 1 of 
the circle S 1 = R/Z over itself defined as T : x i— > 2x mod 1 preserves the 
Lebesgue measure on S 1 . In this section we consider only measure preserving 
maps. 

We say that some property is valid for almost all points of M n if it is valid 
for a subset U C M n of complete measure p(U) = p(M n ). 

A subset U C M n is invariant under the map T if the preimage T~ 1 U 
coincides with U . Thus, a notion of an invariant subset is well-defined even 
when T is not a one-to-one map. The measure-preserving map T : M n — > M n 
is ergodic with respect to \i if any invariant subset has measure or 1. The 
measure-preserving map T : M n — > M n is uniquely ergodic with respect to fi 
if there is no other invariant probability measure. 

Note that if T has a fixed point or, more generally, a periodic orbit (that is 
T k (xo) = Xq for some xq £ M n and some k € N), one can consider an invariant 
probability measure concentrated at the points of the orbit. Thus, such map T 
cannot be uniquely ergodic with respect to any Lebesgue equivalent measure. 
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Now we are ready to formulate the keystone theorem. Consider an inte- 
grablc function / on M n and some point x. Consider an orbit of T of length n 
starting at x. Let us evaluate the values of / at the points of the orbit, and let 

us calculate the "mean value" — (^f(x) + f(Tx) + ■ ■ ■ + f(T n ~ 1 x)j with respect 
to the discrete time k of our dynamical system . . . , T k ~ 1 x, T k x, T k+1 x, .... 

Ergodic Theorem. Let T : M n — > M n preserve a finite measure [i on M n . 
Then for any integrable function f on M n and for almost all point x G M n 
there exists the time mean: there exist the following limit: 



The function f(x) is integrable and invariant under the map f(Tx) = f{x). 
In particular, if T is ergodic, f is constant almost everywhere. Moreover, 



First Return Map 

The following theorem allows to construct numerous induced dynamical sys- 
tems which are closely related to the initial one. 

Theorem (Poincare Recurrence Theorem). For any subset U C M n 
of positive measure and for almost any starting point x G U the trajectory 
x, Tx, . . . eventually returns to U , i.e. there is some n > 1 such that T n x £ U . 

The minimal n = n{x) £ N as above is called the first return time. Accord- 
ing to Poincare Recurrence Theorem integer-valued function n(x) is defined 
almost everywhere in U. Consider the first return map T\jj : U — > U defined 
as T\jj : x i— > T n ( x 'x, where x G U. In other words, the map T\u maps a point 
x G U to the point where trajectory Tx, T 2 x, . . . first meets U. 

Lemma. For any subset U C M n of positive measure the first return map 
T\u : U — > U preserves measure /i restricted to U. IfT: M n — > M n is ergodic 
than T\jj : U — > U is also ergodic. 

The first return time induced by an ergodic map T has the following 
geometric property. 

Kac Lemma. For an ergodic diffcomorphism T : M n — > M n and for any 
subset U G M n of positive measure the mean value of the first return time 
equals to the volume of entire space: 





/ n{x)dn = u(M n ) 
Ju 
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Ergodic Theory: more serious reading. There are numerous nice books on 
ergodic theory. I can recommend a classical textbook of I. P. Cornfcld, 
S. V. Fomin and Ya. G. Sinai [CFSin] and a recent survey of B. Hassclblatt 
and A. Katok [HaKat]. 

B Multiplicative Ergodic Theorem 

In this section we discuss multiplicative ergodic theorem and the notion of 
Lyapunov exponents and then we present some basic facts concerning Lya- 
punov exponents. As an alternative elementary introduction to this subject 
we can recommend beautiful lectures of D. Rucllc [Ru]. A comprehensive in- 
formation representing the state-of-the-art in this subject can be found in the 
very recent survey [BP 2]. 

B.l A Crash Course of Linear Algebra 

Consider a linear transformation A : R™ — > R™ represented by a matrix A e 
SL(n, R). Assume that A has n distinct eigenvalues e Al , . . . , e A ™; let v l7 . . . , v„ 
be the corresponding eigenvectors. Note that since det A = 1 we get \j + ■ ■ ■ + 
A„ = 0; in particular, Ai > and A rl < 0. 

Consider now a linear transformation A N : R" — > R™, where TV is a very 
big positive integer. For almost all vectors v £ R" the linear transformation 
A N acts roughly as follows: it takes the projection v pro j of v to the line 
Vi = V^ec(vi) spanned by the top eigenvector vi and then expands it with a 
factor e 1 . So, roughly, A N smashes the whole space to the straight line Vi 
and then stretches this straight line with an enormous coefficient of expansion 
e^ 1 . (Speaking about projection to Vi = Vec(vi) we mean a projection 
along the hyperplanc spanned by the remaining eigenvectors V2, • • ■ , v„.) 

To be more precise, we have to note that the image of v would not have 
exactly the direction of vi. A better approximation of A N (v) would give us 
a vector in a two-dimensional subspace V2 = Vec(yi,V2) spanned by the two 
top eigenvectors Vi,V2. When A2 > the direction of A N (v) would be very 
close to the direction of Vi though the endpoint of A N (v) might be at a 
distance of order e NX2 from Vi, which is very large. However, this distance is 
small in comparison with the length of A N (v) which is of order e NXl ^ e NX2 . 

Further terms of approximation give us subspaces Vj spanned by the top 
j eigenvectors, Vj = Vec{v\, . . . , vj). Note that starting with some k < n the 
eigenvalues e Xk become strictly smaller than one. This means that the image 
A N (v) of any fixed vector v gets exponentially close to the subspace Vfc-i- 

Going into details we have to admit that vectors v from a subspace of 
measure zero in the set of directions expose different behavior. Namely, vectors 
v from the hyperplane £2 spanned by V2,. . . , v n have smaller coefficient of 
distortion than the generic ones. From this point of view vectors from linear 
subspaces Cj = Vec(\~j, Vj+i, . . . , n) expose more and more exotic behavior; 
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in particular, all vectors from the subspace Ck = Vec(vfc, Vj+i, . . . , v n ) get 
exponentially contracted. 

B.2 Multiplicative Ergodic Theorem for a Linear Map on the 
Torus 

Consider a linear transformation A : W l — > M. n this time represented by an 
integer matrix A £ SL(n, 1) as above. Consider the induced map F : T" — > T" 
of the torus T" = R"/Z™. This map preserves the natural linear measure on 
the torus. Consider the JV-th iterate F N of the map F, where N is a very 
large number. 

Note, that the differential of F in the natural coordinates on the torus is 
represented by the constant matrix D Xo F = A for any xq £ T". Note also 
that the differential of the iV-th iterate of F is represented as a product of 
N differentials of F along the trajectory x , F(x ), F(F(x )) ■ •• , F N ~ 1 (xo) of 
x : 

D X (F N ) = D FN - 1(XQ) F o • • • o D F(xo) F o D X0 F (35) 

Hence, in "linear" coordinates we get D X {F N ) = A N . Thus, the results of 
the previous section are literarily applicable to the local analysis of the map 
F N , where now vector v should be interpreted as a tangent vector to the 
torus T™. In particular, these results have the following interpretation. If we 
consider the trajectory xo, F(xq), . . . , F N (xo) of the initial point and the tra- 
jectory x,F(x), ... , F N (x) of a point x obtained from xo by a very small 
deformation in direction v, then for most of the vectors v trajectories would 
deviate exponentially fast one from the other; while for some special vectors 
they would approach each other exponentially fast. 

Namely, we get a distribution of linear subspaccs Ck in the tangent space to 
the torus such that deforming the starting point of trajectory in any direction 
in Ck we get two exponentially converging trajectories. The subspace Ck = 
Vec(yk, Vfc+i, . . . , v n ) is spanned by the eigenvectors of the matrix A having 
eigenvalues, which are smaller than one. This distribution is integrable; it 
defines a so-called stable foliation. 

There is also a complementary unstable foliation corresponding to the dis- 
tribution of subspaces Vk-i = Vec(vi, . . . , Vfc_i) spanned by the eigenvectors 
of the matrix A having eigenvalues, which are greater than one. Passing from 
the map F to the map F _1 the stable and unstable foliations change the roles: 
stable foliation of F becomes unstable for F" 1 and vice versa. 

When matrix A has an eigenvalue (or several eigenvalues) equal to ±1, we 
get also a neutral foliation corresponding to the distribution spanned by the 
corresponding eigenvectors. 

Exercise. Evaluate the limit 

io g |l^~)HI 

AT^oo N 
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for a tangent vector v G T xa T" having a generic direction. What are the pos- 
sible values of this limit for any tangent vector v e T xa T n 7 Show that vectors 
leading to different values of this limit are organized into a flag of subspaces 
£i Z) £2 Z) ■ ■ ■ Z) £ n , where we assume that all eigenvalues of the matrix A are 
positive and distinct. How would this flag change if some eigenvalues would 
have multiplicities? Would we have a flag of subspaces defined by different 
values of the limit above for the most general matrix A <G SL(n,l) (which 
may have Jordan blocks, complex eigenvalues, multiplicities, ...)? 

B.3 Multiplicative Ergodic Theorem 

Consider now a smooth measure-preserving map F : M™ — * M n on a manifold 
M n . We consider the case when that the total measure of M n is finite, and 
when the map F is ergodic with respect to this measure. 

Consider some generic point xq ■ Let us study, whether we have convergence 
of the limit (36) for tangent vectors v £ T Xo M n in this more general situation. 
We can always trivialize the tangent bundle to M n on an open subset of full 
measure. Using this trivialization we can reduce our problem to the study 
of product of matrices (35). This study is now much more difficult than in 
the previous case since the matrices D F k^ XQ ^F are not constant anymore. The 
following multiplicative ergodic theorem formulated for general mappings of 
general manifolds mimics the simplest situation with a linear map on the 
torus. 

Theorem (Oseledets). Let a smooth map F : M n — > M n be ergodic with 
respect to a finite measure. Then, there exists a collection of numbers 

Ai > A 2 > • ■ • > A fe , 

such that for almost any point x € M there is an equivariant filtration 

E" ~ T x M n = A D £ 2 D • • • D L k D £ k+1 = {0} 

in the fiber T x M n of the tangent bundle at x with the following property. For 
every v 6 £j - Cj+ii 3 = lj • • • > k, one has 

lim llog||pF JV ) :c (v)||=A j 

iV^+oo iv 

The multiplicative ergodic theorem was proved by V. L Oseledets [02]; 
a similar statement for products of random matrices was proved earlier by 
H. Furstenber [Fur]. 

Multiplicative ergodic theorem has several natural generalizations. The 
theorem essentially describes the behavior of products (35) of matrices along 
trajectories of the map F. Actually, matrices D X F are not distinguished by 
any special property. One can consider any matrix-valued function A : M n —> 
GL(m,M.) and study the products of matrices 
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A(F N - 1 (x)) A(F(x))-A(x) 

along trajectories x,F(x), . . . ,F N ^ 1 (x) of F. A statement completely anal- 
ogous to the above Theorem is valid in this more general case provided 
the matrix-valued function A{x) satisfy some very moderate requirements. 
Namely, we do not assume that A(x) is continuous or even bounded. The 
only requirement is that 

/ iog + p^)!!^ < +00, 

where log , (y) = max(log(y), 0). When this condition is satisfied one says that 
A(x) defines an integrable cocycle. The numbers Ai > • ■ • > Afc are called the 
Lyapunov exponents of the corresponding cocycle. 

Exercise. Formulate a "continuous-time" version of multiplicative ergodic the- 
orem when instead of a map F : M n — > M n we have a flow Ft which is ergodic 
with respect to a finite measure on M" . Show that under the natural normal- 
ization the corresponding Lyapunov exponents coincide with the Lyapunov 
exponents of the map F\ obtained as an action of the flow at the time t — 1. 

Consider a vector bundle over M"; suppose that the vector bundle is 
endowed with a flat connection. Formulate a version of multiplicative ergodic 
Theorem for the natural action of the flow on such vector bundle. 

Note that in the latter case the Lyapunov exponents are responsible for the 
"mean holonomy" of the fiber along the flow. Namely, we take a fiber of the 
vector bundle and transport it along a very long piece of trajectory of the flow. 
When the trajectory comes close to the starting point we identify the fibers 
using the flat connection and we study the resulting linear transformation of 
the fiber. 

Note that the choice of a norm in the fibers V x is in a sense irrelevant. 
Consider two norms || || and || ||' and let 

c(x) = min ||v||' C(x) = max ||v||'. 
IMI=i IMI=i 

If 

max (| log(c(a;)|, | log(C(x)|) fi < +00, 

then neither the filtration Ck(x) nor the Lyapunov exponents Xk do not change 
when we replace the norm || || by the norm || ||'. In particular, when M is a 
compact manifold all nonsingular norms are equivalent. 

In general, even for smooth maps F : M — > M (flows F t ) the subspaces 
defined by the terms Ck(x) C V x of the filtration do not change continuously 
with respect to a deformation of the base point x. However, these subspaces 
behave nicely for maps (flows) which have strong hyperbolic behavior (see [Po] 
for a short introduction; a recent quite accessible textbook [BP1] and a sur- 
vey [BP2] describing the contemporary status of Pesin theory). 



1 
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Currently there are no general methods of computation of Lyapunov expo- 
nents other than numerically. There are some particular situations, say, when 
the vector bundle has a one-dimensional equivariant subspacc, or when F t is 
a homogeneous flow on a homogeneous space; in these rather special cases the 
corresponding Lyapunov exponents can be computed explicitly. However, in 
general it is extremely difficult to obtain any nontrivial information (positiv- 
ity, simplicity of spectrum) about Lyapunov exponents. 
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